VOLUME AUGUST 1957 


Canadian 
Journal Physics 


Editor: DUCKWORTH 


Associate Editors: 


ELLIOTT, Atomic Energy Canada, Ltd., Chalk River 
FOSTER, McGill University 

HERZBERG, National Research Council Canada 

LEPRINCE-RINGUET, Ecole Polytechnique, Paris 

SARGENT, Queen’s University 

VOLKOFF, University British Columbia 

WATSON, University Toronto 

WOONTON, McGill University 


NUMBER 


Published THE NATIONAL RESEARCH COUNCIL 


OTTAWA 


CANADA 


CANADIAN JOURNAL PHYSICS 


(Formerly Section Canadian Journal Research) 


Under the authority the Chairman the Committee the Privy Council Scientific 
Industrial Research, the National Research Council issues THE CANADIAN JOURNAL 
and five other journals devoted the publication, English French, the 
results original scientific research. Matters general policy concerning these journals are 
the responsibility joint Editorial Board consisting of: members representing the National 
Research Council Canada; the Editors the Journals; and members representing the 
Royal Society Canada and four other scientific societies. 


EDITORIAL BOARD 


Representatives the National Research Council 


Miller, University Alberta Thomson, McGill University 
Thode, McMaster University Watson (Chairman), University Toronto 


Editors the Journals 


Bailey, University Toronto Elliott, Montreal Neurological Institute 
Cameron, Macdonald College Léo Marion, National Research Council 


Duckworth, McMaster University Murray, University Western Ontario 


Representatives Societies 


Bailey, University Toronto A.C. Elliott, Montreal Neurological Institute 
Royal Society Canada Canadian Physiological Society 

Cameron, Macdonald College Murray, University Western Ontario 
Royal Society Canadian Society Microbiologists 

Duckworth, McMaster University Thode, McMaster University 
Royal Society-of Canada Chemical Institute Canada 
Canadian Association Physicists Thorvaldson, University Saskatchewan 


Royal Society Canada 


officio 


Léo Marion (Editor-in-Chief), National Research Council 
Rosser, Vice-President (Administration and Awards), National Research Council 


Manuscripts for publication should submitted Dr. Duckworth, Editor, Canadian 
Journal Physics, Hamilton College, McMaster University, Hamilton, Ontario. 


(For instructions preparation copy, see Notes Contributors (inside back cover).) 

Proof, correspondence concerning proof, and orders for reprints should sent the Manager, 

Editorial Office (Research Journals), Division Administration and Awards, National 
Research Council, Ottawa Canada. 

Subscriptions, renewals, requests for single back numbers, and all remittances should sent 
Division Administration and Awards, National Research Council, Ottawa Canada. 
Remittances should made payable the Receiver General Canada, credit National 
Research Council. 


The journals published, frequency publication, and prices are: 


Canadian Journal Biochemistry and Physiology Monthly $3.00 year 
Canadian Journal Botany Bimonthly $4.00 year 
Canadian Journal Chemistry Monthly $5.00 year 
Canadian Journal Microbiology Bimonthly $3.00 year 
Canadian Journal Physics Monthly $4.00 year 
Canadian Journal Zoology Bimonthly $3.00 year 


The price single numbers all journals cents. 


\ 
j 2 
: 
j 


CORRECTIO 
Canadian Journal Physics, Vol. 35. 


.20327975 should read —.20327975. 
.25141195 should read .25141495. 


should read 


345. the section 9/2, the numbers printed the 7/2 column should 
moved the 9/2 column. 


N 
| 
Ke 


Canadian Journal Physics 


Issued THE NATIONAL RESEARCH COUNCIL CANADA 


VOLUME AUGUST 1957 NUMBER 


NEW TECHNIQUE FOR THE STUDY 
SCATTER PROPAGATION THE TROPOSPHERE! 


ABSTRACT 


The power spectrum the fluctuations received signal strength near- 
optical U.H.F. circuit has been measured. The sidebands associated with 
these fluctuations can overlap the information-carrying sidebands com- 
munication system. When this happens, these sidebands must taken into 
account determining the signal-to-noise ratio the system. other words, 
the fluctuations then have the characteristics noise, and therefore they are 
called propagation noise paper. Experiments carrier frequency 
500 Mc. have shown that the propagation noise power density usually varies 
with sideband frequency the carrier) for the range 
0.1 Departures from this law have been observed the regions near 
0.1 and The measurement the power spectrum directly offers 
several advantages over the conventional signal strength recording method, and 
these are discussed herein. 


INTRODUCTION 


During the last few years, radio waves U.H.F. have been used consider- 
ably broadcasting and communication systems. result, wave pro- 
pagation U.H.F. through the troposphere has been subject intensive 
studies (Proceedings the Institute Radio Engineers 1955). These studies 
have shown the feasibility long-range beyond-the-horizon communication 
frequencies very much higher than the maximum frequencies usable for 
ionospheric propagation. This long-range propagation has been called ‘scatter 
propagation’ the belief that due scattering from irregularities 
the refractive index structure the troposphere. However, general agreement 
the exact mechanism has not yet been achieved. 

The irregularities which are associated with this long-distance scatter 
propagation affect propagation over short line-of-sight paths also. some 
ways the short paths are more suitable for studies scattering radio waves 
atmospheric inhomogeneities than are the long paths. For example, the 
use short paths and moderately narrow beam antennas involves scattering 
volumes near the surface the earth where atmospheric studies are most 
easily made. addition, the atmospheric conditions are more likely 
uniform over the shorter path. experimental method that especially 
applicable the short-path experiments presented this paper. 


received original form November 13, 1956, and, revised, May 1957. 

Contribution from Defence Research Telecommunications Establishment, Defence Re- 
search Board, Ottawa. Research conducted under Defence Research Board Project D48-38- 
03-24. 
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FACTORS AFFECTING U.H.F. PROPAGATION 


There are five factors which have appreciable effect upon U.H.F. signals 
propagated through the troposphere: diffraction occurs around obstacles 
the path (Bullington 1947, 1955; Crysdale et. 1956); reflections occur 
ground water surfaces (Kerr 1951); refraction takes place because the 
vertical gradient refractive index that normally exists the lower atmos- 
phere (Kerr 1951); scattering reflection may take place from large-scale 
structure the troposphere, such clouds elevated layers (Northover 
1955; Starkey 1956); finally, scattering may caused also small-scale 
turbulent structure the troposphere (Sutton 1953; Booker and Gordon 
1950; Gordon 1955; Villars and Weisskopf 1954, 1955; Fannin 1956; Muchmore 
and Wheelon 1955). 

Other phenomena, such reflections from ionospheric layers, meteor trails, 
and aurora, which influence radio wave propagation lower frequencies, 
can neglected the U.H.F. region. Absorption can neglected also, 
except the highest microwave frequencies. Some miscellaneous effects, 
such reflections from aircraft, usually can identified and accounted for. 

Qualitatively, least, these factors can divided into two classes: first, 
those which produce slow changes signal strength and, second, those which 
produce rapid changes. slow change one whose effect radio system 
can described terms signal strength alone. rapid change, the 
other hand, one which the rate change signal strength well 
the instantaneous signal amplitude must taken into account. This done 
considering the sidebands associated with the fluctuations. these side- 
bands overlap the sidebands that carry the information, then the fluctuation 
signal strength termed rapid one. 

can seen that this classification change signal strength fast 
slow depends the particular system being considered. Modern equipment 
makes possible use sidebands removed only few cycles per second 
from the carrier. Thus, for the sake being definite, fluctuations the 
order lower might considered slow, and those above 
c.p.s. fast. any rate, information the behavior these sidebands 
can used determine whether noise due propagation equipment 
limits the minimum frequency shift usable. 

When the fluctuations signal strength are fast (for particular system), 
the information signals lowest frequency must compete with the fluctua- 
tions. Clearly then, the variations take the characteristics noise; hence 
the fluctuations are called propagation noise this paper. line with this 
outlook, experimental methods common noise studies have been used. 
The power spectrum the fluctuations has been determined experimentally 
means band-pass filters. This technique well known noise studies, 
for example, used van der Ziel (1954). has proved quite applic- 
able the present problem, and offers several advantages over the con- 
ventional signal strength analysis. These will discussed the end this 
paper. 

Consideration the above-mentioned factors radio wave propagation 


A 
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leads the conclusion that the first three produce slow changes signal 
strength, while general the other two produce fast changes. Since there are 
experimental difficulties the use band-pass filters the frequencies 
involved slow signal changes, the method here presented best suited 
the last two factors. short, the direct determination the power spectrum 
propagation noise should permit detailed study those features the 
radio signal received over short path which are due scattering from 
atmospheric irregularities. 


EXPERIMENTAL RESULTS PROPAGATION NOISE 

Experiments propagation noise have been carried out the Ottawa 
valley the Radio Physics Laboratory since January, 1955. Three experi- 
mental paths have been used. The transmitters were located White Lake, 
Meath, and Foymount, distances 31, 60, and miles respectively from 
the receiver site R.P.L. Path profiles four-thirds earth basis are shown 
Fig. Other paths were used for short tests. The equipment used in- 
dicated the block diagram, Fig. principle, the equipment very simple, 
but was found, discussed below, that the equipment had very 
high quality. Standard communications equipment, operating 500 


Fic. Path profiles for the three paths used for most the experiments, based 
four-thirds earth radius allow for normal refraction. 


30 
600 


| 
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RECEIVER RECORDER 


Fic. Block diagram the equipment used the experiments propagation noise. 
For continuous recordings several band-pass filters were used parallel. 


VARIABLE 
BAND-PASS 
FILTER 


TRANSMITTER 


was modified order obtain the lower noise levels needed. the experi- 
ments date, 5-ft. parabolic antennas were used, with gains about db. 
relative isotropic radiators. These were mounted 50-ft. towers, although 
other heights 120 ft. were used for short tests. 

The audio signal was taken from the second detector the receiver. 
This was analyzed the use Krohn-Hite model 330-A band-pass filters 
which have upper and lower cutoff frequencies independently variable be- 
tween 0.02 c.p.s. and 2000 c.p.s., and have minimum usable bandwidth 
one octave. For continuous recordings, several filters parallel were used, 
but short determinations power spectrum often were made using one filter 
with various pass bands turn. The outputs from the filters were recorded 
directly. using Edin pen recorders, run slow chart speeds. the use 
well-known statistical relations for Gaussian noise (Rice 1945), the r.m.s. 
value each filter output was determined from the records. The receiving 
and recording equipments were calibrated either using the front end noise, 
taking advantage the known noise figure, the use crystal diode 
modulator. 

Most the measurements made indicated that the noise the output 
the second detector has power spectrum given approximately 


(1) wf) 


Here w(f) the power density, power per unit bandwidth, and the 
receiver noise power density, the latter being independent frequency. 
Also, independent the dimensions are discussed below. Fig. 
gives one experimental result. For the results shown Fig. each separate 
measurement was made turn with one filter, the pass band always being 
one octave wide. emphasize the fact that filter finite bandwidth was 
used, each measurement depicted horizontal line rather than 
point, where the length the line indicates the bandwidth the filter. Two 
determinations the power spectrum were made, about hours apart, 
with each determination taking about minutes; the results are shown sepa- 
rately the graph. The agreement between the two runs demonstrates the 
stationary character the fluctuations during this time. addition the 
experimental results, Fig. shows curve representing equation (1) with 
380. The constant value w(f) very high frequencies, 


means simply that receiver noise predominates. lower frequencies, where 
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SPECTRUM 1300 HRS, 


wath 


TOTAL NOISE PER UNIT BANDWIDTH 


FREQUENCY, CPS 


Fic. Propagation noise power spectrum mile path, Feb. 16, 1955. The short 
horizontal lines indicate the average power density the filter pass bands, decibels above 
the receiver noise power density. The smooth curve represents the power spectrum given 
equation (1) with 


receiver noise can neglected relation propagation noise, inverse 
square relation holds, 


(3) 


Here has been written bring out the dimensions explicitly. 
Continuous records showed that the exponent equation (3) was not 
always exactly equal This evident from the sample continuous re- 
cordings Fig. Three filters, separated frequency factors five, 


= 
| 
= 
a 
A-0.018 TO 0.08 CRS, 
4 


TIME, HOURS AND DAYS 


Fic. continuous record propagation noise mile path. The three curves 
represent the noise power densities the outputs three band-pass filters indicated, 
decibels above receiver noise. 
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were used here simultaneously. Therefore, the logarithmic plot, the in- 
verse square law were hold the three curves should separated db. 
apparent that the inverse square law fairly well obeyed the higher 
frequencies, but not valid the low frequencies. the low frequencies, 
slightly smaller exponent indicated. Other results have indicated stronger 
dependence than inverse square above c.p.s. Nevertheless, the inverse 
square law very good representation results the region immediately 


above 0.1 


IV. DISCUSSION 


Over short path the received signal quite strong, and the fluctuations 
are small part the total signal. mile path, for example, the 
value the fluctuating part was only few per cent, most, the mean 
signal. Therefore, some means suppressing the steady component must 
used the fluctuations are studied the short path. Differential 
d-c. recorders have been used others but perhaps easier experimentally 
work with filters instead. either case, equipment noise levels, 
ticular transmitter hum and receiver front end noise, determine the minimum 
propagation noise levels that can studied. was found that dynamic 
range least db. was needed for the study propagation noise the 
present paths. 

The experimental observation that the power spectrum propagation 
noise follows inverse square law for wide range frequencies has direct 
bearing the choice equipment. Suppose that signal strength recording 
made with fast differential recorder, for the purpose studying fast 
fading. The recorder sensitivity determined the low frequencies, since 
they occur with the larger amplitudes shown the present results. Since 
the higher frequency components occur with diminishing amplitudes, they 
will not apparent the recording, owing its finite resolution. Such 
recording thus acts filter, which discriminates against the high frequencies 
because their small amplitude. The use band-pass filters, the present 
experiments, permits the use higher gain for the higher frequencies. The 
dynamic range the recording instrument need not then limit the range 
variation the variables interest. 

should noted also that, compared with conventional methods, the 
use narrow band-pass filters reduces the amount computation needed 
the analysis the results. The determination either the power spectrum 
the autocorrelation function from signal strength recording tedious 
task. 

The experimental result here presented indicates predominance low- 
frequency components the spectrum propagation noise. However, for 
physical reasons, obvious that the inverse square law cannot hold down 
fading rate the signal. The fading rate may defined the number 
crossings the mean signal level, which given Rice (1945) 
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The lower integral diverges the lower limit and the upper integral diverges 
the upper limit. Physically, the divergence the lower limit implies 
infinite amount power the sidebands. This can corrected adding 
constant term the denominator, i.e. assuming the power spectrum 


(5 ) = wof 2 (i? +f?) 
or 
(5d) 


where some characteristic frequency for the system. this form, 
represents constant power density frequencies well below This 
power spectrum implies exponential autocorrelation function. 

The difficulty the upper limit remains. related the cusp the 
origin the exponential autocorrelation function. This has been discussed 
others, e.g. Rice (1945) and Muchmore and Wheelon (1955). Physically, 
implies completely discontinuous process and, practice, parasitic 
effects always enter the problem introduce high-frequency cutoff. Mathe- 
matically, this may taken care stronger dependence frequency, 
e.g. power spectrum the form 


or, equivalently, 


This form the power spectrum equivalent the rounded exponential 
autocorrelation function proposed Muchmore and Wheelon (1955). 
Incidentally, the fading rate then geometric mean between the two charac- 
teristic frequencies. 

might added here that deficiencies the experimental equipment, 
e.g. the finite resolution signal strength record, introduce the upper cutoff 
frequency F,, then the observed fading rate partly determined the 
equipment, and not parameter the received signal alone. 

conclusion, can said that the simple inverse square law describes 
most the present results, obtained variety paths and weather 
conditions. Nevertheless some the experiments have indicated upper and 
being the order 0.01 0.1 with the 500 carrier frequency. 
While the form the power spectrum seems independent such things 
path length and weather conditions, seems obvious that the three sugges- 
ted parameters wo, and cannot be. These will discussed the 
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ROTATIONAL TEMPERATURES MEASURED AURORAE 
CHURCHILL, MANITOBA! 


ABSTRACT 


series auroral spectra showing the band have been obtained 
Churchill, Manitoba. Rotational temperatures ranging from 230° 480° 
have been deduced. These are agreement with measurements made Saska- 
toon, Saskatchewan, and indicate that there latitude effect. vertical 
temperature gradient some K./km. suggested the observations. 


INTRODUCTION 


One the methods used determine temperatures the upper atmosphere 
auroral heights measure the intensity distribution the rotational 
fine structure the first negative system the aurora. The validity 
this method has been discussed and demonstrated for the case electron 
excitation the band Vallance Jones and Harrison (1955) and Shepherd 
and Hunten (1955). Their results for Saskatoon indicate temperatures ranging 
from approximately 170° 450° K., depending the type and height 
the aurorae. The results Shepherd and Hunten, particular, point 
vertical temperature gradient some K./km. agreement with the 
measurements made rocket flights (The Rocket Panel 1952). Petrie (1953), 
also Saskatoon, but during period greater auroral activity, measured 
temperatures ranging from 450° 700° These were likely derived from 
spectra obtained near the summit visible aurorae. the other hand, 
Vegard (1951), the basis large number measurements, concluded 
that constant temperature 230°+20° characteristic all types and 
heights aurorae. Furthermore, this result was obtained both Oslo and 
latitudes 60° and 69.7° respectively. 

Vegard (1956) has suggested that the higher temperatures measured 
Saskatoon were result proton excitation the band. This possibility 
based the Solar Ray Bundle theory (Vegard 1956) which predicts that 
the proton flux should greater low latitude and high altitude aurorae. 
Recent observations (Montalbetti and Vallance Jones 1957) confirm that 
proton flux tends associated with the southern extent auroral 
display. 

Since Churchill situated the center the auroral zone, seemed 
desirable make temperature measurements this station and compare 
them with the results obtained Saskatoon. The intensity any aurora 
which occurred Saskatoon during the interval that exposures were being 
made was noted, since this provided indication the position the 
proton flux. 


'Manuscript received May 15, 1957. 
Contribution from Radio Physics Laboratory, Defence Research Board, Ottawa, Ontario. 
address: Physics Department, University Saskatchewan, Saskatoon, Saskat- 
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EXPERIMENTAL 

high-dispersion spectrograph consisting essentially 8X8 in. grating 
and camera having relative aperture was used. step entrance 
slit was employed that certain spectral features could used determine 
the the emulsion. The aurora was imaged the 
slit means 14-in. focal length quartz lens and movable mirror. The 
band N.* was photographed the third order with dispersion 
about A/mm. 

During the months August November, 1954, and November, 1956, 
number satisfactory plates were obtained. attempt was made differ- 
entiate between forms aurorae which occurred during any one exposure 
interval, but was observed that most the displays were either strong 
weak. practically all cases satisfactory plate densities resulted from exposure 
times hours, the length the exposure depending the intensity 
the visible aurorae. 

microphotometer tracing the band obtained the night 
September 30, 1956, reproduced Fig. clear that the band 


Fic. Microphotometer tracing first negative band. 


sufficiently resolved determine the intensity distribution the rotational 
fine structure, lines characterized rotational quantum numbers greater 
than six are Near the band head the overlapping the individual 
lines prevents accurate determination their relative intensities. 
independent methods were used measure the relative intensities the 
individual lines. One was take the 2:1 intensity ratio between 
alternate lines, was done Vallance Jones and Harrison (1955); and the 


other was utilize the spectrum obtained the stepped slit. purposes 


26 
\ 
~ bed 
~~. 
~ 
V 
=> 
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MONTALBETTI: 


calibration, the 3884 band may considered continuum, and hence 
microphotometer tracing taken along the length the image the stepped 
slit gives the relation between plate density and relative intensity. Fig. 
example the calibration curve obtained the two methods for the 
plate shown Fig. 


INTENSITY (ARBIT. UNITS) 


from ratio 
between lines 


from stepped 
slit 


PLATE UNITS) 


Fic. Characteristic curve emulsion plate shown Fig. 


The temperature was determined two well-known methods: 

Method rotational quantum number which corresponds the 
maximum intensity the branch determined. This quantity related 
the absolute temperature the following equation: 


Method Since the relative intensity the components the 
branch varies with the rotational quantum number according the 


relation 


plot A’) vs. vields straight line. The slope this 
line related the absolute temperature 


that was well exposed. almost every case straight line resulted. 
two instances curved lines were obtained which could represented essenti- 
ally two straight lines. The temperatures calculated trom the initial and 
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Fic. Plots against for auroral plates described Table 


final slopes correspond rather conservative upper and lower limits for the 
temperature range shown Vallance Jones and Harrison (1955). 
worth noting that the exposure time for these two plates was the order 
hours, whereas for all the other plates the exposure time ranged from 
hours. number plates were not sufficiently well exposed measure 
temperatures the slope method but was still possible determine the 
position the maximum the intensity distribution, and hence the tempera- 
ture, from them. 
RESULTS AND DISCUSSION 

Table lists the temperatures calculated both methods. and 
correspond the temperatures deduced method and method 2respectively. 
The temperatures listed the bottom the table were calculated method 
only, since these plates were not sufficiently well exposed make possible 
accurate intensity measurements the individual rotational lines. These 
temperatures are probably reliable, since there good agreement between the 
two methods seen from the upper portion the table, and view the 
fact that only very few cases were straight line plots not obtained. 
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TABLE 
Aurorae Churchill Temperature, 
Saskatoon Curve No. 
1954 
Aug. RA, 2-3 245 
RA, PS, 285 275 None 
270 270 None 
HA, PS, 2-3 250 220 
Faint arcs 350 320 
Auroral clouds 290 220 
PS, faint arcs 330 315 
Clouds, arcs 300 315 
1956 
Oct. RB, DS, 1-2 340-660 360 
Nov. RB, HA, 1-2 440 410 
RB, RA, 1-2 480 430 
RB, DS, HA, 1-2 310-490 None 
1954 
Sept. Auroral clouds 315 
Faint arcs 315 
RA, PS, 2-3 230 
Faint arcs 315 


Mean temperature weak aurorae (Intensity 0-2) 
250+23° 


Mean temperature strong aurorae (Intensity 


Observations the occurrence the emissions Saskatoon and 
Churchill (latitude 52.1° and 58.8° respectively) indicated that the proton 
flux tends associated with the southern extent auroral display 
(Montalbetti and Vallance Jones 1957). Similar observations have been 
reported Vegard (1952). Thus the occurrence aurorae over Saskatoon was 
taken evidence that the proton flux was south Churchill. the second 
last column Table the intensity the aurora which occurred over Saska- 
toon during the time that exposure was being made Churchill listed. 
proton excitation the band important, likely that higher 
temperatures would measured Churchill during times weak 
aurora over the southern station. Examination Table shows that 
correlation exists between the temperature measured Churchill and the 
intensity aurora over Saskatoon. Furthermore, these temperatures range 
from approximately 230° 480° K., agreement with the results obtained 
the Saskatoon group. These results suggest that either the excitation 
the band protons unimportant the protons excite the nitrogen 
molecules such manner not affect the thermal equilibrium. 

The intensity the aurora Churchill which occurred during the exposure 
interval listed the third column Table Averaging the temperature 
variation for weak aurora and for strong aurora shows that the temperature 
some 85° greater for the weak auroral displays. According Stérmer 
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(1955) the average difference height between the two types displays 
approximately Thus vertical temperature gradient K./km. 
noted agreement with other observations (The Rocket Panel 1952; Shepherd 
and Hunten 1955). 

mentioned the previous section, the exposure times for the two cases 
when curved lines were obtained was the order hours. thus likely 
that the exposure was sufficient duration obtain spectra very weak 
aurorae superposed the spectra stronger aurorae. The upper temperature 
limit deduced from such spectra very likely corresponds very weak high 
altitude aurorae. 
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HIGH RESOLUTION RAMAN SPECTROSCOPY GASES 
ROTATIONAL SPECTRUM BUTATRIENE 


ABSTRACT 


The rotational Raman spectrum butatriene 
pressure 2cm. was photographed with ft. grating spectrograph. 
analysis this spectrum (based the symmetric top approximation) yields 
the rotational constant the two outer C=C 
bonds butatriene are assumed have the same length the bonds 
allene, namely. 1.309 found that the central C=C bond has length 
value which shorter than that the bonds ethylene 


and allene. 
INTRODUCTION 

The length single bonds known depend primarily 
adjacent bond types (for summary see Herzberg and Stoicheff 1955). Recently 
dependence bond length environment has also been found for 
carbon double bonds, and the present paper reports further results support 
such dependence. 

Evidence for double bond considerably shorter than that 
ethylene was found ketene from its microwave spectrum (Johnson and 
Strandberg 1952) and allene from its infrared and Raman spectra (Overend 
and Thompson 1953; Stoicheff 1955). The lengths and allene are 
respectively and compared with 1.344 for 
This difference approximately 0.03 bond length apparently 
due the fact that allene and ketene, the bonds are adjacent two 
single bonds and one double bond while ethylene the bond adjacent 
four single bonds. butatriene the central bond 
adjacent two double bonds, and the present analysis the rotational 
Raman spectrum butatriene shows that this bond even shorter than that 
allene and ketene. 

EXPERIMENTAL 

The apparatus described paper (Stoicheff 1954) with some modifications 
was used photograph the rotational Raman spectrum; four lamps replaced 
the original two, and cylindrical lens was placed front the photographic 
plate the spectrograph. The sample butatriene was prepared Merck 
and Co., Montreal. mass spectrometer analysis showed that least 95°; 


received May 1957. 

Contribution from the Division Pure Physics, National Research Council, Ottawa, 
Canada. 

Issued No. 4404. 

and Strandberg (1952) obtained the value for the sum the 
C=C and C=O distances ketene. one assumes that 1.162 found for 
the molecule (Herzberg and Herzberg 1953), then ketene. 

value 1.344 ethylene was calculated from the Raman value 
(Romanko, Feldman, Stansbury, and 1954) using and 
Z.HCH 120° obtained from the infrared spectra and (Gallaway and Barker 


1942). 
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was mass 52; the remainder was mainly mass and was probably 
butadiene. 

difficulty was encountered when filling the Raman tube with butatriene 
pressures higher than cm. Hg, the Rayleigh light became extremely intense 
within few minutes introducing the gas into the tube. This phenomenon 
was due the rapid polymerization the butatriene the impurities. The 
difficulty was overcome working with pressures 1-2 cm. Hg, renewing 
the sample every hour, and keeping the Raman tube and filter 
solution surrounding the tube temperature about Fortunately, 
butatriene with its three double bonds has extremely intense rotational 
Raman spectrum and was possible photograph the spectrum even 
this low pressure reasonable exposure times. 

Several photographs the rotational Raman spectrum were obtained 
the second order the ft. grating, two with the cylindrical lens placed 
front the photographic plate and two without this lens. The exposure 
times the former were and hours and the latter and hours. Fig. 
shows the Raman spectrum the Stokes side the exciting line, photographed 
hour. The rotational lines have spacing about 0.52 and are 
fairly well resolved. However, the spectrum falls within the exciting 
line, region which masked grating ghosts, and therefore many the 
rotational lines are superimposed intense background and some are 
blended with the grating ghosts. Second- and third-order iron lines excited 
iron hollow cathode (Stanley and Dieke 1955) were used wavelength 
standards. 

Two plates were evaluated, one photographed hour with the cylindrical 
lens, and the other hours without this lens. The wavenumber displace- 
ments agreed within and the average values are listed 
Table 


TABLE 
OBSERVED WAVENUMBER SHIFTS THE ROTATIONAL RAMAN LINES BUTATRIENE 
| 


“All values are averages the Stokes and anti-Stokes lines measured two plates. 
*Lines blended with grating ghosts. 


4 
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ROTATIONAL ANALYSIS 


The simple appearance the rotational spectrum, with its relatively sharp 
lines and narrow regular spacing, suggests that the rotational lines correspond 
transitions the type that is, the observed spectrum 
resembles that symmetric top. This, course, expected once 
assumed that the four carbon atoms butatriene are collinear. this 
case slightly asymmetric prolate top, the selection rule holds 
approximately and the most intense lines are due transitions 
Therefore, the displacements the Raman lines can expressed good 
approximation the equation 


where The observed line spacing was used deter- 
mine the rotational numbering and the rotational constants were evaluated 
the usual way plotting graph (Fig. against 
The values obtained are 


BUTATRIENE 


graph for the ground state butatriene represents lines blended 
with grating ghosts). 


The rotational constant related the moments inertia the molecule 


MOLECULAR STRUCTURE 
Although the symmetry the butatriene molecule has not yet been deter- 
mined experimentally, there doubt that the molecule planar and 
belongs the point group V,. Its structure and four structural parameters 
are represented Fig. The moment inertia about the axis 


and that about axis through the center and perpendicular the plane 
the molecule 


I, = +2my (rs? +h?)]. 


When assumed that the C—H distance 1.07+.01 and the HCH 


0.5260 
0.5240 
0.5230 
1000 2000 3000 4000 
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Structure and structural parameters butatriene. 
angle (the values found allene Stoicheff 1955), the above 


bonds are assumed equal length one finds value 
1.299 But since the central bond and the outer bonds 
have different adjacent bond types seems preferable assume 
that the bond lengths differ. this case the similarity the environment 
the outer carbon-carbon bonds that allene suggests that the length 
1.309+0.003 found for the bonds allene assumed for 


these bonds. This leads value 1.284+0.006 for the length the 
central double bond butatriene, value considerably shorter 
than allene ketene. 

The known lengths double bonds, obtained from spectro- 
scopic investigations, can summarized follows: 


1.344A 1.309 1.284 


All three values quoted here refer values, that is, they include zero point 
vibration and hence might differ from the equilibrium internuclear distances 


much 0.005 the three values, that for probably the 


most accurately known, and mentioned earlier (Section nearly equal 


values were found allene and ketene. The present value for 


not without objection even though the values and for ethylene 
used this calculation are thought the most accurate available. more 
reliable value can probably obtained study the pure rotational Raman 
spectra and under the conditions high resolution used 
the present research. 

short length for the bond has also been obtained from electron 
diffraction experiments with carbon suboxide Mackle 
and Sutton (1951) found value 1.28+0.03 and more recently Bastiansen 
(1951) has found 1.279 Yet spite the good agreement this value cannot 
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considered definitive since there still some doubt the linearity 
the molecule. 

the present paper, admittedly, several assumptions were necessary 
order obtain value for the central bond length butatriene, 
and would desirable determine the structural parameters independently 
any assumptions. the value for the isotopic molecule were 
known, could lead the independent evaluation the central and outer 
carbon-carbon bond lengths. However, the limited accuracy with which these 
rotational constants can presently obtained would not lead significant 
result. 


grateful Dr. Lossing for the mass spectrometer analysis the 
sample butatriene, and Dr. Herzberg for helpful discussions. 
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ULTRAVIOLET ABSORPTION SPECTRA HCN AND DCN 


HERZBERG AND INNEs? 


ABSTRACT 


The ultraviolet absorption spectra HCN and DCN have been investigated 
below 2000 Four band systems have been found which two, a—X and 
are described and analyzed the present paper. the long wavelength end 
the system both HCN and DCN consists simple progressions 
sharp bands whose fine structure has been measured. shorter wavelengths 
diffuseness (predissociation) sets gradually but somewhat different energies 
HCN and DCN. Both the vibrational and rotational analyses lead indepen- 
dently and unambiguously the conclusion that the upper state, unlike the 
ground state, the molecule non-linear. The geometrical parameters obtained 
from the rotational fine structure are ro(CN) 1.297 
and 125.0°. Two vibrational frequencies the state are 
The stretching frequency has not been found. 

The system has been clearly observed only for DCN. HCN 
apparently strongly predissociated that its presence among the diffuse a—X 
bands difficult establish. the state, the state, the molecule 
bent. The somewhat fragmentary rotational analy sis leads 
and 114.5° when assumed. Only one 
vibrational frequency, has been established the state 
DCN. The vibrational quantum number uncertain since isotope effect 
available establish the position the 

the rotational analysis, particularly the a—X system, the effects 
the asymmetry the molecule the upper state are clearly demonstrated 
the doubling the levels with and and the large contributions 
the asymmetry the terms From the sign the doubling 
species i.e. the electronic eigenfunctions are antisymmetric with respect 
the plane the molecule. appears probable that both transitions a—X 
and correspond forbidden transitions and the linear 
case. The non-linear structure the excited states briefly discussed terms 
electron configurations. 

number interesting features the predissociation the state are 
pointed out. shown that the dissociation products are all probability 


INTRODUCTION 

Hydrogen cyanide exhibits absorption throughout the visible and near 
ultraviolet regions down 2000 Price (1934) and Price and Walsh (1945) 
were the first investigate the absorption the vacuum ultraviolet and 
found extremely strong diffuse absorption bands below 1120 and somewhat 
weaker sharp absorption bands the region 1550 1350 Hilgendorff (1935) 
reported the observation another absorption band system longer wave- 
lengths (2000-1700 which very much weaker than band systems 
and requires absorbing path about cm-atm. several respects this 
HCN absorption appears analogous the near ultraviolet bands while the 
HCN bands near 1500 appear analogous the bands the same region. 
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The fine structure the near ultraviolet bands has recently been 
analyzed detail Ingold and King (1952, 1953) and Innes (1954). The 
analysis showed unambiguously that the upper state the molecule 
has non-linear structure and the dimensions the molecule this state 
were determined from the spectrum. view the similarity the electronic 
structure HCN and and view the fact that few rotational 
structures electronic spectra polyatomic molecules have been analyzed, 
appeared worth while attempt more detailed investigation the HCN 
and DCN bands including, where possible, analysis the rotational 
structure. the present paper, report the results the study the 
first weaker absorption region (2000-1600 where two band systems have 
been found (designated a—X and later paper will deal with the 
second, much stronger absorption region 1500 and below. The principal 
result these investigations that the HCN molecule non-linear 
least three excited electronic states. 


EXPERIMENTAL 


All HCN and DCN spectra were taken the fourth order the meter 
vacuum grating spectrograph described Brix and Herzberg (1954). The 
dispersion was 0.65 A/mm. and the resolving power actually obtained was 
150,000. The spectra were measured against iron lines second order. The 
accuracy non-blended lines estimated +0.05 

Two absorption tubes were used, one the other 100 cm. long. They 
were filled with HCN DCN pressures ranging from 0.01 240 mm. 
Unlike the case because the low vapor pressure HCN, was 
not possible obtain spectra low temperature. 

For the earlier exposures with the first absorption tube, the continuous 
background for the absorption spectra was Lyman flash discharge tube 
viewed end-on described Brix and Herzberg (1954). Some the later 
exposures above 1500 with the second absorption tube were taken with the 
superior continuum discharge described Wilkinson and Tanaka 
(1955). 

The discharge tube was separated from the absorption tube spherical 
LiF lens whose focal length was chosen that light the desired wavelength 
range was focused the grating. LiF prism with one cylindrical face 
(Curry and Herzberg 1934) was placed between the absorption tube and the 
spectrograph and focused the desired spectral region the slit the spectro- 
graph. 

OBSERVED SPECTRA 


the region from 2000 1600 there are several fairly long progressions 
simple bands both HCN and DCN. Fig. shows sections the 
spectrum for both isotopic molecules. the long wavelength end the bands 
are quite sharp, each having two heads and well-developed fine structure. 
This shown the two enlargements Fig. HCN and DCN 
band. one goes shorter wavelengths, first the individual lines become 


4 
q 
a 
q 
q 
fi 
4 
4 
? 


844 CANADIAN JOURNAL PHYSICS. VOL. 35, 1957 


more and more diffuse, then only two diffuse heads are recognizable, and 
finally even these merge into one broad diffuse peak. The diffuseness sets 
earlier and more pronounced for HCN than for DCN. For this reason, more 
sharp progressions are found for DCN than for HCN. particular, will 
shown below, three these progressions extending from 1700 1550 
form separate band system which shall call the system distin- 
guish from the main system which designate a—X. bands the 
system have been identified with certainty HCN. The strong band 
system below 1550 described subsequent paper will referred 

Tables and the wave numbers and assignments the heads all 
the observed bands HCN and DCN, respectively, are given for the region 
2000 1550 the case sharp bands, both and heads are listed; 
for the diffuse bands, only one head or, the head longer pronounced, 
the position the intensity maximum given. These latter wave numbers 
are clearly much less accurate than those the sharp band heads. The 
intensities given the second last column are rough estimates. 


VIBRATIONAL ANALYSIS 

shown Fig. and Table two fairly long progressions are immediately 
obvious the a—X band system HCN. The spacing both progressions 
about 940 and decreases regularly toward shorter wavelengths. 
addition, five weaker progressions with similar spacings can readily dis- 
tinguished. Similar progressions are found DCN. 

One the strong progressions HCN extends the long-wave side 
52278 while the other extends 53767 conclude therefore 
that the second progression one which, addition the vibration 
excited the first, another vibration excited fixed amount. The 
separation 940 successive bands either progression corresponds 
all probability vibration the upper state similar the bending 
vibration the ground state. assume that the first 
progression other vibration excited, that is, might described 
(000). the second progression the additional vibration has the 
the lower state. Thus the second progression may 
and possibly (000) may assigned among the diffuse bands the 
shortward part the band system (see Table 

gression the bending vibration with maximum intensity far removed 
from the first band can accounted for only the assumption that the 
equilibrium angle the molecule the excited electronic state different 
from that the ground state, i.e., that the molecule bent the excited state. 
the molecule were linear the excited state, the first band progression 
the bending vibration should far the strongest, which contrary 
observation the identification the 940 interval with the bending 
vibration correct. 
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Further independent and conclusive evidence for the non-linear configuration 
HCN the state comes from study the weaker progressions just 
the case Two these weaker progressions HCN have fairly 
constant separations (735-750 and 645-605 respectively) from the 
main progression (000). impossible assume the presence 
two further vibrations this magnitude the upper state. One these 
weaker progressions might have been accounted for the progression 
molecule linear both upper and lower state had the correct separa- 
tion from the main progression, which 712 and which should con- 
stant. Since the observed weak progressions not fulfill these conditions, 
different explanation must sought. 

Both the occurrence two rather than one weak progression, whose 
separations from the main progression lie either side the bending fre- 
quency the ground state, and the slight variability the separations 
with can immediately accounted for the assumption made that the 
molecule non-linear the upper state. Under this assumption, HCN would 
nearly-symmetric top the upper state with one very small moment 
inertia, and one would have very coarse rotational structure corresponding 
the small moment inertia. Fig. gives schematic energy level diagram 


onw XK 


on 


Fic. Vibrational levels and coarse rotational structure the upper and lower states 
the bands HCN. 
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showing two successive vibrational levels the upper state with the coarse 
rotational structure, well the three lowest vibrational levels the 
ground state. The main bands must correspond transitions from 
the ground state the levels the upper state since they 
are perpendicular bands (see Section E). The bands the two weaker pro- 
gressions under discussion are interpreted transitions from the first excited 
is, these bands are also perpendicular bands with +1. this way 
the fact understood that the interval not observed 
band separation anywhere the spectrum. Rather, two separations, one 
larger and one smaller than 712 are found and these separations vary 
because the variation the rotational constant with (see 
below). 

addition the two weak progressions just discussed, which may 
designated and there weak progression 
corresponding the second main progression 
Only one member the progression has been found for HCN. 
Finally, two still weaker progressions are observed which can readily inter- 
preted arising from the second excited vibrational level the 
ground state. shown Fig. according the selection rule +1, 
three such progressions may expected: 
The two observed progressions fit the first and last 
these. The first progression has the same upper states the main progression 
and the wave number differences corresponding bands 
should constant (unlike the other progressions) and equal the 
wave number the infrared band, which 1411.42 The observed 
differences the ultraviolet bands are 1411.60, 1411.39, and 1411.42 
which agree very satisfactorily with the infrared value. 

The assignments all observed HCN bands are included Table The 
well represented the formula 
(1) v= §2277.7 9+ 1505.901' — 


The wave numbers calculated from this formula are given last column 
Table Except for the band the agreement with the 
observed wave numbers for all the sharp bands well within For 
the diffuse bands, particularly those for which and heads merge and only 
one peak was measured, the agreement much poorer. This includes all 
bands with and and therefore the value the coefficient the 
term quite uncertain. Attempts improve the agreement were not 
successful possibly account the overlapping system which other- 
wise was not definitely recognized HCN (see below). The rather large 
deviation between calculation and observation for the (000) band 
far outside the experimental error. must ascribed either vibrational 
perturbation wrong possibly the band belongs the 
HCN. 


850 CANADIAN JOURNAL PHYSICS. VOL. 35, 1957 


the upper state the main bands there small amount rotational 
energy due the term (see Section since Thisamount 
has been included formula (1), that is, the formula holds for the 
rotational levels the upper state. The weaker progressions 
Fig. can obtained from formula (1) subtracting 


(see Douglas and Sharma 1953). Since will determined from the 
observed corresponding main bands and bands (see 
Tables and X), calculated wave numbers the bands are not 
and thus drops out. 

DCN the analogues the two main progressions HCN start only 
appreciably shorter wavelengths agreement with expectation account 
the narrower Franck—Condon region. turns out extrapolation lower 
values and ascertaining which extrapolated DCN band corresponds 
the 000-000 HCN band that the first measured DCN band has 
addition the progressions (00°0) and further 
fairly strong progressions are observed shorter wavelengths, three which 
These progressions are much more clearly developed DCN than 
HCN since, mentioned before, the diffuseness the bands sets only 
shorter wavelengths DCN. Therefore, possible have somewhat 
more complete vibrational formula than for HCN. This formula (again 
for the rotational levels the upper -state) 


(2) v= 52404.7 9+ 1505.80,’ — — 4.20572"? 


The assignments the DCN bands and the wave numbers calculated from 
(2) are included Table II. The wave numbers the bands (the 
weak progressions) are again obtained from (2) subtracting 
Checkland and Thompson 1955), and 1138.08 (see Allen, 
Tidwell, and Plyler 1956). 

There one additional fairly strong and very long progression DCN 
which first was thought the progression the a—X 
system, that is, correspond the excitation one quantum the C—D 
stretching vibration. However, such interpretation would have led quite 
anomalous cubic and quartic anharmonic constants because the second 
difference this progression less than half that the other progressions 
and, moreover, constant, unlike the behavior the other progressions for 
which the cubic term (2) quite noticeable. This interpretation was there- 
fore abandoned favor the assumption that this progression due 
another electronic transition (see Section E). weaker but very distinct 
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progression accompanies the main progression and readily inter- 
preted the progression the system. still 
weaker and somewhat fragmentary progression (01'0). 

The main progression can represented the formula 


while the progressions are obtained again 
The numbering (3) quite uncertain since 
the system was observed only for DCN. The reason that definite 
evidence for the system has been found HCN must sought the 
fact that diffuseness the a—X system sets earlier than DCN. This 
makes more difficult ascertain the presence the weaker bands the 
system. addition, the bands HCN are very probably much 
more diffuse than those DCN thus increasing the difficulty detection. 

Neither for the system nor for the a—X system can the zero order 
frequencies the upper state determined since not all the anharmonic 
constants are known and, for the system, even the absolute numbering 
unknown. rigorous discussion the isotope effect therefore not possible. 
But the differences the values (the coefficients the a—X 
systems HCN and DCN are the expected order. The wave number 
difference the (000) bands HCN and DCN 123.3 and 
must equal the difference the differences the zero-point energies the 
upper and lower states. The difference the zero-point energies HCN 
and DCN the lower state (using the data Douglas and Sharma 1953) 
Therefore follows that the corresponding difference the 
upper state must 472.1 this difference, according (1) and (2) 
about 105 accounted for due and leaving 367 for 
Thus the difference between and must about 730 
This value about the same magnitude the ground state, which 
what might have been expected. Thus the assignments the (000) 
bands HCN and DCN are confirmed. 


ROTATIONAL ANALYSIS 

The fine structures bands HCN and bands DCN have been 
measured and analyzed. The wave numbers the lines are given Tables 
IV, and All bands have strong branches addition and 
branches. 

The bands the main progressions etc. 
both the a—X and systems are similar every respect 
bands linear molecule shown the examples Fig. The numbering 
the lines many the bands obtained inspection since the first lines 
each the branches, R(0), Q(1), and P(2), are well resolved. The numbering 
obtained confirmed unambiguously the agreement the combination 
differences with those obtained from the infra- 
red bands measured Douglas and Sharma (1953). This agreement 
shown for the case DCN Table VI. 
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53217.90 54140.15* 
214.04 127.03 136.44 
190.12* 206.34 113.28 128.76 152.71* 
179.50* 199.12 229.74* 102.38 
167.46 190.12* 112.22 
153.72 179.50* 076.60 101.32 
146.40 222.24 095.22 145.28 
145.06 
136.98 129.69 
128.32 
119.14 
53998.62* 
20.13 
53998 
910.59* 917.66* 929.05* | 760 .60* 767.44 779 .06* 
869 .92 890.43 930 .33* 720.00 739.78 780. 20* 
863.38 884.98 929 .05* 713.34 734.16 779. 06* 
817.29 842.74 690.60 
681.38 
905.38 672.46 
789.81 814.91 662.62 
652.14 


P(J) 
52277.65 
255.55* 
250.16* 
9 | 239.35 259.61* 288.41* 
226.71 
220.04 
213.01 
205.54 
197.82 
190.29 
23 

911.24 
920.29 
908.74 


*Blended lines. 


| 
| 
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TABLE (continued) 


WAVE NUMBERS THE LINES THE a—X BANDS HCN (continued) 


or 


524.12 
523.79 462.05 
376.94 
515.03 453.17 372.22 


| | 
| | | 
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TABLE III (continued) 


| 
| (030)°—(01!0) 


(040)°—(0110) 


(050)°—(01'0) 


WAVE NUMBERS THE LINES THE a—X BANDS HCN (continued) 


173.82 224.94 
179.70 

936.13 
930.39 
910.47 
902.71 
894.93 


854 
791.19 
781.73 


HERZBERG AND INNES: 


TABLE III (continued) 


WAVE NUMBERS THE LINES THE a—X BANDS HCN (continued) 


ULTRAVIOLET ABSORPTION 


470.40 
464.70 
530.85 
524.79 
521.69 
0 | 
55318.67 
194.34 
13 | 


16 


855 
| 
17 | i 
= 
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WAVE NUMBERS THE LINES THE BANDS HCN (concluded) 


TABLE (concluded) 


(020)'— (02°0) 


(030)'!— (02°0) 


(040)!—(02°0) 


11 | 677.9 699 .96 738 580.01 601.08 639 .90* 457 .53 478.09 
16 668.52 726.76 59.20 | 
17 660.90 561.15 
(020)*— (0220) (0220) 
6 909 .55 837.42 | 
831.52 
827.98 
814.71 
13 | 809 . 
791.72 
790.74 
785.12 
17 | 784.23 
| 


tUnresolved doublet. 
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TABLE IV. NUMBERS THE LINES THE a—X BANDS DCN 


SPECTRA. 


857 


848.06 854.30 864.15 562.7 572.78 
| 
5 925.40 934.82 950.56 587.77 596.87 613.23 | 231.04* 239.69 256.02 
918.94 931.14* 953.66 | 581.62 593 .30* 616.16 224 .94* 236.09 258.99 
161.60 
156.43* 


*Blended lines. 


| 

q 
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TABLE (continued) 


WAVE NUMBERS THE LINES THE a—X BANDS DCN (continued) 


56727 .88 58170.42 
724.68 58162.86 167.28 174.42 54666.62 
157.62* 165.24 177.74 660.84* 669.72 
093.96 
082.56 178.44 
176.92 
175.34 
171.48 
733.25 434.59 094.38 
730.29 431.23 088.62* 
724.04 424.49 076.75* 
723.08 
719.72 420.54 070.80* 107.28 
412.01 098.65 
407.14 093.77 
401.83 
070.80* 
063.79 


j 


HERZBERG AND INNES: 


TABLE (concluded) 


ULTRAVIOLET ABSORPTION SPECTRA. I 


WAVE NUMBERS THE LINES THE a—X 


859 


BANDS DCN (concluded) 


(120)2— (0110) 
54807.90 


54812.04* 
54742.71 749. 54807 .16 813.81 
734.59 778.70 
729.66 768.98 792.89 
718.94 782.74* 
776.76 
775.98 
771.71 
767.44 
762.86 
759. 26 


Blended with a line of the stronger (120)° band. 


WAVE NUMBERS THE LINES TWO B—X BANDS DCN 


TABLE 


(040)'— (00°0) 


(070)'—(00°0) 
vo = 59642.21 


P(J) R(J) 


57519.55 

57519 .22 523.98 
511.50 517.84 527.22 
508 .72* 516.48 528.64 
505.52 514.96 529.90 
502.10 513.22 530.78 


531.74 
532.58 
533 .02 


510.96 
508 .72* 
505 .96 
502.92 


59642.10 59646.30 


59638. 16 641.43 648.31 
634.54 640.28 649.71 
631.28 639 .04 650.90 
627.56* 637.37 651.88 


652.53 

653 .09* 
653 .09* 
653 .09* 


624.34* 635.41 
620.70* 633.03 
616.33 630.51 
611.92 627 .56* 
607 .28 624 .34* 
602.39 620.70* 
616.82 
612.64 
608.12 
603 .26 


|_| 
0 
1 
2 
3 | 
4 
5 
6 
7 
8 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
| 

0 
1 
2 

3 
4 
5 
7 
8 
9 
11 
12 
14 
15 E 

| 
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TABLE 


COMBINATION DIFFERENCES FoR THE GROUND STATE DCN 


From the a—X bands 

(84.53) (84.56) 84.44 
(94.32) 94.12 


“From Douglas and Sharma (1953). ; 
© Values in parentheses are based on blended lines and are not included in the average. 


Further confirmation the numbering provided the occurrence 
number perturbations both the HCN and DCN bands. They occur 
corresponding values the and branches given band different 
values the branch confirming that the upper state which 
type. many instances the peak the perturbations extra lines appear 
and these are included Tables and IV. 

Just the case (Ingold and King 1953; Innes 1954) the com- 
bination differences for the upper state the HCN and 
and Q(J+1)—P(J+1), show very considerable combination 
defect, far larger than can accounted for /-type A-type doubling. 
Again the cause this splitting must sought the non-linearity the 
molecule the upper states and which has already been shown follow 
from the vibrational structure the band systems. the molecule non- 
linear, nearly symmetric top with one very small moment inertia 
about the a-axis (see Fig. and two nearly equal moments inertia 
and about the and axes. The energy levels are given approximately 


where 


this approximation the levels with are doubly degenerate. Since 


? 
860 
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Fic. Structure and principal axes HCN. The axis perpendicular the plane 
the molecule. 
actually the molecule asymmetric top, i.e., this degeneracy 
split, the splitting being large for low (non-zero) values but negligibly 
small for high values. According Wang (1929) and Lawrance and Strand- 
berg (1951) the splitting for given 


where and refer the upper and lower the two components into 
which level given and split. this type splitting that causes 
large combination defect between the branch the one hand and the 
and branches the other, similar the /-type and A-type doubling 
linear (and diatomic) molecules with Here the levels one 
series for totally symmetric electronic state) are for even 
and for odd while those the other series have opposite parities. 
the ground state HCN, which has only one series rotational 
levels, the even levels are positive, the odd negative. Because the selection 
rule that positive combine only with negative levels, the transitions with 
(F, the upper state totally symmetric) while those with 
and branch) the other. 

both the a—X and systems the sign the combination defect 
indicates that the upper levels the lines form the lower series sub- 
levels, This observation shows unambiguously that the upper electronic 
states and are not totally symmetric (A’) but have the species (are 
antisymmetric with respect the plane the molecule). 
non-linear molecule all vibrational levels are, course, totally symmetric 
and thus the character the upper states cannot due vibration. 

upper state implies that the transition moment the transition 
perpendicular the figure axis the nearly symmetric top, conclusion 
that agreement with the observation that for all observed bands 
bands symmetric top molecules (see also Section D). Unlike 
transitions violating these rules have been observed for HCN DCN. 

(6) introduced into (4) and symmetric splitting assumed, one 
obtains for the rotational levels slightly asymmetric top with 
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These formulae represent the rotational levels the present case with very 
good accuracy since both HCN and DCN, even when they are strongly bent, 
are nearly symmetric tops. According the accompanying paper Polo 
(1957) the deviations the energy levels from (7) account the asymmetry 
the rotator can represented quadratic and higher terms J(J+1) 
addition very small correction terms order the constant 
and linear terms (7) where 


(8) 


Since the present case less than the corrections the constant 
and linear terms can safely neglected. The quadratic term J(J+1) 
for again neglecting higher powers 

(9) 

This term can combined with the centrifugal stretching term 
If, addition, introduce the centrifugal stretching terms and 


(7) 


where 


(10), terms cubic J(J+1) and still higher terms, all given Polo, 
case for 25. The cubic term 


has opposite sign for the and levels while for the quadratic terms the 
asymmetry correction the same for and 

the basis (10) the values and the upper states the 
type HCN and DCN bands can now evaluated standard methods 
from the observed lines (see Herzberg 1945). The combination sums 
and the branch were used obtain and 
respectively well Since and are accurately known 
from the infrared 1.47823, for HCN and 
By” 1.20776, for DCN; see Douglasand Sharma 1953), 
the values obtained this way are more accurate than 
those which one would have obtained from the combination differences 
(without the use The final and values for the 
various vibrational levels HCN and DCN are listed Tables VII and 
VIII. The values are given separately for the and levels since 


difference the apparent quadratic term found Herzberg and Ramsay (1955) for 


HCO due the effect the cubic terms. 
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TABLE VII 


EFFECTIVE ROTATIONAL CONSTANTS AND THE VIBRATIONAL LEVELS THE STATE 
HCN 


Level | -——|——— --—— 
| | 
| | 


Note.—D values in parentheses are assumed. 


few cases differences were found which are entirely outside the limits error 
the determinations. 

consist two subbands, one with the other with the upper 
state. agreement with this interpretation, the fine structure the sub- 

The subbands like the main bands consist simple 
and branches. Here again there combination defect, but much 
smaller than for the main bands and agrees within the accuracy the 
measurements with the /-type doubling found for the level the 
ground state the infrared (Douglas and Sharma 1953). Those sub- 
bands that are well resolved show double double and double 
branch, expected since the upper state has K-type splitting and 
the lower state /-type doubling. 

The bands arising from the second vibrational level (020) the ground 
state consist general three subbands 
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TABLE VIII 
EFFECTIVE ROTATIONAL CONSTANTS B AND D OF THE VIBRATIONAL LEVELS OF THE @ AND B 
STATES DCN 
(020)! | 1.0826 1.1409 1.1117 | 8.7 (8.7) 6.99 | 
a (03 0)! | 1.0801 1.1423 1.111, 8.1 8.1 7.36 0.7 0.7 | Pert. PRJ >8, Q(5) 
a (040)! | 1.0782 1.1471 1.1126 41.1 8.08 3.0 | Pert. Q(13) 
a (050)! | 1.0742 1.1475 1.1108 8.6 (8.6) 7.68 0.9 
a (060)! | 1.0701 1.1504 1.1102 8.2 22 
a (07 0)! | 1.0654 1.1475 1.1064 8.0 8.0 
a@ (1 2 0)° 1.1027 25.5 | 
a (1 20)? 1.1041 12 
a (13 0)° 1.1033 32 | 
| 
a (1 40)° 1.1033 32 31.2 0.8 
(140)! | 1.0696 1.1373 1.1034 14 7.80 -0.1 0.1 
a@ (240)! | 1.0631 1.129; 1.096; 8.7 13 7.54 l.s 
| 
B (04 0)!| 1.0565 1.1099 1.0832 6.60 | 
— 
Note.—D values in parentheses are assumed. 


and which are respectively the types and 
The structure the first type identical with that the main bands. 
That the second similar, since the doubling the branches expected 
this type has not been resolved (see below). II—A subbands were sufficiently 
strong sufficiently free overlapping for their structures analyzed. 
Further confirmation the subband types comes from the number 
missing lines near the band origin. all cases which there sufficient 
freedom from overlapping the first lines the branches are agreement 
with expectation (see Tables III and IV). particular, while the main 
and the last named addition R(1), Q(1), and Q(2) are 
missing. This leads very characteristic gap between the first line the 
branch and the first line the branch (see Fig. 1). 
The rotational levels nearly symmetric top with 
(see (10) above) can represented 


unlike 


those with 


where 


While for the coefficient J(J+1), changed markedly 
compared the symmetric top and differently for the and levels (eq. 
(10)), here changes occur only the quadratic and higher terms and only 
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they differ for the and levels According Polo (1957), 
the quantities are again determined the asymmetry parameter For 


Small corrections the coefficients the principal terms (12) are also 
the order (A—B)e?; but these can neglected entirely compared 
and 

The constants and for 2,3 could again have been evaluated 
the standard way from the observed branches. However, actually some- 
what different procedure was followed since most the bands with 2,3 
are weak and often lines the and branches are not available for corre- 
sponding values. Therefore the appropriate quantities 
known from the infrared were added the wave numbers the lines; the 
resulting term values the upper state were averaged for each and, for 
separately for the and sublevels. These average term values were 
then fitted equation the form (12) yield and values. The 
resulting values are included Tables VII and VIII. The J-independent 

comparing (11) with (13) immediately seen that the average 
level were negligible. Conversely the observed agreement between 
and within the accuracy the determinations (which 
pared and cannot determined from differences the values. Within 
the inaccuracy introduced this neglect, obtain directly from and 
for the values the rotational constants and Their values 
are collected Tables and XII. 

columns and Tables VII and the asymmetric top corrections 
are included calculated according (9) and (14) from the final 
values Tables and XII. Throughout, the values for DCN are 
larger than the corresponding values for HCN since the deviation from the 
symmetric top greater. spite the rather low accuracy the observed 
values they show, expected, the same trend the values. 
particular the relatively high values correspond high values (see, 
e.g., the (000)°, (010)°, (030)° levels HCN and the and (140)° levels 
DCN; note that according (9) and (14) There are, however, 
few cases (e.g., the (000)! and (020)° levels HCN) which much 
larger than presumably owing perturbations. The differences 
give the centrifugal stretching constants D,, which are listed the last two 
columns Tables VII and VIII. The scatter these values considerable 
but, apart from the few values that are clearly perturbed, seen that 
smaller for DCN than for HCN, again expected. 


» 
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Since the same for the and levels accuracy better than 
+0.01 the difference between the must arise 
from difference the values the and levels. shown Polo 
(1957) the formulae given Kivelson and Wilson (1952) for centrifugal 
stretching asymmetric top molecules indeed give different effective 
for and F,, but not easy predict the magnitude this effect. 

For there difference between and which indeed the 
reason for the A-type doubling for Its magnitude 


The observed line splitting the best-resolved subbands after sub- 
traction the splitting due the state fits quite well with that predicted 
from the values (see Tables and VIII). agreement with the negative 
sign the sign found negative the few well-resolved cases 
((010)? and (020)? HCN and (040)? DCN). 

the evaluation and obtain the same time the origins 
the subbands (see Tables I-V), which turn can now used for evalua- 
tion the rotational constants According eq. (12) the separation the 
origins the subbands and given vibrational transition 
(see Fig. 3), that is, the separation the levels and 
according (12), can obtained the separation the 
and subbands. While subbands have not been observed, their 
positions can obtained subtracting the case HCN, AG” 
the case DCN, AG” from the corresponding 
main bands. Similarly the separation the levels and 
which according (12), can obtained the separation 
origins the and subbands respectively, have accord- 
ing the above discussion 


[vo (11) —vo(Z)] = A-—B-D,. 


are listed for the various levels HCN and DCN. addition, Table 
for DCN the values are given since several cases they 
appear preferable For all the levels HCN and DCN 
one finds 


expected from (16) the centrifugal stretching constant positive 
and appreciable. The resulting values A—B and are given the last 
two columns Tables and seen that for HCN, increases 
rapidly with v2; for DCN, this increase much less pronounced. This difference 
behavior similar that (see above). For the levels the effect 
cannot clearly established. was therefore assumed zero. 


| 


—vo(Z)] vo(II) —v(Z) A —B Dr 

Level 

a(0 20.80 21.41 0.20; 
10) 23.31 23.42 0.11) 
20) 25.57 25.82 (0.25) 


ROTATIONAL CONSTANTS 
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TABLE 


ROTATIONAL CONSTANTS A-B DETERMINED FROM SUBBAND 


ORIGINS FOR THE VIBRATIONAL LEVELS THE STATE HCN 


TABLE 


FOR THE VIBRATIONAL LEVELS THE AND STATES DCN 


This value based the head the (020)° band. 


13.66 has been used rather than the less accurate value. 


has been included these average values. 


: }[vo(A) —vo(Z) 4[vo(4) —vo(I1)] vo( Il) —vo(Z) A —B Dr 

a(0 40) 18.04; (17.97) 18.28 0.078 
a(110) (12.7) 12.5 
a(1 20) 14.36; (14.01) (15. 43) 15.3 
a(140) (17.72) (17.56) 18.21 18.37; 
a(150) (20.55) (20.1) 
a(2 40) (17.30) 18.2; 0.24 
(20.35) (20.14) (21.0) 21.2 
a(3 40) (15.64) 16.5 
a(3 50) (19.34) 
40) (13.46) 13.5 
(14.85) (17.1) 14.7 
100) (17.1) 17.0 
120) (19.4) 18.6 
20.6 
150) 20.0 

Values based heads rather than band origins. 


> 
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should emphasized again that the asymmetry corrections the terms 
(12) shown the formulae Polo (1957) are 
entirely negligible the present case. 

The values Tables and when combined with the and 
values Tables and XII lead the values given the second column 
the latter tables. 

The rotational constants vary with fairly normal manner, 
similar that found for and Fig. and this variation, 
function represented graphically for the states HCN and 


-! 

cm ° 

v 
(a) (b) 
15 
By 
1.05 
1.20 


Fic. Dependence the rotational constants A,, the state (a) HCN 
and DCN. Note the difference scale the upper and lower parts the diagrams. 


DCN respectively. While decreases with increasing v2’, and increase. 
The (non-linear) increase very marked (note the difference scale 
may noted that similar behavior the rotational constants found 
the ground states and D.O where also and increase with 
increasing (bending vibration) and where particular the increase 
with very marked. This strong increase with the bending vibration 
clearly connected with the gradual approach situation similar that 
linear molecule. contrast, all three rotational constants show normal 


i 
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TABLE 


ROTATIONAL CONSTANTS, MOMENTS INERTIA, AND INERTIAL DEFECTS THE STATE 
HCN 


(010) 
(020) 
(030) 
(040) 
(050) 
(100) 
(110) 
(120) 27.09% 


Note.—Values square brackets correspond eq. (17). 


— 


— 
| 


TABLE 


ROTATIONAL CONSTANTS, MOMENTS OF INERTIA, AND INERTIAL DEFECTS IN THE @ AND 8 S?ATES 


DCN 
a(110) 2.00 
a(120) 
a(1 40) 19.47, 1.1712 1.0358 1.687 
a(230) 18.1, 1.542 
240) 19.3; 1.1627 1.0299 1.445 24.072 1.659 
a(250) 22.3 
14.8 
15.8 
17.5 
100) 18.1 1.546 
120) 19.7 
20.0 1.399 


square brackets are extrapolated according eqs. (18) and (19). 


| | 
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decrease with increasing although for A(DCN) this decrease within 
the rather large experimental error and for A(HCN) complicated 
term 

For the state (DCN) the increase with linear while and 
both decrease. must remembered that the vibrational numbering not 
established and that all rotational constants are less accurate than for the 
state. 

The following equations for the constants were found: 


The constant terms these expressions are the Ao, Bo, values. Since 
for the state and and for the state cannot determined, not 
possible obtain the equilibrium values A,, B,, 

From the rotational constants Tables and XII the effective 
moments inertia are obtained according 

and similar relations for and These values are included Tables 
and XII. The last column gives the inertial defect 


(21) A = 


well known, the equilibrium position this defect zero for any planar 
molecule, but deviates increasingly from zero with increasing vibration. This 
borne out the data Tables and XII. The value for the (000) 
level the state HCN well DCN far the smallest, would 
expected. this way the vibrational numbering confirmed. 

Since equilibrium moments inertia cannot determined, not possible 
obtain unambiguous equilibrium internuclear distances and angles. How- 
ever, effective values for the lowest vibrational level the state can 
obtained. Since with the two isotopic molecules there are four independent 
moments inertia while there are only three geometrical parameters 
possible get two (or more) sets parameters whose agreement disagree- 
ment gives indication the effect zero-point vibration these para- 
meters. 

The moment inertia non-linear triatomic molecule readily 
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where mx, my, are the masses the three atoms, 
rxy and ryz the two internuclear distances, and the angle. 
this equation applied HCN and DCN, possible eliminate simul- 
taneously and and one obtains 


where and refer HCN and and DCN. 

Using the values for the a(000) levels obtained from equations (17) 
and (18) (i.e., not the directly observed values but the smoothed values) 
one finds 


Using this value for can determine and separately for 
HCN and DCN from the observed values. 
From Hirschfelder’s (1940) equations one finds 

(25) I, Ip => sin’6 
Substituting the final and values for the a(000) state HCN and 
DCN well the value determined above and combining with eq. (22) 
obtain 


The values for HCN and DCN agree fairly satisfactory way. The 
differences are the order expected account the effect zero-point 
vibration. must also remembered that the moments inertia for DCN 
are extrapolated values (see Fig. and therefore not reliable those 
HCN. The following geometrical data for the state may therefore 
adopted: 


the effect zero-point vibration neglected, that the inertial 
thus the geometrical parameters can evaluated without recourse 
this way find 

for HCN: 1.165A, 121° 28’, 


These values deviate somewhat more from the previous set than the difference 
between HCN and DCN found either case. one were calculate the 
moments inertia back from this second set parameters one would find 
exact agreement for and while would come out incorrectly the 
whole amount the inertial defect, i.e. 11% and 14% respectively for 
HCN and DCN. The first set the other hand gives values which 
deviate from the input values case more than 1.2%. This set 
therefore much better compromise between the conflicting values arising 
from the effect zero-point vibration. 
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the ground state HCN 1.0657, 1.1530 and 
180°. Thus, excitation considerable increase r(CN) takes place 
addition the large change the angle. While the latter accounts accord- 
ing the principle for the length the main progression 
the former accounts for the strong excitation one 
quantum. The change r(CH) seems large enough for excited 
the upper state well. The fact that its excitation has not been found may 
related the particular form the normal vibrations the state, 
viz. the fact that the atom and does not move respectively per- 
pendicular and parallel the bond, possibly predissociation (see 
Section F). 

For the state reliable geometrical parameters cannot obtained since 
has only been observed for DCN and since the vibrational numbering 
uncertain. Assuming the numbering correct, i.e. using the extrapolated 
rotational constants for the state, give the parameters for two 


Unless one were assume that the distance much greater than 
follows that the state r(CN) distinctly larger, distinctly 
smaller than the state. 


PREDISSOCIATION 


mentioned Section the bands the a—X system become more 
and more diffuse toward shorter wavelengths. This diffuseness undoubtedly 
due predissociation the upper state. Very little known yet about 
the details predissociation polyatomic molecules. therefore some 
interest study the predissociation case such the present one where 
the nature the rotational, vibrational, and electronic states involved 
known with certainty. 

the main progression HCN very slight lack 
sharpness the lines barely perceptible for For the lines 
are distinctly diffuse, but still well separated from one another. For higher 
ve’ the lines become more and more diffuse merging more and more into one 
another. the and heads are still distinct but beginning with 
the bands are completely diffuse and show longer distinct heads. 
may set the predissociation limit this progression 55900 How- 
ever, more sensitive criterion were used, such the breaking off emission 
bands, the limit might possibly found much lower. 

the progression, the sets somewhat sooner, 
become entirely diffuse that much sooner than for the 
(000) progression. other words, the excitation one quantum the 
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vibration greatly accelerates the predissociation. This appears not un- 
reasonable since stretching vibration while bending vibration 
(see below). 

Another striking difference occurs between the two 
progressions. This clearly shown Fig. la, which shows the 
band only very slightly diffuse while the band com- 
pletely diffuse. Similarly, the band entirely sharp while the 
band shows already slight diffuseness. The diffuseness the 
main band intermediate between the last two. Thus, there 
considerable increase predissociation probability with increasing 
each vibrational level for which predissociation possible. 

For DCN the predissociation limit shifted very appreciably higher 
energies. the progression diffuseness barely noticeable for the 
(070)'—(000) band and becomes distinct only for the (080)'—(000) band. 
The predissociation limit may set while for HCN was 
55900 the progression diffuseness sets about 
the same energy but increases much more rapidly: the band 
the rotational lines are longer separated. 

Part the difference between the predissociation limits HCN and DCN 
can accounted for due the difference zero-point energy the 
ground state, which 568 The beginning diffuseness DCN must 
necessarily appear shifted this amount higher frequencies compared 
HCN. The remaining difference can explained fairly narrow potential 
barrier assumed between the discrete and the dissociated state, barrier 
that can penetrated more easily, i.e. lower energy, than 
For diatomic hydrides such effect has been observed for example for AlH 
(Olsson 1938; Herzberg and Mundie 1940). 

understanding the effect the vibration the predissociation 
appears possible the use considerations Franck, Sponer, and Teller 
(1932). These authors have emphasized that the complicated Lissajous motion 
polyatomic molecule will affect the predissociation phenomenon. They 
concluded that there should strong dependence the predissociation 
probability the particular vibrational motion the upper vibrational 
level which reached light absorption from the ground state. This 
dependence strikingly illustrated HCN. The bending vibration 
clearly one that does not lead directly the pass which transition 
the dissociated state (H+CN HC+N) could occur, and only after 
considerable time will the Lissajous motion lead the image point the 
pass. Therefore, accordance with the considerations Franck, Sponer, 
and Teller, diffuseness sets very gradually the progression. 
With the excitation one quantum the situation changed: pre- 
dissociation sets earlier (at lower and lower energy) since apparently 
the vibration the upper state changes the Lissajous motion such 
way that the pass reached much sooner. would expected that the 
vibration leads the pass even more readily and therefore predissociation 
even more accelerated than for This consideration may account for the 
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non-observation the spectrum since the corresponding bands may 
diffuse that they are not readily observed. 

explanation the strong dependence the predissociation the 
rotational quantum number requires more detailed study the form 
the eigenfunctions involved and the resulting selection rules. similar 
and even more striking dependence the quantum number was found 
(subsequent the observations here described) Herzberg and Ramsay 
(1955) for HCO where the (linear) excited state only the levels with 
are sharp, those with are entirely diffuse, and those with are 
strongly broadened that they are longer observable. 

From the heat formation HCN (—31.2 given Rossini (1948), 
the dissociation energies and (Douglas 1952), and the heat 
sublimation carbon (Gaydon 1953), one finds for the atomic heat 
formation HCN the value 304.0 kcal. This yields, with the values 
D(CH) kcal. and D(CN) 191 


D(HC—N) 223.9 kcal. and D(H—CN) 113 


The first these much higher than the energy the first diffuse band 
the system HCN (159.8 kcal.). Therefore, the products pre- 
dissociation cannot HC+N spite the fact that excitation the C—N 
stretching vibration the upper state favors predissociation. the other 
hand, D(H—CN) well below the predissociation limit and therefore the 
products the limit must H+CN. 

The energy difference between the predissociation limit and 
large enough that the radical formed the predissociation may 
excited the state which 26.42 kcal. above the ground state. The 
observed strength the predissociation indicates that the radiationless 
transition allowed one and since the state all probability singlet 
state like the ground state follows that the state causing the predissociation 
also singlet state. the other hand only one singlet state, presumably 
the ground state, arises from and therefore the state causing 
the predissociation must arise from other words the 
dissociation products the predissociation process are This 
still leaves near the limit kcal. excess energy which may distributed 
kinetic vibrational (or rotational) energy among the dissociation pro- 
ducts. Thus quite possible that there fairly narrow potential hill 
appreciable height such required for the explanation the change 
the limit going HCN DCN. 

For the bent configuration, gives rise two singlet states 
14’ and Since the observed predissociation strong one affecting 
equally all rotational levels likely homogeneous predissociation, 
that is, since the state the state causing the predissociation the 
forbidden for zero rotation. 


value independent the value assumed for the heat sublimation carbon, 
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DISCUSSION AND CONCLUSION 


the preceding sections two excited electronic states HCN and DCN, 
designated and have been studied. has been shown unambiguously 
that unlike the ground state the molecule non-linear both excited states. 
subsequent paper will shown that the upper state the y—X 
system which extends from shorter wavelengths the molecule 
also non-linear. The non-linearity excited states molecules which are 
linear their ground states was first suspected Mulliken (1941) the 
case CS, but first definitely established the case Ingold and 
King (1952, 1953) and confirmed Innes (1954). 

Ingold and King (1953) have discussed the reasons for the bent structure 
its first excited state and have concluded that due the 
fact that change hybridization from (linear) sp? (120° 
bond angle) the energy lowered states which electron excited 
from the bonding orbital (filled the ground state) what the linear 
case would antibonding orbital. The same reasons apply HCN, 
which has the same number electrons therefore not surprising 
find that the first excited state HCN fact non-linear. 

Independently Walsh (1953) had predicted that both HCN and 
should non-linear their first excited states. However, his detailed pre- 
dictions with regard HCN have somewhat modified the light 
the observed spectra here under discussion. 

the ground state HCN has the electron configuration The 
lowest excited states are obtained taking electron from the 
orbital antibonding orbital. This gives rise four excited states 
linear HCN six states non-linear HCN follows (in the order 
increasing energy): 

linear non-linear 


a”? a’ 14’ 


The orbitals the non-linear case are those used Walsh. have indicated 
underlining which states the molecule linear and which non-linear. 
The question whether the observed excited singlet states (which are pre- 
sumably the lowest excited singlet states) can correlated with these 
predicted states. 

The lowest predicted state the second excited state the right 
the above list. appears natural identify the state with this predicted 
state. The observed type the state thus accounted for and the low 


876 CANADIAN JOURNAL OF PHYSICS. VOL. 35, 1957 


intensity the transition can ascribed the fact that for linear 
HCN would forbidden transition 

Walsh has identified the state with the above list, i.e., assumed 
the molecule linear this state. However, the structure the 
system shows conclusively (as will elaborated paper that the molecule 
non-linear the state, which moreover type Since the state 
lies much higher than the state cannot identified with the first state 
the above list but must all probability correspond the state arising 
from linear HCN. This state predicted non-linear and have 
roughly the right excitation energy (it the exact analogue the state 

According its energy one would tempted identify the observed 
state with the third predicted state the right the above list. However, 
according Walsh this state predicted linear while the observed 
state definitely non-linear. difficult believe that the state could 
correspond the last state the list (derived from linear HCN) since 

Walsh’s predictions are based certain assumptions about the variation 
orbital energy with angle (see Fig. paper). order shift 
the equilibrium the state from the linear the non-linear configuration 
some modification Walsh’s curves will required. may also that the 
discrepancy found just consequence the drastic simplifications inherent 
Walsh’s treatment. 

much more elaborate treatment the electronic structure HCN has 
recently been given Iguchi (1955) using both the molecular orbital method 
with configuration interaction and the 
ionic-homopolar resonance. However, the calculations have been restricted 
linear framework and configurations involving electrons only. The first 
10.08 ev. Presumably these correspond the first two states the list 
above. The correlation with the observed states clearly not very good. Cal- 
culations the effect bending the energies the predicted states 
would important. 

addition HCN and there are three other molecules for which 
change from linear non-linear configuration has been established. HCO, 
with one more electron than HCN, is, expected from considera- 
tions, the ground state which the molecule non-linear. Upon exciting the 
most loosely bound electron from the the component the orbital 
linear configuration predicted and has indeed been observed (Herzberg 
and Ramsay 1955). Similarly has been shown Ramsay (1956) 
that, while the molecule non-linear the ground state, linear the 
upper state the bands. mentioned before, Mulliken had already 
predicted non-linear configuration the first excited state. This configuration 
has now been established detail recent work Kleman this laboratory. 

The numerous perturbations that have been observed the a—X system 
both HCN and DCN are probably produced the otherwise unidentified 


2) 
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but expected state arising from 'A. While seemed possible that they 
might have been caused the state agreement levels was found that 
could account for the perturbations this basis. 

The potential energy the HCN molecule the state function 
the angle bend assuming the molecule fixed plane represented 
schematically Fig. There are two equivalent minima separated 


Fic. Potential function for the bending vibration non-linear triatomic molecule. 
The abscissa the angle bend. assumed that the motion restricted one plane. 


maximum for i.e. corresponding the linear configuration. The vibra- 
tional levels near the minima would expected follow regular course. 
Again far above the maximum they would follow such regular course but 
with half the spacing. Here the levels are essentially those linear molecule, 
alternate levels having even and odd values while near the minima all 
values occur for each value the region near the maximum the levels 
would follow rather irregular course. The observed progression 
the perfectly regular for both HCN and DCN. Since 
addition, least for HCN, includes the (000)—(000) band must corre- 
spond the region the neighborhood the potential minima. the other 


5 
4 
3 
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hand the principle would lead one expect the maximum 
intensity the progression correspond the region 
the potential maximum since the lower state the molecule linear, that 
is, the spacings should quite irregular. This apparent lack agreement 
between observation and expectation can readily explained account 
taken the fact that according the preceding discussion the a—X system 
corresponds the linear case forbidden transition that is, the 
Franck-Condon intensities have multiplied factor which for 
and rises rapidly with increasing Similar considerations apply the 
system, which very probably also corresponds forbidden transition 
the linear case. 

alternative explanation the above apparent discrepancy might 
sought the fact that predissociation which broadens the higher vibrational 
levels strong near the potential maximum that there remains only 
continuous absorption whose intensity difficult compare with that the 
discrete bands. However, the effect the forbidden character the absorption 
the linear case must arise unless the interpretation the electronic structure 
discussed above should turn out incorrect. 

the case the near ultraviolet bands analyzed Ingold and 
King (1953) and Innes (1954) the situation very similar. The upper state 
with non-linear equilibrium configuration all probability corresponds 
state for linear configuration which cannot combine with the ground 
state. Again the vibrational levels are perfectly regular and the region the 
potential maximum not reached. However, the shortward end the spectrum 
has not been well investigated. 
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ENERGY LEVELS SLIGHTLY ASYMMETRIC 
TOP 


ABSTRACT 


For the rotational analysis the band spectra slightly asymmetric tops 
convenient have explicit expressions for the rotational energies terms 
the rotational quantum numbers. the tables available present the 
coefficients the power series the asymmetry parameter are only tabulated 
numerically for each individual rotational level. the present work more 
general expressions are derived, giving the explicit dependence and 
the asymmetry parameter (including terms First-order corrections 
for centrifugal distortion terms and have also been included. 


INTRODUCTION 

the rotational analysis the band spectra asymmetric top molecules 
the energies (usually) large number rotational levels can deter- 
mined. order deduce the values the rotational constants convenient 
have expressions giving the explicit dependence the energies the 
rotational quantum numbers. While for general asymmetric top this not 
possible, the case near-symmetric top the energies can expressed 
very conveniently power series and whose coefficients are functions 
the rotational constants. The main terms the series are those 
and while the coefficients all the other terms are the order some 
power, etc., asymmetry parameter and vanish the limiting 
case symmetric top. Fitting the data power series presents essential 
difficulties, provided the convergence fast. The rotational constants can 
determined from the series coefficients. This procedure has been found 
very convenient the analysis the ultraviolet absorption bands HCN 
and DCN (Herzberg and Innes 1957). the upper state both molecules 
have bent configuration and are slightly asymmetric, near-prolate tops. 

There are tables available for the calculation the energy levels slightly 
asymmetric tops (King, Hainer, and Cross 1943), (Lotspeich, Kraitchman, 
and Solimene 1955), (Schwendeman 1957), which give the energies power 
series certain asymmetry parameter. However, the series coefficients are 
only tabulated numerically for each energy level, and are thus inadequate 
for our purpose. the present work more general expressions are derived, 
giving the explicit dependence the energies and Since the effects 
the non-rigidity the molecule are also important, expressions for the 
first-order centrifugal distortion corrections terms and have also 
been included. 


THE RIGID ROTOR 


The rotational Hamiltonian, the case near-prolate top, can written 
most conveniently the form? 
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the momenta are measured units and the rotational constants are defined 
etc., the energy levels are expressed 
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(1) 


non-vanishing matrix elements the angular momentum operators Eq. (1) 


are 
Since the Hamiltonian diagonal and independent the 
determination the eigenvalues can reduced finding those 


(3) 
where the asymmetry parameter, has been 
(4) 


may noticed that not identical with the parameter introduced 
Wang (1929). They are related 
(5) 
The eigenvalues can obtained standard methods perturbation 
theory (Condon and Shortley 1935), treating the term Eq. (3) asa 
small perturbation. The results given below have also been checked 
alternative derivation using the continuous fraction method (King, Hainer, 
and Cross 1943). 

Each rotational level specified the values (prolate value), 
and (oblate value). the near-prolate case while can 
either (lower) (upper). 


The rotational energies are given power series the form 


n=1 


n=2 

The series expansion has been followed and including the fifth power 
this approximation, for the two energy levels associated with 
given are degenerate, and simple general expression for the coefficients 


are several reasons that perhaps justify the introduction another asymmetry 
parameter. Since the parameter always positive, hile always negative. 
The symbol seems preferable this also used denote rotational constants. The 
inclusion the factor the parameter simplifies the numerical ients the power 
series expansions for the energies. The use the parameter King, 
Hainer, and Cross (1943) has more serious disadvantages. The determinant 
even function (or 6), but not The result that the energy expressions for 
are even functions while for change only interchanges the expressions 
for the two levels with the same These simple properties not obtain the case 
and the energy expressions always odd powers 
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Cmn given. For the energies also depend and expressions 
for each value the pair are given. For odd only the levels with 
need given. Those with J—K can obtained from 
them replacing —e. 

The coefficients Eq. (6) are presented Tables and II. The headings 
the rows and columns indicate directly the factors K?, and 
that multiply the coefficient For even the energy even function 
for odd there are both even and odd powers, but the lowest odd power 

TABLE 
COEFFICIENTS IN THE EXPRESSION FOR THE ENERGY LEVELS 


All expressions and tables are given for near-prolate top 
changing and throughout. should remembered that the near- 
oblate case stands for while J—K stands for K,, and for 
2C). 

may seen from the tables, the contribution the energy the 
terms containing high powers becomes more significant for low 
values. This contribution most important for the level and the 
level the same species but much less for the level 
The energies for these three levels are given including 
terms although should remembered that when high order 
terms become significant, the convergence the series expansion usually 
poor. 

The general expression (6) takes simpler form when the relative importance 
the different terms properly assessed, and only the most significant 
ones retained. For sufficiently small usually enough retain only the 
terms that is, only the first terms the summations Eq. 
(6).4 Thus, for example, the energies for the levels (lowest the pair 


*This equivalent the (unmodified) Mathieu function approximation Golden 
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This expression the same form that for symmetric top molecules, 
and the coefficients can determined the same methods. The energies 
for the levels are obtained from Eq. (7) changing into —e. The 
most significant change occurs the and causes the 
large K-type splitting characteristic the levels asymmetric 
tops. 

THE NON-RIGID ROTOR 

The first-order centrifugal distortion corrections terms the rotational 
quantum numbers are given below. The results are based the treatment 
Kivelson and Wilson (1952). The distortion constants are those used 
Nielsen (1951). the limiting case top, the constants 
and vanish, and the usual expression Dg, and obtained. 


For the sign the underlined terms should reversed. 


(1 1 ) F(Jo, y-2) = F(Jo, 5-2) ngiat [ 16D, 


For J—3, the sign the underlined terms should reversed. 
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The above distortion constants are related the used 
and Wilson (1952) follows: 


Re = 6 4 — 27 rryy— AT 


obtaining Eqs. (8) (13) the identification the axes has been 
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DISCUSSION THE ELECTRICAL PROPERTIES 
COMPOUNDS WITH THE NICKEL ARSENIDE STRUCTURE! 


PEARSON 


ABSTRACT 


valence scheme has been proposed for alloys such CrSe, MnTe, and FeS, 
with the nickel arsenide structure, which indicates why they may semi- 
conductors, while chalcogenides and and antimonides and bismuthides 
Cr, Mn, Fe, Co, and with similar structure are metallic. The valence 
scheme thus consistent with the known electrical properties (and also magnetic 
properties) these phases. then predicts that FeSe and CrS may well 
semiconductors and indicates that all other compounds Cr, Mn, Fe, Co, 
and with the NiAs structure should metallic. 


INTRODUCTION 


Although now possible classify and predict semiconductivity 
compounds which not contain transition metal elements (cf. Mooser and 
Pearson 1956, 1957), this not generally possible compounds containing 
the transition metals, because lack experimental data the electrical 
properties these compounds under precisely specified conditions structure 
and composition which base initial analysis. This particularly 
compounds with the chalcogenides (S, Se, and Te), although the case 
the oxides and mixed oxides much theoretical and practical work has 
been done towards understanding the conduction mechanisms. 

The difficulty which arises the transition metal chalcogenide compounds 
lack priori knowledge the valence state the transition metal. 
That say, not know whether the electrons are involved the 
bonding (formation the valence band) and, they are not directly involved, 
whether they exist discrete states band. is, however, well known 
from the semiconductivity compounds such chalcopyrite cf. 
Austin al. 1956) that the valence band (which must contain even number 
electron states*) can quite independent the shell the transition 
metal. transition metal compound which semiconductor may therefore 
contain odd number valence electrons (i.e., and 
electrons). 

an.attempt develop systematic picture the electrical properties 
compounds the transition metals with the chalcogenides and other 
elements, discuss this paper the electrical properties compounds 
with the nickel arsenide structure. propose valence state for the transition 
metal (Cr, Mn, Fe, Co, and Ni) these compounds which consistent with 
the known magnetic properties and which permits qualitative discussion 
and explanation the known occurrence metallic conductivity and semi- 
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conductivity, and hence allows prediction the electrical properties the 
compounds which have not yet been examined. 


II. THE NICKEL ARSENIDE STRUCTURE AND ELECTRICAL 
PROPERTIES COMPOUNDS WITH THIS STRUCTURE 

The NiAs structure can regarded built from close-packed hexa- 
gonal structure atoms the octahedral holes which lie the atoms 
(the holes are precisely octahedral the axial ratio has the value 1.63). When 
phases with this structure occur the As-rich side the stoichiometric 
composition, vacancies occur the transition metal sites, and when phases 
occur the As-poor side, extra transition metal atoms are admitted the 
trigonal-bipyramidal holes where the coordination five. the stoichiometric 
compound where all the octahedral holes are filled the transition metal 
atoms lie lines along the axis that the atoms are also six-coordinated, 
being located the center right-trigonal prism. 

present know three compounds with the NiAs type structure 
which appear exhibit intrinsic semiconducting properties. MnTe has 
high resistivity 5ohm-cm.), negative temperature coefficient 
resistance the paramagnetic range above about 325° K., and large thermo- 
electric power (Squire 1939; Uchida al. 1956; see also Palmer 1954). FeS 
has fairly high resistivity (up ohm-cm.) and negative temperature 
coefficient resistance (Kamigaichi 1956), and have also been 
informed that CrSe appears behave semiconductor having negative 
temperature coefficient resistance and resistivity ~1/40 ohm-cm. (F. 
Hulliger, private communication). These characteristics have since been con- 
firmed our laboratory sample with resistivity 0.05 ohm-cm. The 
magnetic moments the transition metals these substances have frequently 
been taken indicating predominantly ionic properties but the large thermo- 
electric power (~550 100°C.) MnTe found Uchida al. 
leaves doubt that the conductivity electronic (cf. Joffé 1956, quoting 
Justi 1948 and Meisner 1955). 

With one possible exception, all the other compounds with the nickel 
arsenide structure whose electrical properties have been examined appear 
metallic conductors. 


THE VALENCE STATE THE TRANSITION METAL 


The valence states of, say, the and atoms MnTe must satisfy 
the requirement that the valence band contains even number electrons 
per unit cell, the magnetic requirement five unpaired electrons each 
atom (vide infra), the requirement that the bonding predominantly 
covalent, and also the particular coordination the and atoms 
the unit cell. 

postulate the following valence states for MnTe and similar compounds, 
not only are these requirements satisfied, but are led straightway 
qualitative explanation the electrical properties compounds with the 
nickel arsenide structure. 
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(1) First may note that the resonating bonds (bond number satisfy 
the octahedral and trigonal-prismatic coordination the and atoms 
respectively and furthermore that state which does not involve formal 
electron transfer but uses resonating linear bonds excluded the spatial 
distribution atoms about Te. 

(2) These valence states lead five unpaired electron spins each 
atom, which consistent with magnetic moment about found for 
the atoms MnTe Serres (1947), although Uchida (1956) 
found rather lower number 5.10 These states predict four unpaired 
electrons for the atom CrSe, which agrees with Haraldsen and Mehmed’s 
(1938) conclusion that CrSe contains ions (the magnetic moment observed 
for the atoms stoichiometric CrSe was actually 4.6 The expectation 
four unpaired spins for the atom probably not inconsistent 
with the observations Haraldsen (1941) carried out over rather small 
temperature range the paramagnetic region (i.e. above 320° C.). 

(3) These magnetic measurements not distinguish between covalent 
ionic bonding, provided that the former does not involve the electrons. 
However, the measurements Uchida al. (1956), already 
referred to, leave little doubt that the bonding MnTe predominantly 
covalent. 

(4) The electrons the shell the transition metal are not concerned 
the chemical bonding (i.e. they exist discrete states, usually referred 
the shell model, cf. Zener and Heikes 1956) and hence semiconductivity 
may occur compounds having either odd even numbers and 
electrons. 

(5) With this bond scheme the conditions for semiconducting bonds 
(Mooser and Pearson 1956) are fulfilled alternate atoms bonded together 
have complete and the bonds run continuously throughout the 
crystal lattice. 

(6) The final condition for the semiconducting bond—that there should 
bonds—is satisfied MnTe which the axial ratio 
equal 1.62, FeS where c/a 1.70, and CrSe where c/a 1.63. When 
the transition metal atoms along the axis approach each other too closely 
they interact with the formation intermetallic bonds. That say 
band formed instead the discrete states which exist larger values 
c/a, when the transition metal atoms have greater separation. Metallic 
conduction then occurs owing the formation the band while also 
possible that the orbitals the transition metal atoms play part these 
overlap. The value good criterion the probability 
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bonding changes this ratio come largely from changes rather than 
Present experience shows that can expect semiconductivity com- 
pounds Cr, Mn, Fe, Co, and Ni, with the nickel arsenide structure provided 
that c/a not much smaller than the ideal ratio 1.63;* when much 
less, say 1.50, metallic conduction certainly occurs. 

(7) further requirement for there metallic conductivity that 
the composition should not extend towards the transition metal-rich side 
the equiatomic ratio, for the presence metal atoms the bipyramidal 
holes leads bonds between the transition metal atoms. 


(8) This bond scheme should not compromise explanations antiferro- 
magnetism based superexchange these compounds the electrons 
states the chalcogenide atom are available for this purpose, while ferro- 
magnetism could arise either through direct indirect exchange mechanisms. 
Heikes (1955) has discussed conduction mechanisms these compounds 
maintaining that ferromagnetic alignment the shell spins leads 
metallic conduction and antiferromagnetic alignment should lead 
semiconducting properties. Our bond scheme can, however, lead semi- 
conducting properties without the compound being antiferromagnetic 
state. 

The conclusion outlined and above, that electrons the shell 
the transition metal remain independent the bonding process, and 
not lead metallic conduction provided that the transition metal atoms 
not approach each other too closely, supported the behavior 
the intrinsic semiconductor chalcopyrite, already mentioned. There 
are many compounds with the chalcopyrite type structure which not 
contain transition metals and which the chemical bonding well under- 
stood, known that iron contributes three valence electrons 
and forms bonds with the neighboring sulphur atoms. The valence band 
formed independent the five unpaired electrons the shell iron. 


IV. DISCUSSION ELECTRICAL PROPERTIES 


The suggested valence states for these compounds with the nickel arsenide 
structure indicate why MnTe (c/a 1.62) and FeS (c/a 1.70) have semi- 
conducting properties, and why CrSe (c/a 1.63) probably has. They 
suggest explanation for the surprising difference between the antiferro- 
magnetic semiconductor CrSe and the ferromagnetic, metallic CrTe 
(c/a 1.54), for the latter compound with the smaller axial ratio the 
chromium atoms are close enough for the formation band (direct 
bonds), while the former the electrons the chromium 
atoms follow the shell model discrete states. They indicate why CoSe 
(c/a 1.46) and (c/a 1.46) are metallic conductors while FeS 
semiconductor. Furthermore, they indicate why none the compounds 
Cr, Mn, Fe, Co, with are semiconductors, because all 


*Mr. Fischer has pointed out that the ratio should depend also the size the 
which increases maximum with five unpaired electrons manganese. 
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these compounds the axial ratios lie below the value c/a 1.43 and inter- 
action takes place between the shells the transition metal atoms.* 

Surveying compounds with the NiAs structure whose electrical properties 
not yet appear have been measured, see that CrS (c/a 1.66) and 
FeSe (c/a 1.633) should exhibit semiconducting properties. Compounds 
vanadium and titanium (TiS, VS, TiSe, VSe, TiTe, VP, and NbS 
also have axial ratios greater than c/a 1.63; however, the case these 
substances (and possibly also the case chromium compounds) 
means certain, because the fewer valence electrons which they contain, 
that the valence states will similar that proposed for MnTe above. The 
bonds formed the transition metals may involve orbitals and therefore 
investigation their electrical and magnetic properties will awaited 
with some interest. 

also interesting note that the majority the NiAs phases which 
are metallic conductors and the metallic phases with the hexagonal tungsten 
carbide structure have c/a values that agree fairly well with the calculated 
c/a values taken function rg/r4 (Schénberg 1954), where the 
radius the transition metal atom, the radius the non-metal metalloid 
atom, and given 


However, the compounds which are known predicted semiconductors 
and the chalcogenides Ti, and Nb, the observed and calculated c/a 
values not agree, indicating difference the chemical bonding these 
substances and the remaining nickel arsenide-type phases. This difference 
presumably the absence any direct intermetallic bonding the com- 
pounds having c/a 1.6, which have discussed. 


APPENDIX 


The filled valence band ‘intrinsic semiconductor equivalent 
filled Brillouin zone. there have been some diverse statements the 
literature generally the selection Brillouin zones for intermediate com- 
pounds and particular zone for the NiAs structure (Makarov 1943), 
and furthermore the influence variable axial ratio some interest 
compounds which are intrinsic semiconductors, append the following 
remarks the problem, which was resolved discussion with 
Chambers, Mr. Fischer, and others, whom grateful. 

Each successive Brillouin zone crystal structure must contain 
integral even number electronic states, and the volume occupied k-space 
the mth zone must exactly cover all the faces the zone, 
three zones may only coincident point along line k-space and 
not over plane. therefore follows that the integral number states 


*We have heard that sintered CrSb sample has been prepared which showed negative 
temperature coefficient resistance. Fused samples which have prepared were metallic. 
Furthermore, Haraldsen (1949) reported the phase homogeneous only the 
metal-rich side the stoichiometric composition would not, there- 
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given Brillouin zone independent the unit cell volume and the ratios 
the lengths the unit cell edges, although the relative axial ratios may 
important determining the size the energy gap (if one exists) between 
the mth and zones. 


REFERENCES 


609. 

HARALDSEN, 1941. anorg. allgem. Chem. 246, 169, 195. 

HARALDSEN, and 1938. anorg. allgem. Chem. 239, 369. 

H., T., and 1949. Arch. Math. Naturvidenskab. 
50, No. 95. 

1955. Phys. Rev. 99, 1232. 

1956. Semiconductors for thermocouples (Izdatel’stvo Akademii Nauk 
S.S.S.R., Moscow, Leningrad). 

Justi, 1948. und Leitungsmechanismus fester Stoffe (Vandenhoeck and 
Ruprecht, Géttingen), 88. 

T., T., H., and 1956. Phys. Soc. Japan, 


606. 
Makarov, 1943. Doklady Akad. Nauk S.S.S.R. 40, 191. 
Mooser, and Pearson, 1956. Electronics, 629. 
1957. Chem. Phys. 26, 893. 


1954. Appl. Phys. 25, 125. 

ScHONBERG, 1954. Acta Met. 427. 

1947. phys. radium, 146. 

1939. Phys. Rev. 56, 960. 

E., H., and 1956. Phys. Soc. Japan, 11, 27. 
ZENER, and 1956. Private communication: Scientific paper 8-0628-P9 


Westinghouse Research Laboratories. 


LOW TEMPERATURE RESISTIVITY TRANSITION 
ELEMENTS: VANADIUM, NIOBIUM, AND HAFNIUM! 


ABSTRACT 


Measurements the thermal conductivity from 90° and electrical 
conductivity from are reported for niobium, and 
hafnium. Although the vanadium and hafnium are not pure might 
wish, measurements these metals and niobium allow tabulation the 
“ideal” electrical resistivity due thermal scattering for these elements from 
300° down about 20° Ice-point values the electrical resistivity 
are 18.3 for vanadium, 13.5 for niobium, and 29.4 for 
hafnium. Values for the thermal resistivity vanadium and niobium 
are deduced from the experimental results although for vanadium and more 
particularly for hafnium, higher purity specimens are required before very 
reliable study “‘ideal’’ thermal resistivity can made. For the highly ductile 
pure niobium, the superconducting transition temperature, determined from 
electrical resistance, appears close 9.2° 


INTRODUCTION 


number recent papers have presented experimental data the 
electrical and thermal resistivity some transition elements temperatures 
between and 300° (White and Woods 1957a, and also see Kemp 
al. 1955, 1956, 1957). These data are intended fill gaps our knowledge 
the low temperature resistance these elements and allow possible 
assessment the importance different electron scattering mechanisms, e.g. 
d-band scattering electrons thermal vibrations and the possible 
influence interaction very low temperatures (see, for 
example, Kemp a/. (1955) and White and Woods (1957a) for more detailed 
discussion these aspects). 

The present paper deals with results thermal resistivity measurements 
from 90° and electrical resistivity measurements from 300° 
the elements vanadium, niobium (or columbium), and hafnium. far 
are aware, previous low temperature thermal measurements have been 
confined experiments one sample vanadium from 20° 
Rosenberg (1955), and one niobium specimen from 90° K., also 
Rosenberg (1955). Other data the heat conductivity niobium tempera- 
tures below the superconducting transition K.) have been reported 
Mendelssohn and his collaborators (see compilation Powell and Blanpied 
(1954)). While being very useful from other points view, such data below 
the superconducting transition not give much information about the 
interaction between the conduction electrons and the thermal vibrations, 
which are interested. have been fortunate obtaining specimens 
vanadium and niobium rather higher purity—although means 
pure would desire— than those used Rosenberg and this has enabled 


‘Manuscript received April 1957. 
Contribution from the Division Pure Physics, National Research Council, Ottawa, 
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separate the two electron scattering processes, due impurities and 
thermal vibrations, rather more reliable manner. 

There also appears similar dearth experimental information 
electrical resistance these elements low temperatures. Many years ago 
Meissner and Voigt (1930) determined some resistance values rather 
impure vanadium specimen, and Haas and Voogd (1928) measured the 
electrical resistance hafnium temperatures from 90° 


EXPERIMENTAL DETAILS AND SPECIMENS 


The specimens, generally the form rods about cm. long, were 
mounted cryostat (White and Woods 1955) which the thermal resistivity 
was determined with the aid helium-filled gas thermometers attached 
butyl phthalate oil manometer and the electrical resistivity 
amplifier. 

Some details the source and purity the specimens are given Table 
together with their approximate diameters and physical condition. 


TABLE 
DETAILS SPECIMENS 


Specimen WoT 


99.7% pure 0.0005 in. foil 

ElectrometallurgicalCom- received 3.10 135 120.5 2.57 
pany ~99.9% pure rod 


3.55 Annealed 4.83 196 169 2.85 
1300° 
vacuo 
Nb1 Mackay Inc. (N.Y.) 3X1.5 received 156 
Mallory Ltd. 
Corp. (ductile) >99.9% annealedt 
Nb3 1.59 Drawn 0.47; 2.49 
Foote Mineral Co. 5X1.52 received 4.23 116.5 122 3.47 
0.5-1% 


was formed rolling some 0.0005 in. thick foil into cylindrical coil, which the 
effective cross-sectional area was determined weighing 0.0037 cm.? 

was annealed supporting the wire high vacuum and passing electrical 
current through it, first outgassing the wire and later raising almost the melting tem- 


perature. 


The methods mounting the specimens, i.e. attaching current and potential 
leads (thermal and electrical), varied. 

made rolling thin foil into about inches long, had copper 
leads spot-welded and only its electrical resistance was determined. 

and annealed rod, (for which measurements were not completed 
owing faulty connection the leads), were copper plated near their ends 
and copper leads were soldered with Woods’ alloy the plated surface. This 
method mounting seemed cause rather unreproducible results because 


4 
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the plating did not bond adequately the vanadium. Hence for four 
small copper sleeves were placed the rod and tightened means screws, 
good contact being ensured thin layer silver-conducting paste (‘‘Silver- 
from General Cement Manufacturing Co.) between each sleeve and 
the rod; copper leads had previously been soldered the sleeves. 

three the niobium specimens, Nb1, Nb2, and Nb3, electrical 
resistance measurements only were made. The first two these were mounted 
received, copper plating areas near the ends and then soldering copper 
leads the plating with Woods’ metal. was drawn from rod 
ductile niobium (Fansteel Corporation) and annealed and then mounted 
alumina tube with four small brass screws pressure contact with points 
near the extremities. which thermal conductivity measurements also 
were made, was platinized over small region each end and copper leads 
were soldered the platinum coating. 

were able copper plate the ends the hafnium specimen* first 
electropolishing the surface remove the oxide layer; copper leads were 
then attached with non-superconducting solder. 


RESULTS 


Figs. and are shown the results thermal conductivity measurements 
together with data from Rosenberg (1955) vanadium and niobium. Our 
results are not given they are similar pattern those 
but showed some lack reproducibility owing faulty thermal contact 
between the leads and the specimen. 

might expected for metals which heat transported predominantly 
the conduction electrons the thermal conductivity, proportional 
low temperatures (except the superconducting state where decreases 
markedly below where impurity scattering the electrons 
the conduction electrons impurities static imperfections and 
the electrical resistivity produced the same mechanism, then 
these very low temperatures, where the resistivity due scattering the 
electrons thermal vibrations much less than the impurity resistivity, 
Lorenz law should valid, i.e. 


values po, and given Table indicate that this true for 
the niobium specimen. For somewhat greater suggesting that 
impurities reduce the conductivity due the free electrons such extent 
that lattice thermal conductivity becomes significant. the case hafnium, 
for which experimentally find appears probable that 
lattice conductivity quite appreciable. 

*This hafnium specimen was cut from hafnium supplied Foote Mineral 


Company, through the kindness the United States Atomic Energy Commission. are 
very grateful the Pittsburgh Area Office U.S.A.E.C. for facilitating this gift. 


THERMAL 
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TEMPERATURE (°K) 


TEMPERATURE (°K) 
Fic. Thermal conductivity vanadium and hafnium. 


Fic. Thermal conductivity niobium. 
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deduce the values for ideal thermal resistivity, plotted Fig. 
for vanadium and niobium, have assumed the validity Matthiessen’s 
rule for thermal resistivity, i.e. 


IDEAL THERMAL RESISTIVITY 


100 500 


TEMPERATURE (°K) 
Fic. Ideal thermal resistivity niobium and vanadium. 


The values obtained below about 30° for vanadium are unreliable 
owing the relatively high impurity content, i.e. the high value Wo. 

discussed earlier papers (White and Woods 1957a; Kemp al. 1955) 
and the recent detailed review Klemens (1956), expected for 
ideal free electron metal that for 6/10 


the Debye characteristic temperature. Theoretically expected have 
value about for such ideal metal, but has generally been found 
experimentally the range (see, for example, Klemens 1956; 
MacDonald, White, and Woods 1956) for the alkali metals, noble metals, 
and most transition metals. 

compute from the experimental values for vanadium and niobium 
rather difficult, there appears available data heat conductivity 
near room temperature from which W,, might calculated. However good 
estimate may obtained applying the 
law the form 


| 
= 
n 
= 
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the data for the electrical resistivity room temperature. Our electrical 
resistance data indicate that 295° has the value ohm-cm. 
for vanadium, and ohm-cm. for niobium; thus using the equation 
above deduce 

2.75 cm-deg. for vanadium 


and 2.00 cm-deg. for niobium. 


the case vanadium, this value appears least consistent with values 
for 100° and 500° quoted the Metals (Hampel 
1956). order deduce equation (1) must also know the specific 
heat data Brown, Zemansky, and Boorse (1952) suggest that 250° 
for niobium, and the data Clusius and Buhler (1955) give 390° for 
vanadium. 

obtain W,/T? for niobium and for vanadium. Thus 
find finally: 


From the measurements electrical resistivity, have also endeavored 
separate the two components: constant due impurity scattering, 
and due scattering thermal vibrations. For the vanadium and 
hafnium specimens the experimentally measured apparently constant 
below about K.; values the impurity residual resistivity, for 
Fig. has been obtained assuming Matthiessen’s rule valid, i.e. 


For simple monovalent metal which the conduction electrons are 
expect theoretically that where 5.0 for 6/10; this has 
been generally confirmed experimentally for copper, silver, gold, and the 
alkali metals. For the transition metals investigated previously has been 
found have various values between and with indication some 
very pure elements (where can accurately studied low 
that the index may decrease again about very low temperatures; 
such quadratic dependence has been proposed evidence 
electron interaction. 

the three elements discussed this paper, the purity has not been such 
allow detailed study the ideal resistivity temperatures below 
6/20. fact for hafnium and vanadium, already becoming the dominant 
from Matthiessen’s rule occur, cannot expect obtain truly representa- 
tive values the range near say With niobium, particular 
and the purity much better but the high temperature the super- 
conducting transition—we find 9.25+0.03° K.—makes impossible 
obtain directly value po. From detailed series measurements 
electrical resistance between and 17° K., was apparent that the 
resistance was still decreasing quite markedly with decreasing temperature 


. 
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the temperature was reduced below 10° obtain value for po, various 
algebraic polynomials were fitted the values was found that 
expressions both the types and represented the 
experimental values quite well, the values required for the alternative 
expressions being not very different. However the form the curve 
temperatures below 12° 15° does depend quite markedly the 
particular value that adopted. Therefore definite conclusion can 
drawn regarding the temperature dependence below about 12° 
without more certain knowledge po. Nor can this knowledge achieved 
with sufficient certainty applying magnetic field restore niobium 
the non-superconducting state K., any extrapolation magneto- 
resistance data zero magnetic field bound give questionable values 

more detailed analysis the thermal conductivity hafnium has been 
made, apparent that assumptions dubious validity would have 
made concerning the temperature dependence the lattice thermal con- 
ductivity, the ideal electronic thermal resistivity, before separation 
conductivity resistivity into various components could made. 


DISCUSSION 
apparent from the rather low value found for the index (Fig. 


\ 


10.0 


(Nb) 


IDEAL ELECTRICAL RESISTIVITY (micro-ohm cm) 


m=3.8 (Hf) 


100 500 
TEMPERATURE 


Ideal electrical resistivity vanadium, niobium, and hafnium. 
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from the the d-band the lattice vibrations may important the 
transition elements vanadium, niobium, and hafnium. transitions are 
still the dominant cause the ideal electrical resistance low temperatures, 
Mott (1935, 1936) suggested the case room temperature, then 
might expect (Wilson 1938) the other hand might expected 
that sufficiently low temperatures transitions should become very 
infrequent because the small wave number the phonons with which 
the electrons interact; this so, normal s—s transitions would produce 
electrical resistance varying 7*, with the possible addition quadratic 
term representing interactions. The elements discussed 
this paper not seem sufficiently high purity allow detailed 
Nevertheless set values for for each the three elements has been 
compiled from the experimental data and given Table II. Although these 


TABLE 
VALUES ELECTRICAL RESISTIVITY MICRO-OHM-CM. 
(390° K.) 12.3 (250° K.) 21.6 (210° K.) 
295 19.9 14.5 32.1 
273 18.3 13.5 29.4 
250 16.65 12.3 26.5 
200 12.9; 9.8 20.4 
150 8.7 7.0 14.3 
100 4.2; 8.2 
2.3 2.36 
2.3 
0.38 0.56 1.3 
0.2; 0.56 
0.033 0.08, 0.13 
0.01, 0.04 
Index for 
3.35 2.7 3.8 


*Extrapolated using data Potter (1941). 


may regarded somewhat tentative below 30° K., Matthiessen’s rule 
not being strictly valid, they may represent useful addition our knowledge 
the behavior electrical resistance. The room temperature values 
electrical resistivity seem reasonably consistent with those quoted Adenstedt 
(1952) for hafnium and Potter (1951) for vanadium, although the value 
about ohm-cm. quoted for niobium many physical tables seems 
rather low compared with our value for the resistivity. 
wish draw attention the superconducting transition temperature 
niobium for which resistive measurements give 9.25° Our belief 
that this value perhaps more representative for pure niobium than 
value 8.70° (see review Eisenstein 1954) seems confirmed Preston- 
Thomas (private communication, unpublished), who has found 9.2° 
for ductile annealed niobium both magnetic and resistive measurements. 
The agreement between magnetic and resistive values for obtained 
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Preston-Thomas, and considerably lower value the critical field slope, 
contrast with earlier results (Jackson and Preston-Thomas 1950) 
unannealed specimens. 

hoped that the not too distant future much purer samples elements 
such vanadium and hafnium may available for further investigation 
their transport properties low temperatures. 


wish express our thanks Messrs. Broome and Richardson 
for their technical assistance, Mr. Anglin for helping with much the 
experimental work, and Drs. MacDonald and Basinski for 
helpful discussions various times. 
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THE EMISSION SPECTRUM THE MOLECULE! 


ABSTRACT 


system three red degraded bands 3854, 4228, and 3567 has been 
obtained emission from hollow cathode discharge through helium containing 
little phosphorus vapor and hydrogen. Rotational analyses the bands show 
that they are the 0-0, and bands transition the PH* 
molecule. The state regular with small spin splitting (case while the 
state regular with large spin splitting (case a). For the levels 
the state with and the levels with are predissociated. 
Also the lines arising from the levels the state with are found 
extremely weak compared those arising from the levels. Vibrational 
and rotational constants have been determined and the dissociation energy has 
been found 3.06+0.25 ev. 


INTRODUCTION 


Prior the present investigation, spectrum the molecule was 
known. isoelectronic with SiH, i.e., has electronic structure similar 
that the hydrides the carbon group. may therefore expected 
have strong band system around 4000 corresponding the transi- 
tions the CH-type molecules (Mulliken 1932). investigation the 
spectrum (Narasimham 1956), band was observed 3854 which 
had the appearance hydride band. systematic search produced two more 
bands belonging the same system. Vibrational and rotational analyses of. 
these bands showed that they were the 0-1, and 1-0 bands the expected 
transition the molecule. The method excitation and the 
results the rotational analysis the bands are presented the following. 


EXPERIMENTAL 


The spectrum was excited hollow cathode discharge through helium 
containing little phosphorus vapor and hydrogen. The discharge tube and 
side bulb containing phosphorus were enclosed water bath tempera- 
ture about 50° continuous flow helium 2-3 mm. mercury was 
maintained through the discharge tube. The current through the discharge 
was about 400 ma. The spectrum was photographed ft. concave grating 
spectrograph the second order dispersion 0.35 A/mm. Besides the 
new bands the stronger band systems PH, AIH, and were 
emitted with considerable intensity the spectrum. The rotational structure 
the 4228 band was partly overlapped the band the transition 
search for other new systems and the spectral region 
2000-8000 using the same source, proved unsuccessful. 
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ANALYSIS THE BANDS 


The rotational analysis the band system essentially similar 
that the system the SiH molecule given detail Rochester 
(1936). The 0-0 band reproduced Fig. energy level diagram for 
transition given Fig. assumed that the state belongs 
case and the state case The validity these assumptions will 
evident the following analysis the bands. can seen from Fig. 
there are branches which form close pairs A-doublets. The expected 
first lines the branches are indicated the figure. The observed wave 
numbers the lines the 0-0, 0-1, and 1-0 bands are given Tables 
II, and 


2t 2 
eo ca a 
(4) 
(3) 
(2) 
(2) 
2 


Fic. The energy level diagram (case (case transition. The first lines 
the branches are indicated the diagram. 


State 
(a) Determination and 


The rotational terms doublet state are given the modified Hill-Van 
Vleck expression (Mulliken and Christy 1931) 


where the sign refers the levels and the refers the levels 
and A/B,, being the coupling constant between and 
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Equation (1) gives for the rotational term differences the state 
(2) 
The values were obtained from the differences 


and the values from similar differences lines arising from the 
levels the state. Averages these differences were used for the 
and the values equation (2). and were deter- 
mined from the graph +1) against (J+1)? and 
the resulting values are given Table IV. 


TABLE 


state 
By” = 8.385; B," = 8.1456 

state 
Ay = 1.38 A,’ = 0.82 


(b) Determination and 

for levels, gives for the spin splitting, state approaching 
case 


N(N+1) + y(V+4). 


The spin splittings the state were obtained from the following differences: 


and averages their values were used equation (3). and were deter- 
mined from graph against 1/[N(N+1)]. Two positive values 
were obtained for each vibrational level the state. Therefore the 
state regular. 

The two values for the level the state are 3.8 and 0.2 and 
those are 3.875 and 0.125. each case, the value the 
correct one since was found that the lowest rotational levels, and 
the state form pair with the lying lower (Mulliken 
1932). These values lead the spin coupling constants given Table IV. 
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(c) A-Doubling 

the A-doubling the state were appreciable, could noted 
the fact that the A-splitting the branches would not the same 
the corresponding branches, and similarly for other pairs. However, 
these doublet separations agree well within the accuracy the measurements. 
is, therefore, concluded that the A-doubling state negligible. 


State 

(a) Determination and 

The rotational terms for the state are given equation (1) set 
equal and were again determined from equation (2) using the 

differences for the and the values the 
The resulting values and given Table IV. 

(b) Determination the Spin Coupling Constant 

Taking the appropriate terms for and from equation (1) and 
solving for get 


Small terms involving have been omitted from equation (4). 
equal the following differences: 


average value these differences was taken for the 
right-hand side equation (4) and was determined for each 
value. The values obtained were not constant but were found decrease 
slightly with increasing The value obtained from extrapolation 
equation (4) was used for the determination The 
values obtained for the vibrational levels and the state 
are given Table IV. 

(c) Determination the A-Doubling 

Theoretical expressions for the separations, the A-doublets for 
regular states approaching case were given Mulliken and Christy 
(1931) 


and 
(6) : Avan 


where and are the A-type doubling coefficients. Experimental values 
the separations the A-doublets are given directly the observed splitting 
subbands. From graph against was determined. Sub- 
stituting this value equation (6), was determined. The resulting 
values and are included Table IV. 
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Origins the Bands 
From equation (1), follows that the band origins can calculated from 
and branches the expression 


The band origins determined for the 0-1, and 1-0 bands are given 


TABLE 
DESLANDRES SCHEME OF BAND ORIGINS 
: 
1398.76 


IV. PREDISSOCIATION THE STATE 


Predissociations are observed the rotational levels and the 
state. The 0-1 band 4228 was badly overlapped the rotational 
structure the 4241 band AIH and the expected predissociation could 
not observed. Therefore, the predissociations observed the and 
bands only are given below. The predissociations observed the and 
branches the band are indicated Fig. 


Band 

The lines all the branches the 0-0 band arising from the and 
levels the state with were either missing extremely weak. 
the subband, and were the last lines observed. 
Further lines the and branches with higher rotational 
numbers were not observed. the subband, the last strong lines 
the and branches were those from the levels the state with 
12, namely the and lines. However, the lines 
and were also observed, though very weakly. 
further and lines were observed. and branches the last 
lines observed were those with equal and respectively. 
The intensities the remaining branches the transition were 
comparatively weaker and they were not observed beyond P2(123), 
and 


Band 

the 1-0 band, all the lines arising from the levels the state 
were much weaker than those arising from the levels. 

appears that the sudden breaking off the branches the band 
must due predissociation all the rotational levels with 
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electronic state. the case the state, one must conclude that the 
levels are also predissociated even though sharp breaking off the 
branches observed. 

not all evident what electronic states cause this predissociation. 
Since the breaking off not complete one (i.e., the level which 
above the predissociation limit observed) the predissociation probably 
the behavior the and levels. Since the coupling the state 
close case appears that the only significant difference between the 
two types levels lies the fact that for levels and for 
levels Thus appears that the quantum number has some 
significance the predissociation process. Also the predissociation both 
A-components the rotational levels are affected. Thus rotational levels 
even and odd parity are equally affected the predissociation. The problem 
determine which the possible states PH* capable producing 
predissociation having these characteristics. 

The electronic states PH* arising from the lowest lying atomic states 
are listed below: 


P+'D + H2S — 22+, 2A 


States arising from are considerably higher energy. The state 
arising from the ground states atoms must the ground state the 
molecule and the state from must identified with the 
observed one. The only states which, from energy considerations alone, can 
and states. From the rigorous selection rules for predissociation, i.e., 
Oand follows that the state cannot predissociate both 
components the A-doublets and therefore the state not responsible 
for the observed predissociation. Both from consideration the relative 
positions the potential curves and from the fact that significant 
quantum number the predissociation process, but not for the state 
which has good case coupling, follows that most unlikely that 
the state responsible for the predissociation. the two remaining states, 
and each capable predissociating state even though the 
predissociation forbidden one. seems most improbable that the 
state could predissociate both A-components the levels the state 
for and still leave the levels relatively unaffected. Thus 
seems most probable that the ‘II state responsible for the predissociation. 
not all evident, however, that the ‘II state can give the observed 
type predissociation, particularly the observed difference between the 
intensities the lines arising from the and levels the state with 
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DISCUSSION 


has been assumed this point that the emitter the band 
system the PH* molecule. This can substantiated mainly the following 
considerations: 

(a) The observed spectrum bears close resemblance the band 
system the SiH molecule which isoelectronic with the molecule. 

(6) The vibrational and rotational constants obtained for the and 
states the transition are good agreement with those expected 
molecule. 

(c) The fact that the band system excited hollow cathode discharge 
through phosphorus vapor and hydrogen gas provides further evidence that 
the emitter the band system the molecule. 

The electron configurations the molecule will similar those 
the CH-type molecules. The probable low lying electron configurations 
and the resulting electronic terms therefore, are follows: 


follows from the above electron configurations that the ground state 
should state. The lower state the observed band 
system almost certainly the predicted ground state. 

Because the state certainly very shallow one, possible get 
estimate the dissociation energy According the correlation 
the molecular states with the states the separated atoms, the state 
dissociates into the states the separated atoms while the ground 
state, dissociates into atomic states. The dissociation energy 
the state will then given 

state) 


where the energy difference between the and terms 
P+. Douglas (1957) has obtained fair estimate the dissociation energy 
SiH adopting similar procedure. the present case however, 
lack vibrational levels the state with makes the extrapolation 
the dissociation limit the state uncertain. From the relationship 
between a,, B,, and (see Herzberg 1950) value has been 
determined. With the determined, the dissociation energy the 
state has been calculated 7800 The error such determination 
could large, perhaps much +2000 The dissociation energy 
the *II state corresponding this value for the state, 24690+ 
2000 (3.06+0.25 ev.). This may compared with the dissociation 
energies SiH (3.19+0.25 ev.) and (3.47 ev.). 

The predissociation observed the state (discussed earlier para- 
graph) gives definite upper limit for the energy the dissociation. The 
which gives upper limit 27060 (3.35 ev.) for the energy dissocia- 
tion the state. 


\ 
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THE EMISSION SPECTRUM ALUMINUM MONOCHLORIDE 
THE VACUUM 


ABSTRACT 


The spectrum aluminum monochloride excited high frequency dis- 
charge photographed the Schumann region using meter normal incidence 
grating spectrograph the first order. Six brief new systems bands were 
obtained the region addition three characteristic continua 
the longer wavelength side the discrete systems. The analyses five 
these systems indicate common lower level which identified the ground 
state the AICI molecule. For the remaining system, only the upper state 
vibrational frequency derived. All these systems involve six new, hitherto 
unknown, excited electronic states the AICI molecule. The spectrum, when 
excited discharge from high voltage transformer through 
steaming vapor, found consist numerous discrete emission bands 
the region These bands are different from those occurring 
the high frequency discharge. Because overcrowding these bands, owing 
possibly lack adequate dispersion, analysis presented. 


INTRODUCTION 


The study diatomic fluorides the related group III elements boron 
and aluminum has been the subject numerous investigations recent 
vears. The work Chrétien (1950) the emission spectrum boron mono- 
fluoride has revealed the existence numerous singlet electronic states 
the Schumann region. The term scheme this molecule, presented him, 
has shown striking analogy with that the isoelectronic molecule. 
Recent studies the band spectrum aluminum monofluoride (Rowlinson 
and Barrow 1953; Barrow and Rowlinson 1954; Dodsworth and Barrow 1955; 
Naudé and Hugo 1953, 1954, 1955, 1957) both emission and absorption 
cover extensive wavelength region from 1250. This work has 
led the identification large number singlet and triplet excited elec- 
tronic states the AIF molecule. 

The well-known ultraviolet bands AICI lying between and 2550 
were investigated Bhaduri and Fowler (1934), Mahanti (1934), and Holst 
(1934) emission and Miescher (1936) absorption arising from the 
transition Apart from this early work, recent work the 
emission spectrum AICI has been done Sharma (1951). This author has 
studied the excitation the emission spectrum uncondensed 
discharge the presence helium and obtained two new band systems 
the regions and From the structural appearance 
these systems photographed under the high dispersion the first order 
21-ft. grating, interpreted the systems the violet and green regions 

The present investigations were started with view obtaining analogous 
band systems the remaining diatomic halides aluminum, particularly 
the Schumann region. The results the work described this paper have 
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shown the existence six new discrete band systems the AICI molecule, 
which are here provisionally designated systems II, III, and VI, 
the region preliminary report the approximate vibrational 
frequencies these band systems has been recently communicated Current 
Science. addition the above band systems, discrete emission bands 
the region and three characteristic continua the region 
were also reported. 


EXPERIMENTAL 


The new band systems were excited high frequency discharge through 
vapor using watt oscillator. The discharge tube was ordinary 
pyrex tube 35cm. long and mm. diameter, drawn into adapter 
one end and sealed the open end fluorite window. The adapter end 
the discharge tube was connected Cenco High-Vac pump through 
series absorbing towers. pure sample anhydrous aluminum trichloride, 
supplied B.D.H., was used excite the spectrum. The spectrum was also 
excited ordinary conventional discharge tube having two tungsten 
electrodes fused the two side limbs. The terminals the electrodes were 
connected the secondary transformer. glass condenser was 
connected across the secondary the transformer. spark gap was also 
provided series the secondary circuit the transformer. The spectra 
were photographed the Schumann region using meter normal incidence 
vacuum grating spectrograph* having dispersion 17.3 A/mm. the first 
order. The discharge tube, cemented with fluorite window, was attached 
the slit cone the spectrograph obtaining the spectrum the instrument. 
The spectra were photographed Special plates and exposures 
hours’ duration were found sufficient obtain good photographs 
the spectrum. Aluminum lines which appeared the spectrum were 
employed wavelength standards. The wavelengths these lines were taken 
from the tables and converted into vacuum wavelengths. The order 
accuracy involved the measurement the band heads believed 
the visible and ultraviolet regions, Hilger three-prism glass 
littrow and medium quartz spectrographs were employed. 


RESULTS 


The spectrum obtained high frequency discharge the vacuum 
ultraviolet region shown Fig. seen consist the three 
characteristic continua referred above and six discrete band systems 
designated provisionally systems II, III, and VI. The well-known 
ultraviolet bands and the bands the green region (Sharma 
1951) were also obtained when the spectrum was excited the same source 
and are not reproduced here. The spectrum excited condensed discharge 
from transformer, shown Fig. reveals number discrete bands 
the region These bands are different from those reproduced 


*The grating employed the spectrograph was kindly supplied Professor Siegbahn 
Professor Rao. 
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Fig. and occur the longer wavelength side the discrete band 
systems. Some these bands occur the wavelength region covered the 
three emission continua shown Fig. while the remaining bands extend 
towards longer wavelengths \2365. The well-known band system the 
region also observed the condensed discharge from trans- 
former while the visible region none the band systems corresponding 
the transitions the triplet levels are observed. 

all, six new band systems were identified and attributed the 
molecule from examination the spectrum shown Fig. The wave- 
length, wavenumber, intensity, and classification data the bands belonging 
these systems are shown Table 


TABLE 


Wavelength Wavenumber Intensity System Assignment 
2075.4 48183 0,1 
2041.0 48996 0,0 
2019.8 49510 1,0 
2012.2 49697 2,0 
1990.8 50231 3,0 
1964.6 50901 1,2 
1962.0 50968 0,1 
1934.4 51696 1,2 
1930.7 51795 2,1 
1924.4 51964 
1919.3 52102 0,0 
1902.4 52564 0,1 
1871.9 53422 1,0 
1854.0 53938 3,3 
1839.5 54364 1,0 
1822.6 54867 2,0 


System I.—This system consists five bands clearly degraded shorter 
wavelengths. The wavenumber difference between the first two bands 
48183 and 48664 found 481 This being approxi- 
mately equal the ground state vibrational frequency the AICI molecule, 
the band 48664 taken the 0,0 band. The bands the shorter 
wavelength side are assigned the 1,0, 2,0, and 3,0 bands. The 
differences suggest upper state vibrational frequency about 505 

System II.—System consists only two bands, slightly degraded the 
shorter wavelength side. The stronger the two assigned 0,0 and the 
other 1,0. The approximate w,’ value 
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System system consists bands degraded longer wavelengths 
and occurs three groups which are interpreted sequences. The strongest 
tional assignments the bands are given Table and shown Fig. 
The approximate and values are 413 and 482 

System IV.—This also occurs three groups, the strongest them 
52102 being assigned the 0,0 band. all, five bands belong 
this system. The approximate and values are 537 and 477 
respectively. 

System V.—This system consists four bands, the strongest them 
53045 being identified the 0,0 band. The approximate w,’ and 
values are 377 and 481 respectively. 

System VI.—The most intense group this system seen consist 
four closely spaced bands which belong the sequence. The sequence 
degradation the bands towards the shorter wavelengths. The 1,0 and 
2,0 bands are also shown Fig. The approximate w,’ and values 
are respectively 513 and 481 

All the above brief systems with the exception system have for their 
final level the ground state the AICI molecule, the upper states involved 
corresponding various excited upper electronic levels. Although system 
the 0,1 band not observed, likely that the lower state, common 
with all the other systems, the ground state. 

The three continua the wavelength regions 
and appear the spectrum AICI excited the high 
frequency discharge. this region and towards the longer wavelengths 
the spectrum AICI excited condensed discharge from 
transformer consists numerous bands (about 100 number). All these 
bands occur the short wavelength side the well-known ultraviolet bands 
arising from the transition can seen from Fig. the bands 
seem slightly degraded towards the red the longer wavelength side 
while some the bands the shorter wavelength end appear clearly 
degraded towards the violet. The overcrowding the bands the low dis- 
persion employed the present work makes the analysis these bands 
difficult. such, measurements wavelengths, etc., are givenin 
study the absorption spectrum AICI this region may helpful 
throwing light the analysis these bands. worth while take the 
discussion the nature the electronic states the molecule later 
stage. 
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Wave- 
length 


2359. 
2354. 
2351 
2349. 
2347. 
2344. 
2342. 
2340. 
2338. 
2336 
2333. 
2329. 
2326. 
2325. 
2324. 
2322. 
2319. 
2318. 
2315. 
2314. 
2313. 
2312. 


*Degradation for the bands marked not clear. 
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TABLE II. MEASUREMENTS BAND HEADS 


Wave- 
number 


42278 
42476 
42479 
42528 
42571 
42604 
42662 
42699 
42733 
42759 
42792 
42863 
42933 
42977 
42997 
43016 
43063 
43107 
43139 
43189 
43213 
43230 
43251 
43288 
43358 
43393 
43429 
43459 
43476 
43527 
43567 
43605 
43649 
43691 
43731 
43790 
43854 
43900 
43952 
43985 
44057 
44099 
44170 
44218 
44279 
44320 
44368 
44405 
44452 
44490 
44541 


Inten- 
sity 


wr 


RRR EEO Re OR 


Degrada- 
tion 


Wave- 
length 


2198. 
2196. 
2194. 
2193. 
2190. 


Wave- 
number 


45481 
45535 
45562 
45595 
45660 
45739 
45790 
45832 
45876 
45914 
45987 
46041 
46100 
46157 
46268 
46298 
46356 
46416 
46490 
46496 
46579 
46650 
46743 
46803 
46832 
46942 
47006 
47059 
47088 
47194 
47263 
47350 
47423 
47495 
47601 
47664 
47728 
47794 
47872 
47900 
47966 
48024 
48082 
48200 
48319 
48396 
48480 
48532 
48572 
48612 
48683 
48766 
48866 
48907 
48950 
48984 
49022 
49077 
49101 
49157 
49213 
49276 
49356 


Inten- 


sity 


— 


Line 


Degrada- 
tion 


Line 
Line 
Line 
2186.3 
2183.9 
2181.9 
2179.8 
2178.0 
2174.5 Line 
2172.0 
2169.2 Line 
2166.5 
2161.3 
2159.9 
Line 2157.2 
Line 2154.4 
2151.0 Line 
2150.7 
Line? 2146.9 
Line? 2143.6 
2139.3 Line 
2136.6 
2306. 2135.3 
2304. 2133.3 
2302. 2127.4 
2301. 2125.0 
2300. 2123.7 
2118.9 
2295. 2115.8 
2293. 2111.9 
2291. 2108.7 
2288. 2105.5 
2286. 2100.8 
2283. 
2277.9 2092.3 
2275.2 
2269.8 Line 2084.8 
2264.0 Line 2079.8 
2261.5 2074.7 Line 
2256.3 
2252.0 Line 2060.5 
2247.7 
2245.1 
2243.0 2050.6 
2241.2 44619 2046.4 
2236.8 44707 2044.7 
44920 2042.9 
2221.7 45011 2041.5 
2218.9 45067 
2215.4 45139 2037 
2212.4 45200 
2210.6 45237 Line 2034.3 
2207 45294 2032.0 
2202.0 45413 2026.1 


CHRETIEN, 


1939. M.I.T. wavelength tables (John Wiley Sons, Inc., New York). 
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ABSORPTION SUNLIGHT ATMOSPHERIC SODIUM! 


ABSTRACT 


photometer has been built which sensitive light only very narrow 
spectral region near the sodium D-lines. The light scattered sodium vapor 
low density and the scattered light recorded photoelectric spectrometer. 
The latter resolves the two lines and allows the effect stray light elimin- 
ated. With sunlight the source, the radiation detected the bottom 
the Fraunhofer lines and the residual intensity there may measured. 
found that this intensity decreases slightly (relative the continuum) the 
sun approaches the horizon. This effect attributed atmospheric sodium; 
the amount found this method agrees well with that deduced from twilight 
measurements. The thickness the sodium layer atoms/cm.? varies between 
and winter and about April. The measurements also 
give the residual intensities the Fraunhofer lines the incident sunlight, 
averaged over the whole the sun’s disk. These are for and 
(5.06+0.24)% for D2. 


INTRODUCTION 


Measurements the brightness the sodium lines twilight have often 
been used deduce the ‘‘thickness’’ the layer free sodium atoms 
the upper atmosphere. (By mean the number atoms 
vertical column whose cross section There are two serious uncertain- 
ties this procedure: correction must made for self-absorption, and 
the flux incident sunlight not known accurately. Self-absorption 
produced the sodium itself, and has the effect reducing the apparent 
brightness; several attempts have been made calculate the correction, but 
some doubt remains (Donahue, Resnick, and Stull 1956; Hunten 
Chamberlain 1956). The uncertainty the incident flux arises mostly the 
measurement the residual intensity the bottom the Fraunhofer 
lines; the terrestrial lines are much narrower than the solar that only the 
centers the latter are effective exciting the atmospheric sodium. One 
attempt verify the residual intensities found solar physicists has already 
been reported (Scrimger and Hunten 1955); considerably better values are 
found the present paper by-product. The principal result that the 
absorption the atmospheric sodium has been observed directly; the thick- 
nesses found this method are almost entirely independent those found 
twilight and are not subject either the uncertainties already mentioned. 

The only previous attempt this was made Sanford (1950). 
order avoid the difficulty caused absorption lines the light source, 
studied some stellar spectrograms which the sodium lines were shifted 
the Doppler effect. those available, concluded that the most suitable 
was that Virginis, taken May 27, 1940; but showed line that 
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could definitely attributed telluric sodium. The upper limit deduced 
for the layer thickness was atoms/cm.? Now that more known about 
the seasonal variation thickness (Fig. 4), can see that May far from 
the best month for such measurement, and that January, February, March 
much better. This knowledge was partly responsible for our success, along 
with the use apparatus abie thickness 10° atoms/cm.? 

order limit the sensitivity the region the lines, the sunlight 
was scattered sodium vapor into photoelectric spectrometer which 
separated the lines from each other and from the parasitic scattered con- 
tinuum. Their intensity was found decrease slightly, relative the con- 
tinuum, the sun moved towards the horizon. logarithmic plot against 
path length through the sodium layer gave the thickness from its slope, and 
the residual intensity the Fraunhofer line from its intercept. The results 
have already been reported briefly (Scrimger and Hunten 1956). The residual 
intensities are little larger than found earlier and thus slightly better 
agreement with the values found Priester (1953). The thicknesses are 
little larger than those deduced from twilight without allowance for self- 
absorption, showing that the latter reduces the intensity more than 


EXPERIMENTAL METHOD 


The apparatus sketched Fig. was mounted the front the auroral 
spectrometer which has been described Hunten (1953). The spectral slit- 
widths were about and the scan covered the second order once 


Pivoted Mirror Imege 


t adsorbing 
to Spectrometer 
NO watt 4inch diam. 
== Wi nd 
“ia H 


Fic. Experimental arrangement (shown schematically and not scale). When the 
scattering monitor use replaces the sodium cell. 


per minute. image the sun about 0.2 inch diameter was formed 
front the entrance slit, that scattered light came from all parts 
the disk. The sodium-vapor tube was maintained within 0.1° temperature 
between 145° and thermostat, with the window few degrees 
higher avoid condensation sodium, The tube could quickly slid back 
and replaced scattering monitor designed duplicate the resonant 
scattering phenomenon closely possible except that was achromatic; 
this allowed measurement the intensity the light made. consisted 
three slips glass with few scratches each surface produced very 


~ 7 
alic 47 


920 CANADIAN JOURNAL PHYSICS. VOL. 35, 1957 


light grinding; they were mounted that any reflection would away from 
the spectrometer. The whole arrangement could calibrated against light 
from tungsten lamp turning mirror. Samples the four kinds 
spectrum are shown Fig. order avoid possible heating the 
sodium tube the sunlight, infrared absorption cell cuprous 


Fic. pical spectra: sun, sodium cell; sun, monitor; (c) lamp, sodium cell; (d) 
monitor (gain times larger). Each spectrum the result scan taking seconds; 
longer wavelengths are the left. 


chloride solution was interposed the path. Without it, the temperature 
controller showed slight heating effect; with it, effect could observed. 
Systematic error caused this effect believed negligible. 

The most important point the design the sodium tube was minimize 
parasitic scattering white light. This particularly bad with sunlight because 
the light which excites the sodium only bright the continuum; 
notice the difference between the two sources Fig. The tube was designed 
with small diameter near the windows that the path length the light 
the vapor would small. (If the scattered light must traverse considerable 
amount vapor most will scattered and lost; for similar reason, 
the solar image must formed very near the front window.) The diameter 
farther back was increased that the light diverging from the sun’s image 
would not strike the wall until was well out the view the spectrometer, 
and the end was given the conventional horn-shape. All the walls were painted 
black except for the two windows, and baffles were placed severai places. 
the temperatures used the glass not attacked sodium vapor; ordinary 
Pyrex was used, but the entrance window was flat polished material. The 
sodium was introduced from small ampoule (kindly provided Mr. 
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Schmies the General Electric Company) after the tube had been baked 
temperature and pressure mm. The pressure just 
before sealing off was measured ionization gauge. 
three-stage glass diffusion pump was used. 

single measurement consisted the following steps: 


Several scans were made with the monitor sure that the intensity 
the sunlight was not varying. 


Five scans were made sunlight with the resonance tube. 


One more scan with the monitor gave check that the intensity had 
remained constant during Step 


Five scans were made with the monitor and the tungsten source. 
Five scans this source with the sodium tube completed the process. 


was possible exchange the scatterers the 10-second interval between 
scans the spectrometer; one complete measurement thus took minutes 
actual observation. typical sample each type spectrum shown 
Fig. all but one were recorded the same gain. The exception the 
monitor scan the lamp, for which the gain times larger. day’s run 
consisted such measurements distributed between noon and either 
sunrise sunset; will appreciated that this required unusually fine 
weather and this accounts for the rather small number runs reported. 
While would desirable have runs covering whole year, this was 
not felt worth the very considerable labor required, well the tying 
the spectrometer, which has many demands made upon it. 

may useful report another method which was seriously considered 
and for which apparatus was nearly completed. The sunlight first passed 
through monochromator which selected one the lines and eliminated 
most the continuum except for band about The beam was split 
into two parts, one which traversed sodium-vapor cell; they were then 
reunited the cathode photomultiplier tube. chopper was arranged 
pass first one beam, then the other, and attenuator the reference 
beam allowed the two intensities made equal that there was a-c. 
component the output with the sodium tube cold. When was hot 
would eliminate narrow region wavelength from its beam; the result 
would a-c. signal proportional the residual intensity the Fraunhofer 
line. The intensity the continuum would given the d-c. component 
the output. The chief defect this method that depends critically 
exact balance the two beams, and the balance must maintained for 
several hours. Even this requirement could met mechanically, would 
very difficult guard against the possibility sodium condensing 
the windows the tube, which could not distinguished from decrease 
the residual intensity the Fraunhofer line. The method which was actually 
used just sensitive the presence the atmospheric sodium but far 
less uncertain. However, does require the use unusually sensitive 
spectrometer incorporating refrigeration the photomultiplier with dry ice. 
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RESULTS 


The apparent intensity the lines from the tube affected the 
presence the Fraunhofer lines the parasitic light when observing the 
sun. The correction this effect amounts measuring the line intensity 
from the bottom the Fraunhofer line observed the white light under 
the same conditions; other words, adding fraction the white intensity 
the line intensity measured above the continuum. The fraction may 
derived from the monitor traces which Fig. 2(a) example. This pro- 
cedure discussed Hunten and Shepherd (1954) and more fully Harrison 
and Vallance Jones (1956). Let the corrected intensity for given line 
called and the intensity the monitor trace J,; the corresponding intensities 
for the lamp may called and Then the residual intensity the 
Fraunhofer line given 


this includes telluric well solar absorption and will decrease the sun 
moves toward the horizon the former appreciable. Since the telluric 
absorption was found small, possible define effective absorption 
coefficient for the whole resonance line; this further discussed Section 
Therefore plot the logarithm against relative slant path length 
through the sodium layer should give straight line; the slope gives the 
thickness sodium and the intercept extrapolated zero path gives po, the 
residual intensity the solar line alone. fact the runs did give 
straight lines within the error measurement; the other two will consi- 
dered separately. The path length was calculated using height km. 
for the sodium (Hunten and Shepherd 1954); even km. change this 
height would have little effect. The abscissa used was sec where the 
angle between the layer and the line observation. 

practice the accuracy with which and could measured was not 
high desired, because the intensities were small and the shot noise 
the photocurrent was appreciable. The ratio these quantities should vary 
but slowly during the day, however, and plot the ratios (one for each 
line) was made against time. The ratios used were taken from curve drawn 
smoothly through the mass points. This procedure resulted considerable 
reduction the scatter the final measurements. Straight lines were fitted 
each the final plots the method least squares; all the slopes, inter- 
cepts, and errors given were taken from this analysis except noted later. 

Three examples complete runs are shown Fig. the scatter the 
points them perhaps little smaller than the average. Fig. 3(a) illustrates 
the majority the runs, which gave substantial slope, and Fig. the 
last two, which gave small slope. Fig. 3(c) shows one the two which gave 
line far from straight, suggesting non-uniform distribution sodium 
different parts the layer. The point which the sunlight traverses the 
layer moves several hundred kilometers during run, course. Table 
gives the results each run with their estimated probable errors. The first 
columns show the optical depths for each line; the calibration these 
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(b) April 


Residual 


Percent 


seca Thickness) 


Fic. Examples complete runs: Feb. (good run); April (very little absorp- 
tion); (c) Feb. (evidence sodium 


terms layer thickness discussed the next section. The result given 
was derived giving the line twice the weight the D1, since gives 
twice strong absorption, The optical depths and thicknesses are referred 
the vertical direction. 

three occasions was possible remove the attachment from the 
spectrometer and move the latter outside time get twilight run. Unfor- 
tunately, this was such hurry that there some uncertainty whether 
the grating was set the first second order, and this affects the intensity 
factor 0.6. The likeliest results are given Table 
terms the thickness derived with allowance for self-absorption. Along 
with them are the corresponding thicknesses taken from Table 


923 
° ' 2 3 4 5 6 
18 
(a) Feb, 
° » @ 
12 
18 
° e 
18 
(c) Feb. 28 
\ 
\ ae 
5 ° ar 
\ ° ° 
\ 
2 3 5 6 


CANADIAN JOURNAL PHYSICS, VOL. 35, 1957 


TABLE 


SUMMARY EXPERIMENTAL RESULTS (ALL RUNS EXCEPT THAT NOTED WERE MADE THE 
AFTERNOON) 


Thickness, 

Date Error, Error, Ratio Error, 


pr. 
Apr. —.014 0.00 


Mean 0.59 


*Error 


TABLE 


THICKNESSES DEDUCED FROM TWILIGHT MEASUREMENTS ON THE SAME 
DAY AS ABSORPTION RUNS 


Twilight thickness, Absorption thickness, 
Date (1956) atoms/cm.? 10° atoms/cm.? 


Jan. 4.2 (2.5) 8.9 
Feb. 3.1 6.6 
Apr. 1.5 0.95 


CALIBRATION 


Before the layer thickness can derived from these measurements 
necessary know the effective absorption coefficient for the Doppler-broadened 
lines. The value for the peak the line well known (Mitchell and Zemansky 
1934) but must averaged over the whole line. has been shown recently 
(Hunten 1956c) that some cases the effective optical depth 
where the optical depth the center the line; the approximation 
good and this certainly true here. (The presence hyperfine 
structure does not alter this conclusion, since based essentially the 
area the absorption line, which independent any kind structure.) 
However, the scattering tube had considerable self-absorption that 
accentuated the edges the lines this result might not accurate. Therefore 
was thought desirable carry out calibration the laboratory. 

This was done using the light from the tungsten source after letting 
traverse sodium-vapor absorption cell. The vapor density was varied 
varying the temperature the region between 70° and The length 
the cell was 10cm. and was filled the same way the scattering 
cell. Two separate runs were made which the measurements were similar 


924 
a.m.) 
1.14 5.15 
0.95 6.17 
5.90 
.46 
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possible those made the sun; scans the spectrum were 
made each temperature. The vapor densities were taken from Volume III 
the International Critical Tables, and the mean life the excited state 
used was seconds corresponding oscillator strengths 0.33 and 
0.65 for the two lines. The logarithms the relative intensities the scat- 
tered lines were plotted against the thickness the vapor the absorption 
cell; the points showed considerable scatter but they were consistent with 
the theoretical line already discussed. This result reasonable since the 
temperature the scattering cell was kept fairly low and the effects self- 
absorption and line broadening should have been small. 

With the theoretical relationship and temperature 250° for the 
atmospheric sodium, the relationship between effective optical depth and 
layer thickness follows from standard relations (Mitchell and Zemansky 1934 
Hunten The experimental results for the two lines were added 
and divided giving mean effective optical depth referred D1; 
then the layer thickness The results are given 
Table The probable errors given were derived from the errors for the 
individual lines, but where the two results disagreed greatly the error was 
sometimes increased. 


DISCUSSION 


first all necessary show proof that the effect observed indeed 
due sodium. There little doubt that some atmospheric constituent 
responsible; the only other one likely water vapor, which has many lines 
the yellow and particularly near D2. possible that the wings one 
more these might overlap the wavelengths observed sufficiently give 
measurable absorption, has been suggested different connection 
Bricard and Kastler (1944). However, most unlikely that the absorption 
the wavelength would twice that D1, true sodium. 
While the individual results vary widely, Table shows, the average 
absorbed 0.59 times strongly D2. This satisfactory evidence 
that sodium the absorbing agent. 

much stronger piece evidence found the seasonal variation, which 
just opposite for sodium and water vapor. The twilight results reported 
Hunten are shown Fig. The layer thickness remains large 
from December March, but drops less than half April. the other 
hand, the amount water vapor the probably increased 
factor for the two April runs when compared with the others. This 
based the fact that the winter runs were made days whose tem- 
peratures were the range from —35° —12°C., while the temperature 
for the April runs was above freezing. can seen from Table and Fig. 
that the absorption follows the sodium thickness very well, and behaves 
manner incompatible with the idea that water vapor the absorber. 

detailed comparison the absorption results with twilight ones can 
made with the aid Fig. The twilight intensities have been converted 
thicknesses without allowing for self-absorption, using the factor given 
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THICKNESS 


N M M J J s 
MONTH (Beginning October 


Present results compared with those from twilight. allowance for self-absorption 
has been made the latter. The dashed lines give the limits between which all the twilight 
points fall, and the solid line gives the average. 


Hunten the results themselves are from the same paper. About half 
the winter results agree well with twilight, but the other half are higher 
nearly factor showing the effect self-absorption the twilight. However, 
the effect much less than predicted* (references Section 1). The two 
results for April agree well with the twilight; they are bit lower but 
more than expected from the errors measurement. can concluded 
that self-absorption has negligible effect twilight measurements except 
winter when the sodium layer thickest. Even then the correction does 
not exceed factor The three occasions when twilight runs were obtained 
the same day give similar results (Table II). Again, there indication 
self-absorption factor less than winter and fair agreement 
April. 

The average values the residual intensity the Fraunhofer line, 
are given Table III along with others previously found for the whole sun. 


TABLE III 
RESIDUAL INTENSITIES THE TWO FRAUNHOFER LINES (PER CENT) 
Line Present result Priester (1953) Old result (1955) 
5.90+ .46 6.6 5.0 +0.6 
.24 5.8 4.8 +0.6 
Ratio 0.86+ .085 0.88 


The method averaging Priester’s results (1953) briefly described 


*Since this statement was written, Mack has pointed out that the hyperfine structure 
effect the present measurements. Preliminary estimates show that the calculated self- 
absorption agree very well with the results observed here. 


° 
° 
4 
° 
So 
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Scrimger and Hunten (1955), and the other result taken from the same 
paper. Our new results are little higher than our old ones, but still well 
below Priester’s. Since there apparent reason expect systematic errors 
the method, feel that the new result the best available for the inter- 
pretation twilight measurements. The ratio the two residual intensities 
important quantity also twilight work, and the agreement here with 
Priester very good. The intensity ratio the lines twilight can thus 
not exceed 1.72+.17; self-absorption will usually lower still further. 

The two runs which Fig. one are interesting that they appear 
show non-uniform geographical distribution sodium over distances 
the order few hundred thousand kilometers. the other hand, 
the other nine runs suggest that the distribution more often fairly uniform 
over such distances. The possibility patchy distribution sodium results 
inherent error our method finding the layer thickness, but most 
days this error seems less than that caused not very large signal- 
to-noise ratio. The noise largely shot noise the photocurrent, but probably 
variations the atmosphere contribute too. Non-uniform distributions the 
intensity the night airglow are usually observed and may perhaps related 
those found here; Roach has excellent review the work his 
group this field (1955). 


would like thank Lyon for help with the apparatus and measure- 
ments. 
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IONIZATION CURRENTS ARGON AND HELIUM 


ABSTRACT 


Ionization currents alpha-irradiated argon and helium and liquids 
have been measured over electric field range and conduction pulse 
amplitude distributions have been measured over 300:1 range applied field. 
critical examination theories ionic recombination shows that they are 
unable explain the results the present experiment, and new approach 
the problem columnar recombination suggested. 


INTRODUCTION 


Jaffé (1913) produced the first theory recombination columns ions, 
and showed that agreed with his own experimental data ionization 
currents, but the value this agreement limited the fact that the 
currents were near saturation, having 2:1 variation over electric 
field range 100:1. The theory described ionizable media which the positive 
and negative ions had equal masses and was applied carbon dioxide, air, 
and hydrogen gases and methane and hexane liquids, with which Jaffé 
experimented. 

Gerritsen’s (1948) measurements ionization currents alpha-irradiated 
liquids gave reasonable agreement with Jaffé’s theory for nitrogen and 
hydrogen, but not for argon and helium, which the currents were found 
much more strongly field-dependent. Somewhat better agreement for 
argon and helium was obtained with Kramers’ (1952) modification the 
Jaffé theory. Ionization current experiments gaseous argon, helium, 
nitrogen, oxygen, and carbon dioxide Haeberli, Huber, and Baldinger 
(1953) gave agreement with Jaffe’s theory over limited range field. 

From the pulse experiments performed Davidson and Larsh (1950), 
the number electrons, escaping from alpha particle tracks liquid 
argon was estimated function applied electric field. The range field 
was only 7:1 and their questionable correction for electronic noise prevented 
satisfactory analysis. 

considerable improvement signal-to-noise ratio has been obtained 
the present experiments argon and helium liquids. Pulse size distributions 
which are explicable terms the geometry the ionization chamber have 
been measured, the variation pulse size with field being good agreement 
with the observed ionization currents. extension the field range 
ionization current measurements shows that the Jaffé-Kramers 
model for columnar recombination not applicable liquid argon liquid 
helium. 
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THEORY COLUMNAR RECOMBINATION 


Jaffé’s (1913) theory applied the dense columns ionization formed 
alpha particles materials with electron affinities, which the positive and 
negative ions were assumed have roughly equal masses, mobilities, and 
diffusion coefficients and initially have the same Gaussian density dis- 
tribution standard deviation about the particle track. 

This distribution was modified diffusion, transport ions due 
applied field, and the removal ions recombination. Jaffé solved the 
resulting differential equation first neglecting recombination and then 
correcting for it, that ions originally present finally escape: 


where 


these expressions the Hankel function the first kind, order 
zero, the coefficient ionic recombination, the diffusion coefficient 
for both positive and negative ions, their mobility, the angle between 
the direction the applied field and the particle track, and For 
high fields and low fields equation (1) reduces to: 


Kramers’ (1952) solution Jaffé’s differential equation for ionic density 
was obtained including the recombination term but neglecting diffusion. 
For high fields gives equation (2) and for low fields: 


where numerical factor near unity. 


(3) 


n 
0 


EXPERIMENTATION 


The electrode assembly was essentially that described Williams (1957), 
except that grounded metal ring separated the bakelite supports the 
two electrodes. With grounded ring silver paint surrounding the signal 
lead the outside the glass envelope, this completely prevented any 
leakage from the high voltage element the signal electrode (except through 
the liquid between the electrodes). The observation signals smaller than those 
used the mobility measurements Williams (1957) necessitated further 
smoothing the high voltage supply with screened resistance-capacity 
filters. 

Ionization currents were measured with vibrating reed electrometer, 
driving recorder. The average number electrons escaping recombination 


| 
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per alpha particle was determined from the ionization currents, and the 
source strength, particles per second (measured pulse counting experi- 
ments), since 


being the electronic charge. 

Conduction pulses were measured with amplifier selected bandwidth, 
followed discriminator and counter. The collection ions, due 
single alpha particle, the input capacity, the amplifier was complete 
time short compared with the input time constant and produced 
voltage pulse with sharp rise amplitude The pulse decay was 
controlled the input time constant. Thus the full pulse amplitude was 
observed with amplifier restricted bandwidth, the pulse shape being 
interest these experiments. The advantage this restriction band- 
width the reduction noise level, for which optimum conditions were 
obtained with bandwidth kc. 200 kc. 

Estimation noise level and calibration the pulse amplitude measure- 
ments were made connecting pulse generator the amplifier input 
through resistance large compared with the input resistance. this arrange- 
ment was the same the conditions operation, allowed the effect 
amplifier noise pulse amplitude distributions investigated. Theory 
relevant this given the appendix. 


PULSE SIZE DISTRIBUTIONS 


Wadey (1948) showed that alpha particles from sources deposited 
silver the manner the present experiments suffer very little straggling 
energy. therefore assumed that the pulse amplitude distributions 
observed the present experiments were due differences geometry and 
recombination between the tracks alpha particles different angles the 
applied field, being modified electronic noise. 

The geometrical effect arises from the difference between the positive and 
negative ion mobilities observed liquid argon Williams (1957). The 
positive ions were too slow contribute materially the pulses observed 
the present experiments, that the angle, between the alpha particle 
track and the source electrode determined the average fraction the electrode 
separation, traversed the electrons and therefore the apparent charge 
collection per pulse for electrons actually collected. 

For the source electrode negative this expressed by: 


where the range the alpha particles, being less than Equation (5) 
applies ion column whose weighted center from the source 
the Bragg law (Wilkinson 1950). The ranges polonium (5.3 Mev.) alpha 
particles liquid argon and liquid helium are calculated from the ranges 
the gases 0.052 mm. argon and 0.27 mm. helium. Anianson’s (1955) 
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measurements alpha particle ranges liquids confirm that the inverse 
proportionality range and density extends the liquid region. 

alpha particles are emitted isotropically into hemisphere centered 
the source, then the number emitted smaller angles the electrode 
surface Combined with equation (5) eliminate this gives the 
pulse amplitude distribution: 


For reversed field (source electrode positive), equations (5) and (6) become: 
(5’) n 


The superposition noise peaks whose amplitudes are described 
Gaussian distribution pulses having rectangular distribution indicated 
qualitatively Fig. the cumulative probability, being 
given Fig. shown the appendix that observed pulse ampli- 
tude distribution can approximated the signal-plus-noise curves Fig. 


Noise 


Signal Signal 


— x, X3 Xe 


(b) PROBABILITY PULSES AMPLITUDE 


Fic. effect Gaussian noise pulse amplitude distribution. 


Q(x) 
Signol+Noise 
ignol+Noise , 
4 
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then the minimum, maximum, and average values pulse size, X2, 
for the distribution the absence noise are readily determined. When the 
apparatus has been calibrated the X’s correspond minimum, maximum, 
and average values, From the above geometrical considera- 

tions, with the source electrode negative 


and with the source electrode positive 

Further considerations the effect noise pulse amplitude distributions 
are given the appendix, including possible interpretation results obtained 
Davidson and Larsh (1950). 

Deviations the pulse amplitude distributions from the simple addition 
noise the geometrical effect indicate different degrees recombination 
particle tracks different angles the electrodes. Some evidence the 
effect angle upon recombination presented Section 


EXPERIMENTAL RESULTS 


Demonstration the geometrical effect and measure the range 
alpha particles liquid argon were obtained measuring the pulse amplitude 
distributions for number electrode separations, with both field direc- 
tions. The distributions obtained with electrode separations 0.20 mm. and 
0.060 mm. are given Figs. and which indicate that the geometrical 
effect greatest for electrode separations close the alpha particle range. 
From the results Fig. and equations (6) and (6’) the value 0.049 mm. 
obtained for the range, satisfactory agreement with the calculated value 
0.052 mm. Justification for the procedure idealizing the curves Fig. 
straight lines conformity with equations (6) and (6’) and Fig. 
given the appendix. 

The plot pulse amplitude distributions obtained with electrode 
spacing 0.2 mm. argon (Fig. 2), for which the geometrical effect small, 
shows that the form the distributions was essentially independent field 
and therefore that the dependence recombination upon angle was not 
function field above kv./cm. The pulses were well above the noise 
level, the largest being times the level which the noise pulses were 
numerous alpha particle emissions. 

The ratios the minimum maximum pulse amplitudes, deter- 
mined from the slope the distributions and plotted Fig. have been 
corrected for the geometrical effect and indicate recombination differences. 
That this ratio has almost constant value near unity for fields above 
kv./cm. indicates that recombination not significantly dependent upon 
the angle between the track and the field, for fields above this value. 

Confirmation given comparison pulse sizes with ionization current, 
Fig. which the dotted line represents the current calculated from the 
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FIELD 
ATTENUATION 


PULSE SIZE (ARBITRARY UNITS) 


FIELD 

ATTENUATION 
-43 
+43 


PULSE SIZE (ARBITRARY UNITS) 


Fic. amplitude distribution argon with small range effect. Electrode spacing 
0.20 mm. fraction pulses amplitude greater. 

Fic. amplitude distribution argon with large range effect. Electrode spacing 
0.060 mm. fraction pulses amplitude greater. 


pulse measurements correction for the geometrical effect but not the 
recombination effect. The departure the dotted line from the direct current 
results low fields attributed recombination differences which cause 
similar discrepancies between curves obtained pulse measurements 
different electrode separations. The discrepancies are consistent with ionic 
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= 


3.5 


Fic. Ratio minimum maximum pulse sizes argon, indicating the variation with 
applied field the dependence recombination upon the angle between the particle track 
and the field. 


collection dependent upon the component the field perpendicular the 
alpha particle tracks, the Jaffé theory. 

liquid helium 4.2° ionic mobilities are much lower than for electrons 
argon, the positive and negative ion mobilities only differing factor 
(Williams 1957). This had two effects: firstly, the absence distinct 
geometrical effect separated positive and negative ion contributions the 
pulses, and secondly, longer ionic collection time. Since was not possible 
separate the two ionic components, the amplifier was adjusted respond 
equally both. constant value the ratio minimum maximum pulse 
sizes, 0.53, was obtained, independent the field, for fields between 
and 105 kv./cm. This indicates that although there strong angular 
dependence recombination, independent the field over the range 
these measurements. 

lower temperatures the mobilities the ions helium are greater and 
more nearly equal, certainly eliminating the geometrical effect. All measure- 
ments down 1.4° gave the same values 4.2° 

Fig. gives the values mean ionic collection per alpha particle for both 
ionization current and pulse measurements helium. unlikely that the 
difference between the 4.2° curve and the 1.6° and 1.4° values 
significant, that the recombination ions appears identical helium 
and II. 

the values energy per ion pair for alpha-ionization gaseous argon 
and helium, given Bortner and Hurst (1954), are applicable the liquids, 
the ionization currents were about half saturated the highest fields used 
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Fic. The average number ions collected per alpha particle argon function 
the applied field. 

Fic. The average number ions collected per alpha particle helium function 
the applied field. 


935 
x 
ef 
3 ° 
° 
e 
5 
Electrode 
os* 
2 


936 CANADIAN JOURNAL PHYSICS. VOL. 35, 1957 


(being limited breakdown). The beginning saturation for both liquids 
appears Figs. and both showing departure high fields from the 
roughly linear dependence current upon field low fields. 

Between kv./cm. and helium shows anomalously strong 
dependence current upon field, not evident argon. This has not been 
explained. 

The possibility some electron multiplication argon must not dis- 
regarded, the r.m.s. energy free electrons liquid argon under applied 
field 100 kv./cm. estimated from the results Williams (1957) 
ev., which exceeds the first ionization energy the argon atom. 


DISCUSSION 


The theories Jaffé (1913) and Kramers (1952) can compared with the 
results this experiment, but good agreement not expected since 
these theories clearly assume that the positive and negative ions have equal 
masses. 

examination Jaffé’s conditions for high and low fields and 
shows that the high field region was inadequately covered the 
results Fig. permit comparison with equation (2), but that the linear, 
low field region Fig. corresponds Jaffé’s low field condition. sub- 
stituting the Einstein relationship e/kT into the formula for (equation 
(1)), inserting known values, and expressing volts/cm., find: 


Jaffé and Kramers assumed that the length ion column about 
that liquid argon the order cm., giving: 


from which may take the high field condition 
20,000 volts/cm. 


and the low field condition 
200 volts/cm. 


Fig. shows that, within the limits the experiment, the ionization current 
argon linearly dependent upon field over the entire range Jaffé’s low 
field condition, but the simplification equation low fields (equation 
(3)) does not indicate linear relationship. The Kramers theory does give the 
required linear variation (equation (4)), which somewhat surprising since 
was obtained neglecting the diffusion ions, which, effective all, 
must most noticeable low fields. However, the Jaffé formulae give the 
observed decrease the angular dependence recombination with increasing 
field. 

Both theories give strong dependence ionization current upon the 
ionic mobility and/or diffusion coefficient. There was test for this the 
argon experiments since there was method varying the mobility inde- 
pendently the field, but changes temperature the helium experiments 
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gave changes mobility factor about That the ionization currents 
were observed independent temperature therefore evidence against 
the validity the theories. 

The greatest objection the Jaffé theory its mechanism for charge 
separation. liquid argon about 10° ion pairs per alpha particle escape 
recombination applied field 10° The Jaffé-Kramers model 
thus requires the applied field pull apart two cylinders charge with 
linear charge density ions/cm. (the length the cylinders being 
cm.). The applied field can only compete favorably with the mutually 
attractive force between such cylinders when they are already separated 
least cm., which greater than the diameter which they are con- 
sidered have cm.), that the initial separation remains unexplained. 

That the electrons released the ionization liquid argon not combine 
with neutral atoms form heavy negative ions well established (Williams 
1957). This leads suggestion for new model ionic escape from alpha 
particle tracks, which makes use calculations Erskine (1954). 

Erskine showed that electrons ejected from atoms near the track 
alpha particle are given higher kinetic energies than those from atoms the 
edge the column initial ionization. They are attracted towards the 
alpha particle and emerge the other side the track direction almost 
perpendicular it. observed Auger (1926) and Chadwick and 
Emeleus (1926) are the most energetic such electrons.) Since the same 
process occurs simultaneously for many electrons, expanding cylinder 
charge formed, with the outer electrons moving fastest, and therefore being 
partially screened from the inward attraction the positive charges, which 
for this purpose are considered stationary. this way the faster electrons, 
which have more chance escaping anyway, are held less firmly the 
positive ions. The size the initial column ionization hardly affects this 
argument, since only one atom across, the fastest electrons are still 
outermost. 

The initial expansion the negative charge cloud does not constitute escape 
the fast electrons since they will suffer many collisions before becoming 
free the positive ion column; does, however, provide possible mechanism 
for overcoming the potential energy which would hold together two ion 
columns the Jaffé type. 

Subsequent motion the electrons not easily described, but kinetic 
theory (Williams 1957) shows that they lose their kinetic energy only after 
many collisions, which time they may have moved the side the charge 
cloud which favorable their escape externally applied field. 

Ionic escape more effective fields perpendicular the ion column than 
fields parallel it, leading angular dependence the Jaffé theory. 

The requirement that charge collection independent temperature 
mobility satisfied this model. 
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APPENDIX. NOISE CONSIDERATIONS 


The noise voltage any electronic system has been shown Rice (1944) 
have Gaussian fluctuations. The probability the noise voltage being 


between and any instant 
2 


where the standard deviation the noise voltage. This also the form 
the probability noise peaks having amplitudes between and 

The superposition random fluctuations (noise) series pulses all 
the same amplitude, produces distribution with cumulative proba- 
bility having amplityde greater 


(7) P(X) 


Plotted against cumulative probability paper, Q(X) appears straight 
line, the gradient which determined 

check that this distribution was realized was obtained plotting 
cumulative probability paper the pulse amplitude distribution obtained 
the superposition noise series uniform pulses from signal generator. 
The straight line which resulted gave the standard deviation the noise 


pulses. 
The superposition Gaussian noise pulse distribution: 


which indicated Fig. and integrates the cumulative distribution: 


which shown Fig. 


(9) 
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the standard deviation, the noise known, and hence and 
can determined. 

necessary condition for this analysis apply the pulse amplitude 
distributions described was their symmetry. Unsymmetrical distributions 
cannot treated this way, but the departure from symmetry not 
great, was the case with all the measured distributions, the approximation 
quite good. 

For large signal-to-noise ratios the noise correction particularly simple 
which for 1.5 equal 1/D within 2%. was thus determined 
graphically indicated Fig. extending the tangent Q(X) 
correction was not applied results which gave small value D/c. 

Two further effects noise must considered. Rice (1944) showed that 
the number noise peaks per second exceeding amplitude (greater than 

where and are the upper and lower frequency limits the amplifier band- 
pass. The rate emission alpha particles from the source was 500/sec. 
equation (10), noise peaks greater occur this rate 
amplifier with bandpass cycles/sec. (as used these experi- 
ments), that pulses smaller than about were lost the noise. 

addition their modification pulse sizes, noise peaks appeared the 
distributions their own account, and their amplitude distribution was 
determined experiment with alpha-induced pulses. This distribution 
was subtracted from each the pulse amplitude distributions the first 
process noise correction, but only affected the smaller pulses. 

demonstration that the noise peaks seen the discriminator were 
described equation (10), was plotted against producing straight 
line, which the gradient gave value agreement with that found 
from equation (7). 

The last effect that noise peaks superimposed upon the tails pulses, 
and may regarded correction the assumption that only one noise 
peak could affect each pulse. The bandwidth the amplifier could not 
indefinitely restricted without affecting the pulse amplitudes, that the 
pulses did not completely decay within the dead-time the discriminator. 
This allowed noise peaks too small counted alone exceed the dis- 
crimination level superposition the pulse tails. Since the smaller noise 
peaks were much more numerous than those capable being counted alone, 
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this was significant effect experiments with poor signal-to-noise ratios. 
adequate correction has been found for it, but was found negligible 
for signals exceeding times the standard deviation the noise peaks. Measure- 
ments were made over 300:1 range field without using pulses smaller 
than this. 

the work Davidson and Larsh (1950) the largest signals were approxi- 
mately times the noise level, which appears have been taken the level 
which noise pulses and signal pulses were equally numerous. According 
the above discussion would have become important for signal-to- 
noise ratios (defined this way) about 2:1. Davidson and Larsh noted that 
the distribution pulse amplitudes became much sharper when the signal- 
to-noise ratio was less than 2:1, and this appears have been modification 
‘‘tail noise’’ which caused double counting pulses. the present experi- 
ments, the lower noise level made possible avoid this region. 
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THE DECAY THE POPULATIONS METASTABLE ATOMS 
AND IONS FROM THE SAME D-C. DISCHARGE NEON! 


ABSTRACT 


Measurements electron density and the population the metastable 
state have been made following cessation d-c. discharge neon gas. Deter- 
minations these populations were made under identical conditions accurately 
determined times after termination the discharge and decay curves were 
plotted for pressure range from 0.5 12.0 mm. Hg. low pressures, where 
diffusion the principal factor the decay, both curves are exponential and 
show the same mean life within the limits error. Values the diffusion con- 
stants obtained from these curves are consistent with the results previous 
workers both fields. higher pressures the metastable decay remains expon- 
ential while the reciprocal the electron population changes linearly with 
time. general, the two decays are unrelated. Only the highest pressure 
studied, and this case for only the first half millisecond, was measurable 
effect positive ion recombination found the metastable population curve. 


INTRODUCTION 


Accurate measurements the decay the populations the levels 
the configuration neon were made the University Toronto 
Eckl (1951, 1953) following cessation d-c. discharge the gas 
low pressure. The decay the population the state was found 
exponential agreement with existing theory. the other hand was 
shown that the decay the population the *P» state could represented 
the sum two exponential terms. The population the non-metastable 
state was found measurable and behave similar manner that 
the state, while that the non-metastable state was too small 
followed accurately. form interchange populations between the 
four levels seemed offer satisfactory explanation for the double exponential 
decay curves the and states. 

seems possible that the decay curves the populations these states 
might affected either the lifetime imprisoned resonance radiation, 
the presence and recombination positive ions and electrons. The 
present work was undertaken investigate the latter these, particular 
measure under identical conditions the same gas the decay the popu- 
lation positive ions and the metastable states. These measurements are 
made possible the recent development techniques for the accurate 
determination free electron populations selected times after the cessation 
the discharge producing them. These methods are described Biondi 
(1951) and others. The measurements imprisoned resonance radiation 
require the development special apparatus. This development progress 
but results have not yet been obtained. 


1Manuscript received original form September 1956, and, revised, May 13, 1957. 
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METHOD 


Under suitable conditions (see, for example, Eckl 1953) metastable atoms 
revert lower either diffusion the walls and electrodes 
collision with another particle the gas. The rate the first these processes 
inversely proportional pressure, while that the second directly 
proportional pressure. Three-body collisions one metastable atom with 
two normal atoms would occur rate proportional the square the 
pressure. Such process has not been found occur appreciable rate 
inert gases room temperature (Phelps and Molnar 1953). Theory based 
the two processes reversion above, largely due Zemansky (1929), 
predicts that the decay population metastable state should exponen- 
tial and the form 


where the density metastable atoms, the initial density, the time, 
and the mean life particle the metastable state. given the 
relation 


where the diffusion coefficient metastable atom the normal gas 
unit pressure, the frequency destruction metastable atoms 
two-body collisions unit pressure, and the characteristic diffusion 
length the container. the case cylindrical container, given 


To 


where and are the length and radius the cylindrical container. 

The decay the population metastable level determined measuring 
the absorption, shown lines the neon spectrum ending that level, 
accurately determined intervals following the cessation the exciting dis- 
charge. For such line incident intensity and transmitted intensity 
can shown from the theory above that 


This form valid under certain conditions fully discussed Eckl (1951) 
and adequately fulfilled the experimental apparatus used this work.* 

When the electrical excitation the discharge tube ceases there remains, 
addition the metastable atoms, population positive ions and electrons 
distributed throughout the unexcited gas. Because the relatively enormous 
Coulomb forces these ions and electrons form plasma essentially zero 
net charge per unit volume. The diffusion the charged entities subject 
the additional condition over-all neutrality referred ambipolar 
diffusion. 


*The validity the above expression for all absorption levels and other conditions the 
experiment has been discussed length Eckl (1951). Much this discussion was removed 
from his published paper (1953) for lack space. The authors not feel that repetition 
this discussion appropriate. 
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The method measurement electron and hence positive ion densities 
makes use the alteration the resonance frequency microwave cavity 
enclosing the plasma under investigation. Slater (1946, 1950) has shown that 

where the angular resonant frequency the cavity the absence 
free electrons, the change angular resonant frequency due the 
presence free electrons, the number free electrons per unit volume 
function spatial coordinates, the electric field strength the 
resonant mode used, the permittivity free space, and are the 
electronic charge and mass, and and are the volume the cavity and 
the discharge tube respectively. This relation shown valid provided 
that wo. The integrals may evaluated the spatial distribution 
the electric field and electron density are known. For the (or 
mode cylindrical cavity radius 


Thus the field entirely axial and uniform along the length the cavity. 
Solutions for terms the population for such cavity for 


distributions characteristic both diffusion decay and recombination decay 
are indicated Biondi (1951). 


EXPERIMENTAL 


The method and apparatus for the measurement the population atoms 
metastable states has been described Eckl (1953). Only brief outline 
included here. simplified block diagram the apparatus shown Fig. 
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Fic. Simplified block diagram apparatus for metastable state measurements. 
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All triggering and timing intervals are fixed stable audio-oscillator 
operating 2000 cycles per second. Low pressure gas the absorption tube 
excited sharp rectangular pulse 500 microseconds’ duration, having 
repetition frequency which submultipleof the original oscillator frequency. 
any integral multiple 500 microseconds after excitation has ceased the 
emission tube may fired with pulse microseconds’ duration. The 
light from the emission tube collimated and passes through the absorption 
tube where some the light absorbed, raising metastable atoms 
higher levels. The absorption tube diaphragmed limit the light beam 
the central part the tube. The degree absorption dependent the 
population the metastable level the time the emission tube flash. 

After passing through the absorption tube, the light focused the slit 
monochromator which focuses single spectrum line photomultiplier 
tube. The resulting voltage pulse applied through stepped calibrated 
attenuator one input differential amplifier. the second input 
this amplifier similar voltage pulse applied, this pulse being produced 
light from the emission tube which passes through red filter directly 
second photomultiplier tube. The filter serves limit the wavelength response 
the reference channel the same region used the absorption channel. 
After several stages amplification and gating reduce noise the two pulses 
are applied the grids two carefully matched pentode amplifiers each 
with large integrating condenser its plate circuit. Thus the output the 
differential amplifier proportional the difference amplitude the 
two input pulses averaged over many cycles reduce the effect random 
fluctuations the number electronsemitted from the photomultiplier surfaces 
during any pulse. The d-c. output voltage, instead being read galvano- 
meter was done Eckl (1953), applied d-c. negative feedback 
galvanometer amplifier whose output indicated Type Sensitive 
Research microammeter. This system thoroughly stabilized the negative 
feedback and can read much more quickly and easily. The circuit used 
essentially that described Montgomery (1947) with few modifications. 

taking readings large attenuation inserted the absorption tube 
channel means the calibrated attenuator. With the absorption tube 
excitation turned off, the two-channel amplifier balanced uncalibrated 
gain controls give zero output. The absorption tube excitation then 
turned and the balance restored nearly possible reducing the 
attenuator setting the absorption tube channel. Departure from exact 
balance read the output meter and converted decibels attenuation 
interpolation between steps the attenuator. Thus values 
are read directly three figures from the attenuator setting and meter 
deflection. Repetition this process for number time delays the 
emission tube pulse after the end the absorption tube excitation yields 
complete decay curve population against time. 

The apparatus used the electron density measurements shown 
simplified form Fig. The pyrex discharge tube, referred the absorption 
tube the metastable measurements above, contained wholly within the 
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Trigger out 
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Coupler 


Modulator 


7078 
Klystron Experimental 
Generator 


Simplified diagram apparatus for electron density measurements. 


experimental microwave cavity. The gas this tube excited the same 
500 microsecond d-c. voltage pulse described above. The delayed pulse, 
used fire the emission tube for the metastable state measurements, triggers 
saw-tooth wave which swiftly varies the frequency low level probing 
microwave signal voltage through the expected resonance frequency the 
cavity. This signal coupled into the cavity through piece RG-21/U 
attenuating cable and matching device. The reflected wave from the cavity, 
acting termination the transmission line, sampled the directional 
coupler, and then detected, amplified, and applied the vertical plates 
triggered oscilloscope. Because the coupling loop the cavity adjusted 
present very nearly matching load the resonance frequency the 
cavity, the amplitude the reflected wave decreases the microwave 
frequency approaches resonance, reaching minimum value when the applied 
frequency coincides with the resonance frequency the cavity containing 
population electrons. 

determine the frequency which this minimum reflection occurs, 
second resonance minimum from calibrated reference wavemeter cavity, 
also coupled the system shown Fig. superimposed the oscillo- 
scope display. This second minimum brought into coincidence with the 
first tuning the wavemeter cavity while watching the oscilloscope display. 
The frequency resonance determined from the wavemeter cavity setting. 
time marker generator used measure the small delay time between 
the initial trigger pulse for the saw-tooth and the occurrence the resonance 


Trigg 
m 
tL, 
Junction 
Ref. Wavemeter 
Detector 
Preamp 
arker 
Oscilloscope 
Trigger 


946 CANADIAN JOURNAL PHYSICS, 35, 1957 


minimum. Since the delay the trigger pulse known, the time after cessation 
the discharge which resonance occurs the known frequency may 
accurately determined. The resonant frequency the cavity with discharge 
the tube found similar manner, the frequency difference giving the 
electron density. Repeating these measurements for series delay times 
the trigger pulse permits the plotting complete decay curve for the 
electron (or positive ion) population. 


APPARATUS 


The microwave generator 707B reflex klystron with electronically 
regulated voltage supplies. The frequency modulated applying saw- 
tooth wave the reflector electrode the usual way. The circuits used are 
conventional. The reference wavemeter reflection type resonant cavity 
having excess 25,000. was calibrated against Marconi Crystal 
Calibrator. The detector uses 1N23B crystal diode. This followed 
12AY7 preamplifier and cathode follower output feeding commercial 
amplifiers adequate gain and bandwidth. The directional coupler 
commercial unit having nominal coupling —30db. and minimum 
directivity ratio over the frequency range used. 

The microwave cavity and pyrex discharge tube are shown detail 
Fig. The cavity was designed resonant the mode 


Vacuum 
System 


Fic. Cross-sectional view experimental cavity and pyrex discharge tube. 


the neighborhood 3000 megacycles per second. this frequency coaxial 
line can used for short microwave connections without too much loss. 
Fixing the resonance frequency fixes the diameter the cavity, operation 
this mode assumed, but the length still arbitrary. The length actually 
chosen compromise value permitting long absorption paths for the meta- 
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stable measurements while still preventing excitation higher and undesirable 
modes. aid suppressing these undesirable modes, the feeding loop 
placed the center the length the cavity. 

The cavity was constructed that could easily dismantled permit 
outgassing the electrodes the discharge tube. The use split brass 
cylinder with carefully fitted overlapping end flanges allowed this 
done. Two rigid clamps along the length the cavity were added prevent 
warping. permit metastable state measurements made under identical 
conditions those for electron density measurements, windows ports allow 
the passage light from end end. These were made 1.5 cm. internal 
diameter and cm. long. The effect that window surrounded 
waveguide beyond cutoff. Microwave energy loss reduced negligible 
value. 

The pumping line and one electrode are brought out through hole centered 
the split the cavity one end, and the other electrode through hole 
symmetrically placed the opposite end. These openings are surrounded 
halved cylinders which are clamped together assembly give additional 
rigidity. Sleeves can used adjust the length these cylinders minimize 
the loss microwave energy along the concentric line formed the electrode 
support and cylinder. The probe signal coupled into the cavity loop 
coupling the magnetic field the central plane the cavity. 

spite the sources loss introduced the necessity using de- 
mountable cavity, the was about 4000 without the discharge tube. This 
fell considerably with the tube mounted. large part the additional losses 
appeared due the tungsten electrodes necessary for d-c. excitation 
the discharge. These electrodes were circular rings tungsten wire mounted 
nearly possible plane perpendicular the axis the cavity and 
the electric field the resonant mode. Small departures from exact positioning 
this plane allow small coupling the microwave field and permit energy 
escape along the electrode connections. This escape was minimized 
above. 

major difficulty arose owing sputtering metal from the electrodes 
the exciting d-c. discharge. These sputtered metal films formed conducting 
paths parallel the microwave electric field and when certain critical 
thickness was reached caused the fall quite suddenly values which 
the resonance minimum became broad that its frequency could not 
determined with sufficient accuracy. The discharge tube had then 
removed from the cavity for cleaning. spite this difficulty, was con- 
sidered advisable retain the d-c. excitation the gas that the metastable 
state results would strictly comparable those Eckl (1953) and others. 

Treating the discharge tube cylinder the dimensions shown 
Fig. the characteristic diffusion length 0.697 cm. preparing for 
operation the discharge tube was pumped down and maintained pressures 
registering less than mm. ionization gauge. The elec- 
trodes were outgassed bright red heat induction heater and the walls 
the glass system were flamed for long periods. The neon used was spectro- 
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scopically pure and furnished liter pyrex flasks. Impurities remaining 
introduced the handling process were removed exposing the gas 
calcium vapor. With these procedures little trouble was experienced with 


contamination. 


RESULTS AND DISCUSSION 


Decay curves for both electron and metastable density were obtained for 
the pressure range from 0.5 12.0 mm. Hg. previously observed and 
predicted theory, the decay the population the metastable state 
exponential all pressures. slight initial departure from linearity 
indicated the lower pressure curves, Fig. the first point showing 
consistent tendency high. This was also observed Eckl (1953). Since 


0.2 
0.1 


0.05 


Fic. Representative curves for metastable state decay. 


seemed possible enhance the effect using very slightly impure neon, 
believed that this may due the presence light impurity, possibly 
hydrogen, which rapidly adsorbed after the discharge ceases. This small 
effect entirely different order magnitude from the double exponential 
decay shown the level Eckl’s observations and repeatedly checked 
during the early part this work. the highest pressure, however, the 
first point low shown Fig. This effect also was observed Eckl 
(1951). This increase population the level high pressures the 
first half millisecond after the discharge ceases certainly real; will 
discussed below. 

The decay the metastable state was chosen for the comparison with 
the electron density decay because the absorption the spectral lines ending 
the state sufficiently large that measurements are possible much 
lower excitations than the case with the state. this work the excitation 
had kept rather low avoid rapid sputtering the walls the tube 
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with consequent lowering the cavity Eckl (1953) has established the 
relationship between the decay the populations the metastable and 
states and this relationship was checked repeatedly for various excitation 
levels the early part this work. Hence was judged sufficient establish 
the relationship the electron population decay the decay population 
the single metastable state. 

The variation the mean life the metastable state with pressure 
shown Fig. The curve described the equation 


where seconds and mm. Hg. This may compared with 


Mean Life sec. 


Metastable Decay 
Electron Decay 


Pressure 


Variation mean life with pressure. 


derived Zemansky (1929). this expression the first term represents 
destruction the metastable state diffusion the walls and the second 
represents destruction the state collision with normal atom. 

The constant term the experimental expression was suggested Eckl 
(1951, 1953) one number empirical attempts represent his results 
for the pressure dependence the mean life this state. Inspection the 
results Phelps and Molnar (1953) indicates that their experimental values 
would fit curve with such constant term somewhat better than the fit 
shown for the theoretical form the equation. The added constant term 
would seem indicate that appreciable amount spontaneous reversion, 
independent pressure, occurs from the metastable the normal state. 
Theoretical estimates the transition probability indicate that such 
reversion should negligibly small but these estimates are admittedly crude. 
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Reversion due external illumination found Meissner and Graffunder 
(1928) inadequate explain the effect because the large intensities found 
these workers necessary cause appreciable effect. the present 
work, the absorption tube was enclosed the microwave cavity, reducing 
ambient light minimum. Since the results the present work and 
Eckl (1953) and Phelps and Molnar (1953) agree representing the decay 
the population this state any one pressure straight line logarithmic 
relation, would seem unlikely that experimental error could take this form. 

This deviation from the simple theory may also indicated the intro- 
duction non-integral powers pressure into the expression suggested 
Eckl (1953) and Phelps (1955). 

The coefficient the term directly proportional the pressure shows that 
the frequency destruction metastable atoms two-body collisions 
pressure only per second. Comparison this value with the 
frequency two-body collisions calculated for neon pressure 
shows that only about one collision effective destroying meta- 
stable atom. The use squared term the pressure does not improve the 
fit the curve the experimental points was found Phelps and Molnar 
(1953) and Eckl (1951, 1953). 

Using the value quoted above, the diffusion coefficient meta- 
stable atom normal gas mm. pressure may calculated 136 
This essential agreement with results Eckl (1951), Grant and Krumbein 
(1953), and Phelps and Molnar (1953). 

Curves Av, which proportional electron population, are plotted 
log scale against time Fig. The curves chosen correspond the meta- 
stable curves Fig. low pressures the curves electron population 
show exponential decay, indicating that the electrons and ions are recombining 
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Fic. Representative curves for electron density decay. 
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after diffusion the walls and that volume recombination negligible 
these pressures. These low pressure results could not extended very 
short delay times because the lowering the the cavity and change 
coupling the cavity the transmission line due the very high popula- 
tion free electrons. The largest frequency shift used the results was 
megacycles per second, less than one per cent the resonance frequency 
the absence free electrons, that the approximation amply 
fulfilled. 

For the low pressures for which the electron population decays exponentially, 
the mean lives are marked Fig. along with the metastable mean lives. 
The close agreement the two mean lives these pressures indicates that 
the effective diffusion coefficients for the metastable atoms and positive ions 
and free electrons approach each other within experimental error the 
pressure becomes small. This first seems somewhat surprising, but Biondi 
(1954) has shown that, given the extremely high ratio the mass the 
neon gas atoms the mass the electrons, the ambipolar diffusion coefficient 
might expected approach the normal diffusion the positive 
ion the normal gas. This coefficient turn might expected approximate 
the coefficient for the metastable atom. 

pressures mm. and higher the logarithmic plot electron density 
against time starts show definite curvature which becomes more pronounced 
for higher pressures. the highest pressure, mm. Hg, the log plot (Fig. 
again appears nearly linear out time milliseconds. This appearance, 
however, due the fact that population itself changes only small 
factor the range measurement. intermediate pressures, for which the 
curves are not shown, curvature the log plot very marked. Comparison 
the two log plots mm. for metastable states Fig. and for electron— 
ion population Fig. shows that relation between the two decays exists 
beyond time millisecond. 

For all pressures 5mm. and higher which were studied, was found 
that plot the reciprocal against time yielded straight line. This 
illustrated Fig. and characteristic volume recombination. For this 
process the appropriate differential equation 


which the electron (or positive ion) density and the recombination 
coefficient. Solving this and evaluating the derivative gives 
=at and 


This shows that the rate recombination very large when small and 
rapidly falls off increases. Since least part these recombinations 
the body the gas may expected result radiation cascade, such 
recombinations could act increase the population the metastable state. 
This effect should very large small delay times and decrease very rapidly 
This effect adequate explain the initial increase population 
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Representative curve reciprocal electron density against time high pressure. 


the metastable state mm. pressure shown Fig. Unfortu 
nately, was pointed out above, was impossible carry the determination 
electron population back zero delay time under the same conditions 
because. the very large population electrons. 

The initial concentrations (i.e. the time the exciting discharge ceases) 
both for electrons and positive ions, and for metastable states, may found 
extending the appropriate straight line plots back the zero time inter- 
cepts. Expressions for these initial populations are given Eckl (1953) and 
Biondi (1951) for the case metastable states and for electrons respectively. 
For the discharge conditions used this work, the two initial populations 
are the same order, each being about one 10’ the atoms present. 

Recombination coefficients may also calculated for the cases which 
this process the predominant cause disappearance free electrons, that 
for pressures above Hg. From the slopes the linear plots the 
reciprocal the electron population against time, found have values 
5.16, 8.0, and 12.0 mm. pressure respectively. The values are 
the same order those obtained Biondi and Brown (1949) and Holt, 
Richardson, Howland, and McClure (1950). Since the pressure range this 
work does not overlap that the previous workers, exact comparison cannot 
made. 


CONCLUSIONS 


all pressures the range studied the decay the population the 
metastable state neon exponential. low pressures the atoms are 
removed from the metastable state principally diffusion the walls. 
high pressures removal from the metastable state due principally collisions 
the metastable atom with unexcited atom the body the gas. The 
efficiency such collisions low, the order one spontaneous 
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reversion term seems necessary explain why the experimental values 
mean life near its maximum for pressure variation fall consistently below 
the theoretical curve based the processes above. 

low pressures the population electrons and positive ions decays 
exponentially ambipolar diffusion give recombination the walls 
the vessel. The ambipolar diffusion coefficient for this process is, within the 
limits experimental error, the same that for the metastable atoms. 
higher pressures the decay electron population falls off that the reciprocal 
the population plotted against time gives straight line, this being charac- 
teristic volume electron-ion recombination. this region pressure 
relation exists between the decay the populations positive ions and 
metastable states, except the first half millisecond the highest pressure 
used. Here there evidence that the population the metastable state 
increased the recombination positive ions and electrons. Recombination 
coefficients have been evaluated the higher pressures. 

There indication that the presence positive ions after termination 
the discharge can offer explanation the double exponential decay 
found for the metastable state. 
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INFRARED INTENSITIES AND BOND MOMENTS 
SULPHUR DIOXIDE! 


ABSTRACT 


New values have been obtained for the intensities the three fundamentals 
units The normal coordinate calculation gives values 
for the bond moment derivatives and for the molecular dipole moment which 
demonstrate that the amount charge moving all three vibrations greater 
than that inferred from the equilibrium charge distribution. The involvement 
unshared electrons suggested. 


INTRODUCTION 


Integrated absorption coefficients for the infrared bands sulphur dioxide 
vapor have been published Eggers, Hisatsune, and Van Alten (1955) and 
Morcillo and Herranz (1956, 1957). The measurements reported 
here were made independently the period between these two publications. 
They are presented resolve the disagreement between the two previously 
published sets, and revise the interpretation these intensities terms 
dipole moments and moment derivatives for the sulphur dioxide molecule. 

has been customary attempt relate the observed infrared intensities 
the equilibrium distribution charge the molecule through the assump- 
tion that the formal atomic charges the equilibrium configuration remain 
fixed their respective atoms these atoms move molecular vibrations. 
Tests the applicability this concept have been made difficult the 
experimental problems involved obtaining accurate intensity data and 
the scarcity data adequate accuracy for the potential constants even 
the simpler polyatomic molecules (Callomon, McKean, and Thompson 1951; 
Hornig and McKean 1955). Thus, despite the fact that numerous studies 
have already demonstrated the failure the fixed-charge model (Hornig and 
McKean 1955) worth while provide accurate intensity values and 
moment derivatives for the sulphur dioxide molecule. This because sulphur 
dioxide presents almost unique combination convenience measurement 
with the availability accurate molecular geometry and potential constants 
(Crable and Smith 1951; Kivelson 1954; Polo and Wilson 1954), and should 
therefore useful the quantitative testing more detailed theories 
vibrational intensities. further attractive feature that all three phases 
sulphur dioxide may conveniently studied examine the influence 
intermolecular forces the internal charge distribution; such studies will 
form the basis subsequent paper. 

The paper Eggers, Hisatsune, and Van Alten (1955) contains all the 
relevant procedure for the reduction the observed intensities dipole 
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moment derivatives. Unfortunately few errors were inadvertently introduced 
into their publication (Hisatsune 1956). The normal coordinate treatment 
has been repeated using their notation and symmetry coordinates, that 
(1) So = 
with the result that the transformation from the symmetry coordinates 
the normal coordinates takes the form 
(2) 
where the numerical coefficients the are units The definitions 
the symmetry coordinates together with the assumptions 
charge model vield the relations 


These relations differ from those Eggers, Hisatsune, and Van Alten (1955) 
the signs the coefficients, and factor the equation for The 
latter error was called our attention Polo. 


EXPERIMENTAL 

Intensities were determined the method Wilson and Wells (1946) 
total pressures and atmospheres, using nitrogen the broadening 
gas. Partial pressures sulphur dioxide between 0.8 mm. and mm. were 
used with 10.0cm. path length. Glass-bodied cells were used avoid 
adsorption effects noted with metal cells. The effect incomplete mixing 
the two gases was investigated, and effect was found the pressures 
used. The influence total pressure the apparent integrated absorption 
coefficient was investigated for several partial pressures and atmos- 
phere total pressure was found give complete broadening over most the 
range partial pressures employed. 

The nitrogen used was dried slow passage through liquid-oxygen-cooled 
traps followed reheating room temperature. This procedure avoided 
the loss observed when premium grade nitrogen was used without 
drying. Sulphur dioxide from the Matheson Chemical Co. (of stated purity 
was found condense yellowish liquid, and contain few per 
cent some gas not condensable liquid oxygen temperature. was there- 
fore subjected one two vacuum fractionations prior being evaporated 
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into the evacuated cell. singly distilled sample was analyzed mass spec- 
trometry and shown not less than 99% pure, and have detectable 
sulphur-containing impurities. The amount sulphur dioxide the cell was 
determined using differential oil manometer, and was confirmed quantita- 
tive chemical analysis. The analysis involved complete slow passage the 
cell contents through liquid air traps series, conversion the trapped 
sulphur dioxide sulphate with ammoniacal hydrogen peroxide, and gravi- 
metric determination the sulphate barium sulphate (Vogel 1951). The 
chemical analyses showed scatter +5%, but agreed with the manometric 
pressure within 2%. 

Spectra were recorded the NaCl and KBr regions Perkin-Elmer 
Model double-beam infrared spectrophotometer with matched cell con- 
taining nitrogen the reference beam. Spectral slit widths were chosen 
less than 0.65 the half-width one rotational wing order minimize 
the extrapolation slope (Ramsay 1952) while permitting reduction the 
noise level 0.2% transmission with scanning speed per minute. 
These conditions high performance could not achieved the low frequency 
tail the band which lies the extreme the range the KBr prism. 
This tail was therefore estimated symmetry from the high-frequency tail; 
the area involved was about the total band area. Band integrations 
were performed using Simpson’s Rule with sufficient number intervals 
keep the sums the even and odd ordinates agreement within 0.5% 
for all but the weakest bands. Background and stray light corrections were 
applied required. 

The integrated absorption coefficients were evaluated both from the limiting 
slope the origin plot band area against partial pressure and from the 
intercept plot the apparent integrated absorption coefficients against 
partial pressure. The latter plots are shown Fig. felt that the intensities 
are accurate within the uncertainty the chemical analysis this 
being the largest source systematic error encountered. The over-all repro- 
ducibility, taken the standard deviation points about the least squares 
line the least favorable case, was +4.5%. 


RESULTS 
Table presents comparison the intensities obtained this investigation 
with the previously published values. The agreement between this work and 


that Morcillo and Herranz essentially within the limits accuracy, while 
the results Eggers, Hisatsune, and Van Alten are roughly 50% low. This 


TABLE 


Band Ref. (4) Ref. (9) This work this work 
1151 226 350+20 317+8 58.1 
519 257 360+30 376+20 63.6 
1361 1606 2560+70 +165.2 


INTEGRATED ABSORPTION COEFFICIENTS FOR SULPHUR DIOXIDE 
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(litres 


Partial Pressure 
2 


Fic. Apparent integrated absorption coefficients the sulphur dioxide fundamentals 
partial pressure, 25° and total pressure 760 mm. 


discrepancy has been discussed with Professor Eggers, who has conceded that 
some unsuspected bias, such included nitrogen their sulphur dioxide 
sample, may have prejudiced their determination partial pressures mass 
spectrometry. any event felt that, terms the errors which 
the methods are susceptible and view the two agreeing sets, the higher 
values must considered correct within their limits error. 

The derivatives the total molecular dipole moment with respect 
distortions the manners specified the symmetry coordinates S;, 
are presented Table with the bond moment u,, and its derivatives 
calculated from the relations (3) above. 

The fact that the observed intensities are proportional the square the 
corresponding gives rise alternative values for and 
depending whether the signs and are the same different. 
these signs are the same, the predicted total molecular dipole moment for 
becomes +3.14 Debyes they are different. These 
values are compared with the static dipole moment 1.59 Debyes 
obtained from the pure rotational spectrum (Crable and Smith 1951) from 
dielectric constant measurements (Brown 1956). 


DISCUSSION 
The disagreement between the ‘‘dipole moment” derived from the infrared 
intensities and the dipole measured the usual way 
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TABLE 
MOMENTS CALCULATED FOR SULPHUR DIOXIDE 
(Values are Debyes; all others are 


Relative signs 


Moment and Calculated value 

Same +3.62 
Different +1.69 

Same +1.80 
Different 

+6.44 
Same 
(from and +2.38 
Same 
(from and Different +3.12 
+5.27 


(from v3) 


much greater than any the errors involved measurement, that 
quite clear that these two quantities are not the same physically. 

The failure the fixed-charge model further emphasized comparison 
its other predictions with the results Table II. molecular orbital 
treatment sulphur dioxide Moffitt (1950) assigns formal charges 
+0.48 electrons the sulphur and —0.24 electrons each oxygen. 
This charge distribution consistent with the observed static dipole moment. 
Thus the assumptions the fixed-charge model should find not only 
agreement the values, but should find 0.24 electrons 
about e.s.u. independent the symmetry class the normal 
vibration. None the alternatives Table are good agreement with 
this prediction. The signs the derivatives may also considered. 
formal charges give the positive sign for the coordinate systems used 
here, and predict the positive sign for and the negative sign for 
that the signs and should also taken being different. 
This choice, however, leads large disagreement the value 
from and with that from studies other molecules conclusions 
signs and magnitudes derivatives have occasionally been based 
the consistency argument that the magnitudes calculated from vibrations 
different symmetry ought agree. The accumulated evidence presented 
Hornig and McKean (1955) calls this procedure into question; the present 
evidence demonstrates that the procedure inconsistent with the hypothesis 
used the calculations. 

The over-all conclusion from these findings that the amount charge 
redistributed the molecular distortions sulphur dioxide much greater 
than the formal charges the equilibrium distribution. This suggests that the 
concepts rehybridization, orbital following, and the contributions unshared 
electrons which Hornig and McKean have applied successfully and 
HCN are needed for discussion the intensities. 
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has been shown that unshared pairs can make large contributions the 
static moments (of the order Debye unit), and that they are very easily 
polarized (Smyth 1955). The molecular orbital calculations Moffitt for 
sulphur dioxide and related compounds show that the sulphur unshared 
electrons occupy orbitals which are directed approximately tetrahedrally with 
respect the S—O bonds, and are pulled towards the oxygen atoms. Thus 
the unshared pair moment the equilibrium configuration opposes the S—O 
bond moment. found that the configuration the sulphones, 
nearly tetrahedral while the S—O bond order almost unaltered from 
its value and the dipole moment increases from 1.59 Debyes sulphur 
dioxide more than Debyes for dialkyl sulphones. 

assume that the lone pair orbitals follow the motions the oxygen 
atoms but otherwise retain the concepts the fixed-charge model, qualitative 
consistency can given the signs and relative magnitudes the moment 
derivatives sulphur dioxide. Thus the lone pair moment, ordinarily 
opposition the bond moment but reduced its effect displacement 
the lone pair orbitals towards the oxygens, would further decrease the 
S—O bonds stretched slightly, leading augmented positive change 
for distortions type Similarly, the opposition the lone pair moment 
the bond moment would increase for increases the bond angle, leading 
augmented negative value Distortions according would 
cause both bond and lone pair moment changes the same direction, leading 
large positive values 

Although the above discussion does not advance our understanding 
infrared intensities, does call attention directions for further work. 
Moreover the provision well-established set intensities and moment 
derivatives for simple polyatomic molecule should valuable the testing 
more detailed theories molecular vibrations. The sensitivity the infrared 
intensities the details the internal charge distribution suggests their use 
the investigation intermolecular forces, since even weak perturbations 
may expected produce observable changes intensity. 
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QUANTUM MECHANICS WITH INDEFINITE 
Gupta? 


ABSTRACT 


The formalism quantum mechanics vector space with indefinite 
metric developed simplified way, and the properties the indefinite 
vector space are discussed detail with particular attention vectors neg- 
ative and vanishing norms. shown how the present formalism simplifies 


the use the indefinite metric physical applications. 


INTRODUCTION 


Usually quantum mechanics use vector space with positive definite 
metric, because well known that indefinite metric gives rise the 
difficulty negative probabilities (Dirac 1942; Pauli 1943). But has been 
shown recent years that vector space with indefinite metric can also 
used quantum mechanics without any difficulty, provided that 
ensured means supplementary conditions that the states negative 
probabilities are unobservable. this way has been found possible 
carry out the quantization the electromagnetic field 
way (Gupta 1950). indefinite metric has also been used carry out the 
quantization the vector meson field (Gupta 1951) and the gravitational 
field (Gupta 1952), and introduce auxiliary fields quantum electro- 
dynamics consistent way (Gupta 1953). 

the above applications the indefinite metric two sets quantities 
were introduced (called the Hermitian conjugates and the adjoints), which 
somewhat destroyed the elegance the treatment and also made unneces- 
sarily complicated. Therefore, shall here develop the formalism quantum 
mechanics with indefinite metric simplified way, which involves only 
one set quantities. fact, the present formalism identical appearance 
with Dirac’s formalism (Dirac 1947) the usual quantum mechanics 
long keep our treatment abstract form and not make use 
the matrix representation operators. this paper shall also discuss 
some detail the properties the vector space with indefinite metric 
and its physical interpretation. shall pay special attention vectors 
vanishing norms, because the role such vectors quantum mechanics 
has not been fully before. Further, will shown how the present 
formalism enables simplify the use the indefinite metric physical 
applications. 

shall refer vector space with indefinite metric the indefinite 
vector space, while accordance with the usual practice vector space with 
positive definite metric will called the Hilbert space. 
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THE INDEFINITE VECTOR SPACE 
Let consider vector space infinite dimensions, which the com- 
ponents vectors are general complex numbers. Let call this space 
the indefinite vector space, and let denote any vector this space 
the symbol With the indefinite vector space let associate dual 
space such that corresponding any vector there always exists vector 
the dual space. Let assume that the relationship between the two 
spaces such that and are any two vectors the indefinite vector 
space, then 


where number, and denotes its complex conjugate. follows from 
(2) that 
(3) real number. 


The only difference between the above assumptions and those the usual 
treatment that now the quantity not positive definite, but according 
(3) can positive, negative, zero. 

Since the present treatment very similar the usual one, will 
convenient retain the usual notation and terminology for the various 
quantities. Thus, the vectors and will called the Hermitian conjugates 
each other, and the real number will called the norm the vector 
But, now any vector satisfying the more general condition 


(4) 


will said normalized, because cannot change the sign the norm 
vector multiplying with any number. Further, two vectors and 
are such that their product vanishes, they will said orthogonal. 
set vectors are normalized well orthogonal each other, they 
will said form orthonormal set. 

can also define operators the indefinite vector space quantities 
which operating any vector change into some other vector the same 
space. Thus, operator, and are vectors the indefinite 
vector space and the dual space respectively, and therefore quantity 
the form number. operator will called the Hermitian 
conjugate for any arbitrary vector have 


(5) 


whence follows that 


can easily shown that the Hermitian conjugate the complex conjugate 
product operators and vectors obtained reversing the order 
all the factors and then replacing each factor its Hermitian conjugate. 
operator such that 
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will called Hermitian operator. 

important note that the self-adjoint operators the earlier treat- 
ments the indefinite metric correspond the Hermitian operators the 
present treatment, while the Hermitian operators the earlier treatments 
have analogue the present treatment. 


BASIC POSTULATES QUANTUM MECHANICS WITH 
INDEFINITE METRIC 

The basic postulates mechanics with indefinite metric can 
stated follows: 

(1) The dynamical state system any time can represented 
vector indefinite vector space. 

(2) The dynamical variables for system can represented suitable 
operators, which operate the state vectors the indefinite vector space. 

(3) The result single measurement dynamical variable for system 
cannot general predicted with certainty. But, the measurement 
dynamical variable for system certain state carried out large number 
times, the average result measurement the so-called expectation value 
always given 
(8) (a) 


where the operator represents the dynamical variable, and normalized 
vector representing the state the system. 

Now, represents any quantity physical its expectation value 
must real, i.e. must have 


Since the above relation must hold for any arbitrary vector follows 
that Hence, all quantities physical interest are represented 
Hermitian operators also the present treatment. 


PROPERTIES HERMITIAN OPERATORS 


Although evident from (6) and (8) that the Hermitian operators have 
real expectation values even the present treatment, some properties the 
Hermitian operators are now slightly different from those the usual treat- 
ment. 

the eigenvalue that 


Taking the complex conjugate the above equation, and remembering that 

(11) 


From (10) and (11), then get 
(12) 
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Since according the present treatment the norm may vanish, does 
not necessarily follow from the above equation that Thus, the 
present treatment all the eigenvalues Hermitian operator need not 
real. However, can still infer from (12) that all the eigenvalues 
Hermitian operator corresponding eigenvectors non-vanishing norms 
are real. 

Similarly, cannot now establish general the orthogonality all the 
eigenvectors Hermitian operator. But, can show the usual way 
that all those eigenvectors Hermitian operator which have non-vanishing 
norms form orthogonal set vectors the indefinite vector space, pro- 
vided that all the eigenvectors corresponding the same eigenvalue are 
chosen suitable way. 

shall now assume that all the eigenvectors Hermitian operator 
have non-vanishing norms, the orthogonal set eigenvectors formed 
them complete set, i.e. any vector the indefinite vector space can 
expressed linearly terms these orthogonal vectors. Thus, the eigen- 
values Hermitian operator are discrete and all its eigenvectors have 
non-vanishing norms, can express normalized vector 


(13) 


where the are numbers, and the are orthogonal eigenvectors corre- 
sponding the eigenvalues a,. When these eigenvectors are also normalized, 
have 

where the plus and the minus signs correspond the eigenvectors positive 
and negative norms respectively. Now, let any arbitrary function 


(15) f(a), = f(an)¥n, 
and 


which gives, using (14), for the average result measurement 
the state 


But, the average result measurement f(a) may written 


m 


where some value for other than the eigenvalues a,, and and 
are the probabilities obtaining and a,,’ respectively the result 
measurement Since the function arbitrary, find com- 
paring (17) and (18) 


n 

m n m n 
n 
n = 
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which shows that vector with negative norm represents state occurring 
with probability. 


REPRESENTATION OPERATORS MATRICES 
norms, can represent any operator the frame reference the 
vectors matrix, whose elements are given 
(20) 


also follows from (20) and (6) that Hermitian, 
(21) 


which shows that Hermitian operator represented Hermitian matrix 
also the present treatment. But, shall now obtain important result, 
which slightly different from the corresponding result the usual treat- 
ment. 


where, already mentioned, the plus and the minus signs correspond 
vectors positive and negative norms respectively. Then, any arbitrary 
vector can expressed terms the expansion the form 
(13), and have 


But, according (22), 

24) 

that (23) reduces 

Since the vector arbitrary, follows that 

Hence, have 


UNITARY TRANSFORMATIONS 
usual, shall call any operator that satisfies the relation 
unitary operator. carry out transformation from set vectors 
and operators another set according the relations 
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where and W’, denote the corresvonding quantities 
the two sets, and unitary operator, then the relations (30) will 
said describe unitary transformation. 

obvious that all the usual properties unitary transformation 
remain unchanged the present treatment. For instance, can easily 
shown that pair vectors operators are Hermitian conjugates 
each other, they preserve this property under unitary transformation. More- 
over, any algebraic relation involving vectors and operators remains unchanged 
under unitary transformation, and can pass over the usual way from 
one representation system another means such transformation. 

is, however, interesting note that the Hermitian operators the 
present treatment correspond non-Hermitian operators the earlier 
treatments the indefinite metric, and therefore transformation the 
form (30) corresponds non-unitary transformation the earlier treat- 
ments. 


REMARKS THE USE INDEFINITE METRIC 


Owing certain difficulties physical interpretation, the formalism 
quantum mechanics with indefinite metric can used only under special 
circumstances. shall now discuss these difficulties physical interpretation 
see what limitations they impose the use the indefinite metric. 

have already come across difficulty physical interpretation 
Section where was shown that vector with negative norm represents 
state occurring with negative probability, which physically meaningless. 
Another difficulty physical interpretation connected with vectors 
vanishing norm. For, vector has vanishing norm, can multiply 
any number obtain another vector which also has vanishing 
norm. This gives rise the difficulty that although both and are nor- 
malized vectors representing the same state, they lead different expectation 
values for any dynamical variable except when this expectation value 
vanishes. 

The above difficulties physical interpretation are connected with vectors 
negative vanishing norms. Therefore, indefinite metric can used for 
quantum mechanical system without any difficulty only when the states repre- 
sented vectors negative vanishing norms the resulting treatment are 
physically unobservable. 


THE FIELD OSCILLATORS 


the quantization fields integral spin, come across two types 
which may called the normal and the abnormal 
field oscillators. The usual practice treat these field oscillators using 
vector space with known metric, which either positive definite metric 
special type indefinite metric. However, shall show that both 
types field oscillators can treated analogous way using vector 
space with unspecified metric, and then postulating the following physical 
condition: 
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Every physical system nature has state lowest energy, and the vector 
representing this state has positive norm. 


Normal Field Oscillator 
The Hamiltonian normal field oscillator the form 


(31) Ea*a, 


where positive number, and the operator and its Hermitian conjugate 
satisfy the commutation relation 
(32) aa*—a*a 


shall determine the eigenvalues and the eigenvectors this field oscillator 
using the condition stated above. 


(33 


Then, using (31) and (32), get 

which shows that and are also eigenvalues provided that 
and follows that there must exist vector such that 
(36) 


For, were not so, would follow repeated application that 
the eigenvalues have lower limit, and this would contradict our 
physical condition. 

The eigenvalue for the eigenstate Wo, satisfying (36), 


and repeated application (35) find 


are all eigenvalues corresponding the eigenvectors 


(39) Wo, a*Wo, a*"Vo, @ 


respectively. 

Further, according the physical condition, has positive norm, 
that can normalize 


Then, can express the normalized eigenvector corresponding the 
eigenvalue and its Hermitian conjugate 


as 


and the metric the vector space the system given 
(42) VV, Guns 


which positive definite. 


= 
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Abnormal Field Oscillator 
abnormal field oscillator has Hamiltonian the form 


(43) 


where positive number, and the operator and its Hermitian conjugate 
satisfy the commutation relation 
(44) aa*—a*a 


Using the physical condition, and following exactly the same procedure 
the case the normal field oscillator, can easily show that the eigen- 
values are again given 
(45) 


and the corresponding eigenvectors are 
(46) Wo, a*Wo, a*" Wo, 


where satisfies the relation 
(47) 


Further, view the physical condition, can normalize 
(48) Wo*Vo 


which enables express the normalized eigenvector corresponding 
the eigenvalue and its Hermitian conjugate 


But, owing the different commutation relation (44), the metric the 
vector space the system now given 


which indefinite. 

have seen how the present formalism enables deal with the normal 
and abnormal field oscillators analogous way. The same approach can 
also used the various physical applications the indefinite metric, 
which have mentioned Section For instance, using indefinite 
vector space with unspecified metric for the radiation field, and then 
applying our physical condition all the radiation field oscillators, easily 
arrive the indefinite metric our earlier treatment the quantization 
the radiation field. Moreover, since the physical condition Lorentz invariant, 
directly ensures the Lorentz covariance the resulting treatment. 
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THE CUMULATIVE YIELDS THE KRYPTON AND XENON 
ISOTOPES PRODUCED THE FAST NEUTRON FISSION 
Th22 1 


ABSTRACT 


The relative fission yields for mass chains ending stable krypton and xenon 
isotopes have been measured for the fast neutron-induced fission Isoto 
dilution techniques were used determine the krypton/xenon ratio assist 
obtaining the absolute fission yields for these mass chains. The absolute yields 
were determined the use two methods, both giving results which were 
excellent agreement. The fine structure observed for the curve 
compared with that heavier fissile nuclides. 


INTRODUCTION 


During the past few years considerable amount work has been done 
determine the mass distribution the products from the neutron-induced 
fission and contrast, work has been reported 
the fission yields The principal contribution date has been the 
radiochemical determination fission products Turkevich and Niday 
(1951). The study neutron-induced fission made difficult because 
the high fission threshold 1.10 Mev. (Haxby, Shoupp, Stephens, and 
Wells 1940) and the low fission cross-section, which has maximum value 
0.7 barns Mev. (Hughes and Harvey 1955). The combination these 
effects makes necessary irradiate the with high flux fast neutrons. 

determination the krypton and xenon yields for the neutron fission 
would considerable importance for better understanding the 
fission process. The krypton and xenon mass regions exhibit fine structure 
deviation from smooth curve which was first observed Thode and 
Graham (1947) for fission. This effect has been qualitatively explained 
the basis shell structure Glendenin (1949), Wiles (1950) (see also 
Wiles al. 1953), and Pappas (1953). test such postulates advan- 
tageous obtain many curves possible. The relative position 
the familiar asymmetric fission yield curve depends upon the fissile nuclide. 
Present data indicate that the heavy mass hump remains about the same 
position for all fissile nuclides. Conversely the position the light mass hump 
function the mass the fissile nuclide. The relative positions the 
krypton and xenon yields are such make the ratio xenon/krypton good 
parameter for measuring this effect. obtain such information study the 
fission yields was undertaken. 


EXPERIMENTAL 


Conditions 
The source fission spectrum neutrons was transformer the N.R.X. 
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(Chalk River) Reactor. The transformer consisted partially hollowed, 
natural uranium fuel element and occupied lattice position where the 
thermal neutron density was maximum. The cylindrical tube was in. long 
with outside diameter 1.125 in. and wall thickness 0.125 in. The 
neutron distribution within the transformer consisted the thermal and 
fission spectra. Three samples were irradiated for weeks (irradiation 
and two samples for weeks (irradiation B). 

The samples were irradiated the form thorium metal machined 
cylindrical rods 0.125 in. diameter. The thorium was sheathed cadmium 


reduce the reactions 


fission 


produced the otherwise intense thermal neutron flux. The samples were 
then placed hexagonal aluminum containers permit air-cooling during 
the irradiation. 

Analysis the thorium showed the natural uranium content less 
than This low concentration removes the possibility appreciable 
contribution from uranium fission. The thorium oxide content was about 
25% and there were traces iron possibly arising from the fabrication process. 


Extraction Procedure 

Two methods were employed extract the fission product rare gases from 
the irradiated thorium metal. The first procedure, devised Arrol, Chackett, 
and Epstein (1949) for the extraction rare gases from irradiated uranium 
metal, was found applicable thorium. The thorium was dissolved 
saturated solution cupric potassium chloride and the released gases were 
passed through drying train. The purification the rare gases was achieved 
the use calcium furnace. 

The second method consisted heating the thorium about 1000° 
and reacting with chlorine. When the reaction was completed the thorium 
chloride had all sublimed cool region the vacuum system, which en- 
sured complete release the rare gases. After this reaction the excess chlorine 
was removed hot iron wire. The final purification was achieved the 
use calcium furnace. Although both methods yielded the same results the 
chlorine reaction has the advantage minimum handling. This reduces the 
likelihood contamination with normal krypton and xenon. 


Isotope Dilution Technique 

Extraction the rare gas fission products permits mass spectrometric 
analysis the relative fission yields for the krypton and xenon isotopes. 
information can obtained concerning the quantity krypton and xenon 
the krypton/xenon ratio. obtain such information necessary 
the extracted fission gas with known quantity krypton and 
xenon that has different isotopic composition. The most desirable condition 
spike the fission gas with rare gas isotopes which are not formed 
fission cumulative isotopes. Pure samples and were 
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prepared the irradiation potassium bromide and potassium iodide. 
These isotopes are formed very low abundance primary fission products 
(Kennett and Thode 1956) and not interfere with the cumulative fission 
pattern. This spiking mixture was calibrated against standard mixture 
krypton and xenon prepared Linde Air Products. 

The spiking was done admitting known aliquot the prepared mix- 
ture during the chlorine reaction. Once the spike has been added and thoroughly 
mixed in, 100% chemical yield longer necéssary. Therefore, advan- 
tageous add the spike the reaction vessel early possible. 


Mass Spectrometry 

The fission product rare gases were analyzed with in. directional 
focusing mass spectrometer employing electron multiplier for the detector. 
The output from the multiplier was integrated and amplified vibrating 
reed electrometer and then permanently recorder. The resolving 
power the instrument was 625, which permitted doublet separation the 
hydrocarbon background from the rare gas isotopes. Analysis normal 
krypton and xenon gave results within the accepted values determined 
Nier (1950). 

Normal contamination the fission gases can detected the presence 
and change the ratio found for the spiking mixture. 
measurement these quantities correction can easily made for the 
normal contribution the observed fission pattern. 


RESULTS 
The observed mass spectrometric ratios for the krypton isotopes are given 
Table The values have been corrected for decay and branching and 


TABLE 


THE RELATIVE FISSION YIELDS THE KRYPTON ISOTOPES OCCURRING THE FAST NEUTRON 
FISSION 


Mass spectrometric ratio 


Irradiation 


Irradiation 


*Corrected for decay and 
tNormalized 10.00. 


has been normalized 10.00. The errors quoted for each sample are 
the standard deviations obtained from the measurement spectrograms. 
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The values for individual samples and the average values for the two irradia- 
tions are seen good agreement. 

Table summarizes the observed relative fission yields for xenon where 
has been normalized 10.00. With the exception the results 


TABLE 


RELATIVE FISSION YIELDS THE XENON ISOTOPES OCCURRING THE FAST NEUTRON 
FISSION 


Mass spectrometric ratio 


Irradiation 
Sample 3.00+0.01 10.00 10.70+0.03 
Sample 10.00 10.73+0.02 
Sample 10.00 10.80+0.02 
Average value 10.00 10.74+0.04 
Irradiation 
Sample 10.00 10.45+0.07 
Sample 3.04+0.02 10.00 10.50+0.03 
Average value 3.02+0.02 5.30+0.03 10.00 10.48+0.04 


*Normalized 10.00. 


from both irradiations are good agreement. Irradiation with twice the 
irradiation time will exhibit larger burn-up resulting in- 
crease the yield. The amount burn-up function the 
irradiation time and also the thermal neutron flux within the sample. 
combining the results for both irradiations the value for the thermal 
neutron flux inside the thorium was found, and proved less than 
neutrons/cm.?/sec. Calculations made with this flux indicated that 
irradiation had negligible burn-up and that the observed yield for this 
sample represents the true yield for the 136 mass chain. 

The irradiation conditions for both and are tabulated Table III. 


TABLE III 
IRRADIATION CONDITIONS FOR THE FAST NEUTRON FISSION 


Irradiation Irradiation 


Irradiation time (days) 
Integrated power (megawatt days) 
Irradiation time (days) megawatts 39.1 17.9 
Thermal flux (neutrons/cm.?/sec.) trans- 

former position for megawatts 8.20 10" 


The irradiation time has been calculated the basis full power operation. 
The value for the thermal neutron flux the transformer position has also 
been referred full power operation. 

Isotopic dilution was done for samples and irradiations and 
respectively. The results these experiments are summarized Table IV. 
this table has been again corrected for decay and branching. The 
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TABLE 
ISOTOPE DILUTION RESULTS FOR 


Nuclide Irradiation sample Irradiation sample 
8.10+0.02 
17.04+0.18 
11.79+0.10 
3.17+0.04 
11.98+0.06 
22.49+0.08 10.42+0.04 
24.23+0.12 10.92+0.05 
Total volume 65.48+0.42 30.05+0.18 
Xenon/krypton ratio 0.981+0.011 


Average value xenon/krypton ratio 0.994+0.013 


*Corrected for decay and branching. 
for neutron capture 


xenon/krypton ratio obtained finding the ratio xenon 
krypton volume’’. correction was made for the sample 
compensate for the neutron capture The agreement the xenon/ 
krypton ratio for irradiations and exceptionally good considering the 
wide difference irradiation conditions. the use the average xenon/ 
krypton ratio Table and the results given Tables and II, the abun- 
dance the krypton and xenon yields for the fast fission can 
calculated. The results this calculation are given Table 


TABLE 


THE RELATIVE FISSION YIELD ABUNDANCES 
THE RARE GAS ISOTOPES OCCURRING THE 
FAST NEUTRON FISSION 


Isotope Relative abundance (%) 
6.42+0.04 
11.74+0.08 
12.50+0.08 
9.24+0.06 


*Corrected for decay and branching. 


DISCUSSION RESULTS 


The xenon/krypton ratio not only ties together two regions the 
vield curve but also represents important parameter for determining the 
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displacement the curve function the fissile nuclide. 
examination the curves for and (Katcoff 
1957) shows that the heavy mass hump remains centered about approximately 
the same mass number while the light mass hump moves lower masses 
the fissile nuclide decreases mass. The xenon yields fall near the top the 
heavy mass hump while the krypton yields lie the low mass side the light 
hump. Any displacement the light mass hump results marked change 
the krypton yields and also the xenon/krypton ratio. The value the 
ratio changes from 4.38+0.08 for (Blades, Fleming, and Thode 1956) 
0.994+0.013 for Such change indicates mass unit displacement 
the light hump. For more quantitative interpretation this and other 
effects need, course, some knowledge absolute fission 

The problem absolute fission yields can, for this case, approached 
two ways. The first fit the relative krypton and xenon results the 
curve determined Turkevich and Niday (1951). The second 
calculate the fission spectrum flux within the transformer and thereby 
determine the number fissions the thorium samples. both methods 
give consistent yields then reasonable weight can associated with the 
calculated flux. 


Absolute Fission Yields from Curve Fitting 

The radiochemical technique has inherent errors considerable magnitude 
that Turkevich and Niday (1951) were only able obtain the general nature 
the curve. Fig. shows the fission yield curve with the 


\ 
100 120 140 


MASS NUMBER 


mass yield curve for fast neutron-induced fission The mass spectro- 
metric values are shown and the radiochemical results Turkevich and Niday (1951) 


by I. 
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mass spectrometric and radiochemical data plotted with the associated 
errors. The high precision the mass spectrometer ratios clearly shown. 
obtain the absolute fission yields required that both the krypton and 
xenon results fit the curve and that all the yields for each hump add 
Other than for the fine structure and the agreement with 
the curve very good. The absolute fission yields used Fig. give this 
excellent agreement are listed Table 


TABLE 
ABSOLUTE FISSION YIELDS OF KRYPTON AND XENON OCCURRING IN THE FAST NEUTRON FISSION 
Fission yield (%) Isotope Fission yield (%) 
1.99+0.01 


*Normalized the radiochemical curve Turkevich and Niday (1951), 
give the best fit with all krypton and xenon values. 


Neutron Multiplicity 
The curve shown Fig. symmetric about mass number 115.5 suggesting 


that the neutron multiplicity approximately for The most probable 
masses are and 140 also found Turkevich and Niday (1951). 


Fission Yields from Flux Calculations 

calculate the number fissions that have occurred the irradiated 
necessary know the fission cross section and the fission spectrum 
flux. The fission cross section function neutron energy and 


convenient define 


el7 


(1) 


a 


el? 


el? 


0 0 
where the number neutrons having energy between and 
and the fission cross section for neutrons energy Because the 
fission threshold 1.10 Mev. (Haxby 1940) the range integration 
for the numerator 1.10 Mev. The function for fission spectrum 
neutrons has been determined experimentally Watt (1952). The value 
for has been given Hughes and Harvey (1955). Graphical 
evaluation Equation gave value 0.075 barns for 

Fig. shows the geometric configuration the transformer. assuming 
sphere cross section the origin possible calculate the number 
fission neutrons which pass through per second. The number fast 
neutrons which pass through the sphere per second from annular uranium 
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where the number fission neutrons formed per per second, 
the number uranium atoms per and the absorption cross section 


Fic. The geometric configuration the transformer. 


natural uranium. The quantity the distance the fission neutrons travel 
through the fuel element wall and given 


(3) 


The total flux then 


expansion the exponential, equation (4) can integrated with the 
result for the fast flux 
0.472w. 


The quantity slowly varying function but with the boundary 
conditions for this case can treated constant Equation (4). 
evaluate the average thermal flux within the transformer must known. 
Calculations Case, Hoffmann, and Placzek (1953) for self-shielding 
corrections indicate that the average thermal flux within the transformer, 
Is 


where the unperturbed flux that position. megawatts power 
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for the transformer lattice position. Using 
the above relations the fast flux found 


about 21% the peak unperturbed thermal flux the same position. 
The contribution from the solid fuel element can shown negligible. 

From the spiking results given Table the number fissions per gram 
can calculated for irradiation and This has been done and 
the results are given Table VII, column where the yield per 


TABLE VII 
THE CALCULATED AND FITTED YIELD THE MASS CHAIN FOR FAST FISSION 


No. fissions/g. Yield mass chain (%) 
Average value.of (%) 6.00+0.06 


cent the mass chain. The number fissions per gram can also 
determined use the relation 
(5) fissions/g. (no. atoms/g.), 


where the irradiation time megawatts power. The value Equation 
(5) for and given Table VII, column Equating columns and 
value for can found. The value found this way for irradiation 
Aand given column The value for obtained from the 
curve given column 

The agreement between the two methods remarkably good and gives 
support both the normalization and the calculated fission spectrum flux. 


General Characteristics Fine Structure the Mass—Yield Curve 

Fig. shows the krypton yields for (Fleming, Tomlinson, and 
Thode 1954), (Blades, Fleming, and Thode 1956), and (Wanless and 
Thode 1955) fission. apparent that the curve this region 
moves higher yields the mass the fissile nuclide decreases. This in- 
dicates that the light mass hump actually displaced. The mass displacement 
these curves for fixed fission yield approximately equal the mass 
difference between the fissile nuclides. This implies that the heavy mass hump 
remains relatively constant position. Closer examination Fig. shows 
the well-known fine structure the and mass chains. The high yield 
all four cases probably arises because preference for this chain the 
original fission process cross-chain branching from the mass chain. The 
behavior the yield depends the fissile nuclide. This yield high for 
normal for and low for and The deviation the yield 
from the smooth curve can correlated with the value 
the fissile nuclide. This suggests that branching may occur 
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FISSION YIELD (%) 


MASS NUMBER 


Fic. The krypton yields for neutron-induced fission (Fleming al. 1954), 
(Blades al. 1956), and (Wanless and Thode 1955). 


the chain since the primary yield such prompt neutron emitter would 
related the Z/A value the fissile nuclide. The nuclide question must 
have fairly high primary yield and possess low threshold for neutron 
emission. These conditions are fulfilled This would indicate that the 
fine structure the mass region can explained the Glendenin (1949) 
boil-off effect since has neutrons. Thus the fine structure most 
similar that since both have about the same value, while the 
magnitude the krypton yields the highest since 232 the smallest 
mass number examined date. 

Fig..4 shows the 131, 132, 134, and 136 yields for (Fleming, 
Tomlinson, and Thode 1954), (Blades, Fleming, and Thode 1956), and 
(Wanless and Thode 1955) fission. The yields are lower than for 
the other cases and the 136/131 ratio larger indicating steeper 
yield curve. Both these facts suggest that the heavy mass hump may 
displaced higher masses than for the other three cases. The pattern, 
although having lower yields and, therefore, different region the 
yield curve, still shows the same general shape the other cases. The 
curve has satellite the side the main hump. This increased 
fission probability appears combination structural preference for 
the neutron shell the initial fission process and 
effect. 
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FISSION YIELD (%) 


MASS NUMBER 


Fic. The xenon yields for neutron-induced fission (Fleming al. 1954), 
(Blades al. 1956), and (Wanless and Thode 1955). 
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NOTES 


MULTIPLE SCATTERING MEASUREMENTS AND MESONS 
STOPPING NUCLEAR EMULSIONS* 


JUDEK 


Some previous work the author (Judek 1957) was designed test the 
accuracy the Constant Sagitta multiple scattering method (Dilworth, 
Goldsack, and Hirschberg 1954) for mass determination particles coming 
rest the emulsion. Particularly long track lengths were used, and measure- 
ments made for protons, deuterons, and tritons, indicating that the expected 
theoretical accuracy was attainable. Measurements have now been made 
long tracks and mesons stopping the emulsion. Tracks mesons 
and mesons their electron decays. They were selected from 
batch Ilford plates, thick, exposed cosmic radiation 
balloon and from stack stripped emulsions, 400 thick, 
exposed secondary radiation from carbon target bombarded the Bev. 
proton beam the Brookhaven 

The scattering scheme the Géttingen tables (Fay, Gottstein, and Hain 
1954) giving the mean sagitta, 0.5 for mesons was used, with the 
simplification that the cell length was varied steps instead (It 
was estimated that the error introduced this approximation (Judek 1957, 
Appendix was less than 0.1% for the scheme and therefore negligible.) 
The measured range each track was about giving about 250 inde- 
pendent cells, and was wholly contained one emulsion layer. Measurements 
were made Koristka M.S.2 microscope. The noise level correction was 
determined for nine meson tracks the method elimination using 
measurements full and half cell lengths (see Judek 1957, Appendix 
and this gave noise level 0.13+0.01 account the rather large 
scattering for and mesons, especially near the end the range, was 
necessary realign the tracks frequently. The effect imperfect alignment 
and also steepness tracks the value was carefully estimated. 
was found that under the experimental conditions, i.e. with the average 
inclination tracks, relative the direction the stage movement, 
3/100 radians figure obtained from the mean value the ratio: sideways 
displacement the track/cell length), the error due alignment (Judek 
1957, Appendix was only about 0.2% when averaged over the whole 
length track. The correction for steepness (taken proportional 
cos (Dilworth, Goldsack, and Hirschberg 1954), where the angle 
dip), however, was found the order for the last 1000 from 

are indebted the U.S. Office Naval Research for the cosmic ray exposure 


balloon flight, and the Brookhaven National Laboratory for the cosmotron 
exposure. 
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the end the range, and averaging about for the whole length 
the tracks. This estimate based separate measurements the angles 
dip successive small (100 portions along the whole lengths five 
meson tracks selected random. 

Fig. shows the distributions the values the mean sagittae (corrected 
for the noise level only), for individual tracks mesons and five 


Tr- mesons 


Number of tracks 


Fic. 


mesons. The dotted curves represent the Gaussian distributions centered 
the mean values 0.514 and 0.578 for the and meson groups respec- 
tively and with widths corresponding standard deviations calculated 
the basis the formula number independent cells per 
track). The experimental values the standard deviations, calculated from 
the individual values were found and for 
the and meson groups respectively, which good agreement with the 
above theoretical value 5%. Only partial resolution between the and 
mesons has been obtained the basis measurements cm. residual 
ranges, which expected from statistical considerations. 

The mean value 0.514+0.008 for mesons about higher 
than the expected value 0.500 correction for steepness tracks would 
reduce the experimental value probably about indicated measure- 
ments described above, thus giving good agreement between the theoretical 
and experimental values. The value the mass ratio mesons deduced 
the usual way from their respective average values (see Judek 1957, 
Appendix IV) was found 1.30+0.07, which again good agreement 
with the expected value 273/207 1.32. 

Fig. shows the variation the values with residual range deduced 
from observations all and meson tracks. order combine the 
measurements and mesons all the values were normalized 
dividing them the respective average values i.e. Do, for each group. 


‘ 
mesons 
0-5 
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Residual range (mm.) 


Fic. 


The dotted line the best line through the points and has 
slope less than 0.5% change per cm. residual range. can 
seen that there significant variation with residual range. 


C., J., and 1954. Nuovo cimento, 11, 113. 
Fay, H., K., and 1954. Nuovo cimento, Suppl. 11, 
1957. Nuovo cimento, 886. 
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LETTERS THE EDITOR 


Under this heading brief reports discoveries These 
reports should not exceed 600 words and, for any issue, should not liter thin six 
the first day the month issue. proof will sent the authors. 


Determination the Hydrogen Positions X-Ray Diffraction 


The structures and are isomorphous and the atomic coordinates, 
except those the hydrogen positions, have been known for many years (Aminoff 1919). Even 
for the hydrogen atoms, plausible locations had been proposed Bernal and Megaw (1935) 
the basis symmetry and electrostatic considerations. However, the infrared absorption 
spectra (Mara and Sutherland 1953) and (Petch and Megaw 1954) were 
found show unexpected complexities the region 2.3 3.5 Mara and Sutherland 
concluded that the unit cell must much larger than that found from X-ray 
analysis and that the positions the hydrogen atoms, proposed Bernal and Megaw, are 
incorrect. Their conclusions, correct, should apply well. The unit cell dimen- 
sions and space groups and accepted prior the infrared work, were 
confirmed (Petch and Megaw 1954) X-ray work, and strong evidence supporting the hydro- 
gen positions proposed Bernal and Megaw was preliminary report 

1955) further X-ray study The final results this work are reported ere. 

The specimen, cylinder with the dimensions 0.14 mm. 0.8 mm. and its axis parallel 
a-axis, was ground from crystal prepared the hydration larnite 
Professor Tilley, Department Mineralogy and Petrology, Cambridge University. The 
relative integrated intensities (Okil) reflections were measured with Geiger-counter 
X-ray diffractometer using radiation. The counter readings were corrected for dead- 
time losses and for background counts. The balanced filter technique (Ross 1928) was used 
eliminate the effect all white radiation except that within the wave-length range 
(the wave-lengths the and absorption edges), which was not 
and required separate correction. Absorption effects were small comparison with 
other possible errors neglected this work. The usual corrections were made for the 
Lorentz and polarization effects. The observed structure factors were put absolute 
basis scaling against the calculated values the final stages the analysis. 

The technique the F.)-synthesis, where the F,’s and the are the observed and 
the calculated structure factors used locate the hydrogen atoms. The 
F.)-synthesis was used because greatly reduces the effect Fourier series termination 
errors and shows the contributions hydrogen atoms very clearly when the are evaluated 
the basis model structure not containing the hydrogen atoms. For the first syntheses 
the Hartree scattering factors for and thermal corrections, were used 
the calculations the but the final stages refinement the Hartree values for the 
scattering factors were abandoned favor semiempirical values. The use these semi- 
empirical scattering factors, corrected allow for the anisotropic thermal vibrations the 
calcium and oxygen ions, resulted much better agreement between the and the 

projection the (2110) plane the electron density the asymmetric unit 
from which the calcium and oxygen ions have been subtracted shown Fig. This pro- 
jection the electron density, which was obtained with the final two-dimensional 
synthesis which the contributions the hydrogen atoms the were excluded, has 
its most prominent feature large peak superposed undulating background. For 
number reasons, believed that this peak represents electron associated with hydrogen 
atom and cannot due accumulation errors. Firstly, the electron density 0.9 
this peak times larger than the standard deviation the electron density. Secondly, 
shown Fig. the peak disappears completely leaving the background unchanged when 
the hydrogen atoms are included the model structure (on which the are based) the 
positions suggested the peak Fig. The probability accumulation errors pro- 
ducing not only such large peak, but also peak with shape identical that hydrogen 
atom projection, seems vanishingly small. Thirdly, the value the reliability factor, 
|F.|| decreases from 0.023 0.011 when the hydrogen atoms are 
also included the model structure the positions suggested this peak. Taking the 
position this peak the electron density the hydrogen position the O—H bond parallel 
the c-axis, originally suggested Bernal and Megaw, with length The 
final atomic coordinates (in special positions P3m1) obtained this work are 


listed below: 
(0, 


with 0.2330+0.0004 and 0.395+0.008. 
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cos 30° 


FIG.1 


@ cos 30° 


Fic. 1. A projection on the (2T10) plane of the electron density in the asymmetric unit of Ca(OH): from which 
the calcium and oxygen ions have been subtracted. The contours are drawn at every 0.le. A 2 with the zero 


contour omitted and the negative contours dotted. 
Fic. 2. As in Fig. 1 except that the hydrogen atom has also been subtracted. 


Busing and Levy (1957) have recently located the protons neutron diffraction 
and they report that the O—H bond has direction reported here but length 
projection along the c-axis Because electron-bonding effects there 
reason expect that the point maximum density for electron associated with 
hydrogen atom crystal must coincide with the proton position and this case, 
apparent that there considerable shift. This emphasizes the importance stating explicitly 
what one means atom position” discussions this subject. 

should added that Mara and Sutherland (1956), the basis study the infrared 
absorption spectrum the temperature liquid nitrogen, have concluded that 
possible explain the spectrum (in principle least) without invoking larger unit 
cell changing the equilibrium positions the hydrogen atoms from those proposed originally 
Bernal and Megaw. 

complete account this work will published Acta Crystallographica. 
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The Effect Incident Atomic Velocity the Structure Evaporated Silver Films 


recent letter Holland (1957) has suggested certain objections our evidence (Aziz 
and Scott 1956) for the effect incident atomic velocity the structure evaporated 
silver films. The main point which Holland raises that the latent heat condensation 
silver atoms large compared with the kinetic energy associated with their incident velocity 
and when this energy released will mask any differences due incident velocity. How- 
ever, the structure evaporated films depends the mobility the initial, individual 
atoms the substrate since this mobility determines the surface concentration nuclei 
which subsequent condensation occurs. The latent heat condensation applying 
does silver atoms condensing most favorable sites crystallite bulk silver not 
relevant the question single silver atom moving over glass plastic substrate. There 
what might considered latent heat for condensation the vapor atoms two- 
dimensional gas the substrate but small value compared with Moreover, the 
latent heat heat adsorption then shared among several the adjacent atoms 
the substrate. The adsorbed singlet atoms are motion the surface and collide with one 
another form immobile Frenkel doublets (Frenkel 1924 and Cockroft 1928) with release 
energy These doublets act centers condensation for other atoms. The sum 
i.e., the energy required remove one member the metallic doublet from the surface, 
was found many cases and particular for cadmium glass (Appleyard 1937) 
order magnitude smaller than are concerned, however, only with for the 
energy released when the atom becomes immobile will have bearing the concentration 
nuclei and hence the degree continuity the film structure. then probably only 
few per cent and therefore will not obscure the effect the differences the kinetic 
energy the atoms. 

Holland also mentions two unknown factors which may have influenced our results. Certainly 
there are unknown factors the formation evaporated deposits. However, the two factors 
referred to—substrate temperature and residual gas impurities— are not likely important 
our experiments. From studies the influence substrate temperature film structure 
and from direct measurements the substrate temperature under the conditions used, 
clear that the small temperature differences (5° which may exist could not 
responsible for the observed effects film structure. The influence the impurities the 
gas more difficult assess but since liquid air trap was employed the evaporating 
chamber very unlikely that any significant quantity water vapor would present 
during the evaporations. 
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